DFTT 10/03 



Random matrix theory and symmetric spaces 

O 
O 
(N 

^ ■ M. Caselle^ and U. Magnea^ 

OO ■ ^Department of Theoretical Physics, University of Torino 

^ and INFN, Sez. di Torino 

^ ' Via P. Giuria 1, 1-10125 Torino, Italy 

OH' '^Department of Mathematics, University of Torino 

Via Carlo Alberto 10, 1-10123 Torino, Italy 
^ ' caselle@to.infn.it 

o ■ 

' magnea@dm.unito.it 

o ■ 
B 

pH ■ Abstract 

o ■ 

In this review we discuss the relationship between random matrix theories and symmetric 
i. spaces. We show that the integration manifolds of random matrix theories, the eigenvalue 

' distribution, and the Dyson and boundary indices characterizing the ensembles are in strict 
correspondence with symmetric spaces and the intrinsic characteristics of their restricted 
root lattices. Several important results can be obtained from this identification. In par- 
ticular the Cartan classification of triplets of symmetric spaces with positive, zero and 
negative curvature gives rise to a new classification of random matrix ensembles. The re- 
view is organized into two main parts. In Part I the theory of symmetric spaces is reviewed 
with particular emphasis on the ideas relevant for appreciating the correspondence with 
random matrix theories. In Part II we discuss various applications of symmetric spaces 
to random matrix theories and in particular the new classification of disordered systems 
derived from the classification of symmetric spaces. We also review how the mapping from 
integrable Calogero-Sutherland models to symmetric spaces can be used in the theory of 
random matrices, with particular consequences for quantum transport problems. We con- 
clude indicating some interesting new directions of research based on these identifications. 
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1 Introduction 



The study of symmetric spaces has recently attracted interest in yarious branches of physics, 
ranging from condensed matter physics to lattice QCD. This is mainly due to the gradual 
understanding during the past few years of the deep connection between random matrix 
theories and symmetric spaces. Indeed, this connection is a rather old intuition, which 
traces back to Dyson P and has subsequently been pursued by seyeral authors, notably by 
Hiiffmann j2] . Recently it has led to seyeral interesting results, like for instance a tentatiye 
classification of the uniyersality classes of disordered systems. The latter topic is the main 
subject of this reyiew. 



The connection between random matrix theories and symmetric spaces is obtained simply 
through the coset spaces defining the symmetry classes of the random matrix ensembles. 
Although Dyson was the first to recognize that these coset spaces are symmetric spaces, 
the subsequent emergence of new random matrix symmetry classes and their classification 
in terms of Cartan's symmetric spaces is relatiyely recent [HI H El El IZ| • Since symmetric 
spaces are rather well understood mathematical objects, the main outcome of such an 
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identification is that several non-trivial results concerning the behavior of the random 
matrix models, as well as the physical systems that these models are expected to describe, 
can be obtained. 

In this context an important tool, that will be discussed in the following, is a class of inte- 
grable models named Calogero-Sutherland models [S]. In the early eighties, Olshanetsky 
and Perelomov showed that also these models are in one-to-one correspondence with sym- 
metric spaces through the reduced root systems of the latter 0. Thanks to this chain of 
identifications (random matrix ensemble - symmetric space - Calogero-Sutherland model) 
several of the results obtained in the last twenty years within the framework of Calogero- 
Sutherland models can also be applied to random matrix theories. 

The aim of this review is to allow the reader to follow this chain of correspondences. To 
this end we will devote the first half of the paper (sections |21 through [7j) to the necessary 
mathematical background and the second part f sections ISl through I1U|) to the applications 
in random matrix theory. In particular, in the last section we discuss some open directions 
of research. The reader who is not interested in the mathematical background could skip 
the first part and go directly to the later sections where we list and discuss the main results. 

This review is organized as follows: 

The first five sections of Part I (sections |2Hni) are devoted to an elementary introduction 
to symmetric spaces. As mentioned in the Abstract, these sections consist of the material 
presented in p^, which is a self-contained introductory review of symmetric spaces from 
a mathematical point of view. The material on symmetric spaces should be accessible to 
physicists with only elementary background in the theory of Lie groups. We have included 
quite a few examples to illustrate all aspects of the material. In the last section of Part I, 
section [71 we briefly introduce the Calogero-Sutherland models with particular emphasis 
on their connection with symmetric spaces. 

After this introductory material we then move on in Part II to random matrix theories 
and their connection with symmetric spaces (section |H)). Let us stress that this paper 
is not intended as an introduction to random matrix theory, for which very good and 
thorough references already exist [^|2Hll2Zll2Hll2ni- In this review we will assume that the 
reader is already acquainted with the topic, and we will only recall some basic information 
(definitions of the various ensembles, main properties, and main physical applications). 
The main goal of this section is instead to discuss the identifications that give rise to the 
close relationship between random matrix ensembles and symmetric spaces. 

Section M is devoted to a discussion of some of the consequences of the above mentioned 
identifications. In particular we will deduce, starting from the Cartan classification of 
symmetric spaces, the analogous classification of random matrix ensembles. We discuss 
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the symmetries of the ensembles in terms of the underlying restricted root system, and 
see how the orthogonal polynomials belonging to a certain ensemble are determined by 
the root multiplicities. In this section we also give some examples of how the connection 
between random matrix ensembles on the one hand, and symmetric spaces and Calogero- 
Sutherland models on the other hand, can be used to obtain new results in the theoretical 
description of physical systems, more precisely in the theory of quantum transport. 

The last section of the paper is devoted to some new results that show that the mathe- 
matical tools discussed in this paper (or suitable generalizations of these) can be useful for 
going beyond the symmetric space paradigm, and to explore some new connections between 
random matrix theory, group theory, and differential geometry. Here we discuss clustered 
solutions of the Dorokhov-Mello-Pereyra-Kumar equation, and then we go on to discuss 
the most general Calogero-Sutherland potential, given by the Weierstrass P-function, and 
show that it covers the three cases of symmetric spaces of positive, zero and negative cur- 
vature. Finally, in the appendix we discuss some intriguing exact results for the so called 
zonal spherical functions, which not only play an important role in our discussion, but are 
also of great relevance in several other branches of physics. 

There are some important and interesting topics that we will not review because of lack of 
space and competence. For these we refer the reader to the existing literature. In particular 
we shall not discuss: 

• the supersymmetric approach to random matrix theories and in particular their clas- 
sification in terms of supersymmetric spaces. Here we refer the reader to the original 
paper by M. Zirnbauer while a good introduction to the use of supersymmetry 
in random matrix theory and a complete set of the relevant references can be found 
in [121; 

• the very interesting topic of phase transitions. For this we refer to the recent and 
thorough review by G. Cicuta [T^ : 

• the extension to two-dimensional models of the classification of symmetric spaces, 
and more generally the methods of symmetric space analysis |14j : 

• the generalization of the classification of symmetric spaces to non-hermitean random 
matrices PH] (see however a discussion in the concluding section ITT|l: 

• the so called q-ensembles [TB] : 

• the two-matrix models J7j and multi-matrix models ^Hl and their continuum limit 
generalization. 
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The last item in the hst given above is a very interesting topic, which has several physical 
applications and would indeed deserve a separate review. The common feature of these 
two- and multi-matrix models which is of relevance for the present review, is that they all 
can be mapped onto suitably chosen Calogero-Sutherland systems. These models represent 
a natural link to two classes of matrix theories which are of great importance in high energy 
physics: on the one hand, the matrix models describing two-dimensional quantum gravity 
(possibly coupled to matter) and on the other hand, the matrix models pertaining 
to large QCD, which trace back to the original seminal works of 't Hooft 1201 ■ 
particular, a direct and explicit connection exists between multi-matrix models (the so 
called Kazakov-Migdal models) for large QCD [21] and the exactly solvable models of 
two-dimensional QCD on the lattice [22] • 

The mapping of these models to Calogero-Sutherland systems of the type discussed in this 
review can be found for instance in [23|. The relevance of these models, and in particular of 
their Calogero-Sutherland mappings, for the condensed matter systems like those discussed 
in the second part of this review, was first discussed in |24j. A recent review on this aspect, 
and more generally on the use of Calogero-Sutherland models for low-dimensional models, 
can be found in |25j . 

We will necessarily be rather sketchy in discussing the many important physical applica- 
tions of the random matrix ensembles to be described in section |H1 We refer the reader 
to some excellent reviews that have appeared in the literature during the last few years: 
the review by Beenakker [20] for the solid state physics applications, the review by Ver- 
baarschot 27J for QCD-related applications, and [2H1 120] for extensive reviews including a 
historical outline. 
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Part I 



The theory of symmetric spaces has a long history in mathematics. In this first part of the 
paper we will introduce the reader to some of the most fundamental concepts in the theory 
of symmetric spaces. We have tried to keep the discussion as simple as possible without 
assuming any previous familiarity of the reader with symmetric spaces. The review should 
be particularly accessible to physicists. In the hope of addressing a wider audience, we 
have almost completely avoided using concepts from differential geometry, and we have 
presented the subject mostly from an algebraic point of view. In addition we have inserted 
a large number of simple examples in the text, that will hopefully help the reader visualize 
the ideas. 

Since our aim in Part II will be to introduce the reader to the application of symmetric 
spaces in physical integrable systems and random matrix models, we have chosen the 
background material presented here with this in mind. Therefore we have put emphasis 
not only on fundamental issues but on subjects that will be relevant in these applications 
as well. Our treatment will be somewhat rigorous; however, we skip proofs that can be 
found in the mathematical literature and concentrate on simple examples that illustrate 
the concepts presented. The reader is referred to Helgason's book [20] for a rigorous 
treatment; however, this book may not be immediately accessible to physicists. For the 
reader with little background in differential geometry we recommend the book by Gilmore 
pi] (especially Chapter 9) for an introduction to symmetric spaces of exceptional clarity. 

In section |2l after reviewing the basics about Lie groups, we will present some of the most 
important properties of root systems. In section IHl we define symmetric spaces and discuss 
their main characteristics, defining involutive automorphisms, spherical decomposition of 
the group elements, and the metric on the Lie algebra. We also discuss the algebraic 
structure of the coset space. 

In section 0] we show how to obtain all the real forms of a complex semisimple Lie algebra. 
The same techniques will then be used to classify the real forms of symmetric spaces in 
section El In this section we also define the curvature of a symmetric space, and discuss 
triplets of symmetric spaces with positive, zero and negative curvature, all corresponding 
to the same symmetric subgroup. We will see why curved symmetric spaces arise from 
semisimple groups, whereas the fiat spaces are associated to non-semisimple groups. In 
addition, in section^lwe will define restricted root systems. The restricted root systems are 
associated to symmetric spaces, just like ordinary root systems are associated to groups. 
As we will discuss in detail in Part II of this paper, they are key objects when considering 
the integrability of Calogero-Sutherland models. 

In section E] we discuss Casimir and Laplace operators on symmetric spaces and men- 
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tion some known properties of the eigenf unctions of the latter, so called zonal spherical 
functions. These functions play a prominent role in many physical applications. 

The introduction to symmetric spaces we present contains the basis for understanding the 
developments to be discussed in more detail in Part II. The reader already familiar with 
symmetric spaces is invited to start reading in the last section of Part I, section [71 where 
we give a brief introduction to Calogero-Sutherland models. 



2 Lie groups and root spaces 

In this introductory section we define the basic concepts relating to Lie groups. We will 
build on the material presented here when we discuss symmetric spaces in the next section. 
The reader with a solid background in group theory may want to skip most or all of this 
section. 



2.1 Lie groups and manifolds 

A manifold can be thought of as the generalization of a surface, but we do not in general 
consider it as embedded in a higher-dimensional euclidean space. A short introduction to 
differentiable manifolds can be found in ref. and a more elaborate one in refs. [SBj and 
[SI] (Ch. III). The points of an A^-dimensional manifold can be labelled by real coordinates 
{x^, ...,x^). Suppose that we take an open set Ua of this manifold, and we introduce local 
real coordinates on it. Let ipa be the function that attaches N real coordinates to each 
point in the open set Ua- Suppose now that the manifold is covered by overlapping open 
sets, with local coordinates attached to each of them. If for each pair of open sets Ua, Up, 
the fuction ifja o ?/'^^ is differentiable in the overlap region Ua^iU^, it means that we can 
go smoothly from one coordinate system to another in this region. Then the manifold is 
differentiable. 

Consider a group G acting on a space V. We can think of G as being represented by 
matrices, and of \^ as a space of vectors on which these matrices act. A group element 
g & G transforms the vector v & V into gv = v' . 

If G is a Lie group, it is also a differentiable manifold. The fact that a Lie group is a 
differentiable manifold means that for two group elements g, g' G G, the product {g,g') G 
G X G ^ gg' G G and the inverse g g~^ are smooth (C7°°) mappings, that is, these 
mappings have continuous derivatives of all orders. 
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Example: The space R" is a smooth manifold and at the same time an abelian group. 
The "product" of elements is addition (x, x') x + x' and the inverse of x is —x. These 
operations are smooth. 

Example: The set GL{n, R) of nonsingular real n x n matrices M, detM 7^ 0, with 
matrix multiplication (M, A^) — > MN and multiplicative matrix inverse M —>■ is a 
non-abelian group manifold. Any such matrix can be represented as M = e^**'"^* where 
Xi are generators of the GL(n, R) algebra and are real parameters. 



2.2 The tangent space 

In each point of a different iable manifold, we can define the tangent space. If a curve 
through a point P in the manifold is parametrized by t e R 



where P = (a;^(0), x^(0)), then A = (A^, A^) = {x^{0), x'^{0)) is a tangent vector 
at P. Here i;°(0) = ■^x'^{t)\t=o. The space spanned by all tangent vectors at P is the 
tangent space. In particular, the tangent vectors to the coordinate curves (the curves 
obtained by keeping all the coordinates fixed except one) through P are called the natural 
basis for the tangent space. 

Example: In euclidean 3-space the natural basis is {ex,ey,ez}- On a patch of the unit 
2-sphere parametrized by polar coordinates it is {e6),e^}. 

For a Lie group, the tangent space at the origin is spanned by the generators, that play the 
role of (contravariant) vector fields (also called derivations), expressed in local coordinates 
on the group manifold as X = X°'{x)da (for an introduction to differential geometry see 
ref. Ch. 5, or jMj)- Here the partial derivatives da = form a basis for the vector 
field. That the generators span the tangent space at the origin can easily be seen from the 
exponential map. Suppose X is a generator of a Lie group. The exponential map then 
maps X onto e*^, where t is a parameter. This mapping is a one-parameter subgroup, 
and it defines a curve x{t) in the group manifold. The tangent vector of this curve at the 
origin is then 



x"(t) = a;'^(0) + \H 



a = 1, ...,X 



(2.1) 




t=o — 



X 



(2.2) 
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All the generators together span the tangent space at the origin (also called the identity 
element). 



2.3 Coset spaces 

The isotropy subgroup G„p of a group G at the point f o G ^ is the subset of group elements 
that leave vq fixed. The set of points that can be reached by applying elements g & G to 
vq is the orbit of G at fo, denoted Gvq. If Gvq = V for one point vq, then this is true for 
every v ^ V. We then say that G acts transitively on V. 

In general, a symmetric space can be represented as a coset space. Suppose if is a subgroup 
of a Lie group G. The coset space G/H is the set of subsets of G of the form gH, for g E G. 
G acts on this coset space: gi{gH) is the coset {gig)H. We will refer to the elements of the 
coset space by g instead of by gH, when the subgroup H is understood from the context, 
because of the natural mapping described in the next paragraph. If g ^ H, gH corresponds 
to a point on the manifold G/H away from the origin, whereas hH = H [h E H) is the 
identity element identified with the origin of the symmetric space. This point is the north 
pole in the example below. 

If G acts transitively on V, then V = Gv for any v E V. Since the isotropy subgroup 
Gvq leaves a fixed point vq invariant, gGv^vo = gvo = v E V, we see that the action of the 
group G on V defines a bijective action of elements of G/Gyg on V. Therefore the space 
V on which G acts transitively, can be identified with G/G^g, since there is one-to-one 
correspondence between the elements of V and the elements of G/Gy^. There is a natural 
mapping from the group element g onto the point gvQ on the manifold. 

Example: The SO (2) subgroup of SO (3) is the isotropy subgroup at the north pole of 
a unit 2-sphere imbedded in 3-dimensional space, since it keeps the north pole fixed. On 
the other hand, the north pole is mapped onto any point on the surface of the sphere by 
elements of the coset S0{3)/ S0(2). This can be seen from the explicit form of the coset 
representatives. As we will see in eq. fj3.2U|) in subsection 13.51 the general form of the 
elements of the coset is 



M = exp 



( 



C 
-C^ 



) 



( 



VI - x^x 



X 



) 



(2.3) 



where C is the matrix 
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(2.4) 



and t^, are real coordinates. I2 in eq. (j2.3p is the 2x2 unit matrix. For the coset space 
S0{?>)/S0{2), M is equal to 



M = exp(^^fL,J, Li 
The third S'0(3) generator 












( 












■ ■ 












\ 


-1 


/ 









/ 



(2.5) 



/ 1 
-10 

V 



(2.6) 



spans the algebra of the stability subgroup 50(2), that keeps the north pole fixed: 





exp(t3L3) I 





1 



(2.7) 



The generators Li {i = 1,2,3) satisfy the 5*0(3) commutation relations [Li,Lj] = \eijkLk. 
Note that since the Li and the are real, C*^ = C'^ . 

In ()2.3|1 . M is a general representative of the coset 50 (3) /SO (2). By expanding the 
exponential we see that the explicit form of M is 



M 



1 + {t^) 



(tl)2 + (t2)2 



,2 



(cosV(ti)2+(t2)2_i) .,n2 (cosv/(t^)^ + (f^)^-l) ^i sinV(ti)2+(t2) 

(tl)2 + (t2)2 -L "t" l*^ J (il)2_,.(t2)2 t _ /77T^iT775^ 

^i sin^(fi)2 + (f2) 



, COS J(tl)2 + (t2)^ 

^(tl)2 + (t2)2 \j\ ) ^ \ ) 



(2.8) 



Thus the matrix X 



X 

y 



is given in terms of the components of C by (cf. eq. 1)3.211) ^ 
12 



X 



X 

y 



l sinV(tl)2+(t2)2 



(2.9) 



Defining now z — COS'^ {t^y + (i^)^, we see that the variables x, y, z satisfy the equation 
of the 2-sphere: 



+ y2 + = 1 



(2.10) 



When the coset space representative M acts on the north pole it is easily seen that the 
orbit is all of the 2-sphere: 







/ . 


. x\ 




/ 


X 


\ 




M 


{']- 




■ y 


{']- 




y 




(2.11) 






\ • 


■ w 




\ 


z 


J 





This shows that there is one-to-one correspondence between the elements of the coset and 
the points of the 2-sphere. The coset SO{3)/SO{2) can therefore be identified with a unit 
2-sphere imbedded in 3-dimensional space. 



2.4 The Lie algebra and the adjoint representation 

A Lie algebra G is a vector space over a field F. Multiplication in the Lie algebra is given 
by the bracket [X, Y]. It has the following properties: 

[1] liX,Y e G, then [X, Y] e G, 

[2] [X, aY + pZ] = a[X, Y] + p[X, Z] for a, P e F, 

[3] [X,Y] = -[Y,X], 

[4] [X, [Y, Z]] + [Y, [Z, X]] + [Z, [X, Y]] = (the Jacobi identity). 

The algebra G generates a group through the exponential mapping. A general group 
element is 



M = exp (l^t'^i^ ; f eF, XieG 



(2.12) 
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We define a mapping adX from the Lie algebra to itself by adX : Y — > [-'^, 1^]. The 
mapping X adX is a representation of the Lie algebra called the adjoint representation. 
It is easy to check that it is an automorphism: it follows from the Jacobi identity that 
[adXj,adXj] = ad[Xi,Xj]. Suppose we choose a basis {Xi} for G. Then 



8^dXi{Xj)^[Xi,Xj]^C^.Xk 



(2.13) 



where we sum over k. The Cfj are called structure constants. Under a change of basis, they 
transform as mixed tensor components. They define the matrix {Mi)jk = Cff. associated 
with the adjoint representation of X^. One can show that there exists a basis for any 
complex semisimple algebra in which the structure constants are real. This means the 
adjoint representation is real. Note that the dimension of the adjoint representation is 
equal to the dimension of the group. 



Example: Let's construct the adjoint representation of SU{2). The generators in the 
defining representation are 



1/1 
2 I -1 



J4 




(2.14) 



and the commutation relations are 



[J3, J±] =±J±, [J+, J_] =2J3 



(2.15) 



The structure constants are therefore 6*3^ 



~^+3 — — ^-3 — ^1 - 

and the adjoint representation is given by {M^)^^ = 1, (M3) = —1, (M_|_)_|_3 

(M_,_)3_ = 2, (M_)_3 = 1, (M_)3_|_ = —2, and all other matrix elements equal to 0: 



-1, 



M3 



/ \ 
1 

\o -1) 



/ 2 \ 
-10 

V / 



/ -2 \ 


VI 0; 



(2.16) 



These representation matrices are real, have the same dimension as the group, and satisfy 
the SU(2) commutation relations [M3,M±] = ±M±, [M+,M_] = 2M3. 
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2.5 Semisimple algebras and root spaces 

In this paragraph we will briefly recall the basic facts about root spaces and the classifica- 
tion of complex simple Lie algebras, to set the stage for our discussion of real forms of Lie 
algebras and finally symmetric spaces. 

An ideal, or invariant subalgebra I is a subalgebra such that [G,I] C I. An abelian ideal 
also satisfies [I, I] = 0. A simple Lie algebra has no proper ideal. A semisimple Lie algebra 
is the direct sum of simple algebras, and has no proper abelian ideal (by proper we mean 
different from {0}). 

A Lie algebra is a linear vector space over a field F, with an antisymmetric product defined 
by the Lie bracket (cf. subsection 12. 4j) . If F is the field of real, complex or quaternion 
numbers, the Lie algebra is called a real, complex or quaternion algebra. A complexification 
of a real Lie algebra is obtained by taking linear combinations of its elements with complex 
coefficients. A real Lie algebra H is a real form of the complex algebra G if G is the 
complexification of H. 

In any simple algebra there are two kinds of generators: there is a maximal abelian subal- 
gebra, called the Cartan subalgebra Hq = {-f^i, -f^r}; [Hi, Hj\ = for any two elements 
of the Cartan subalgebra. There are also raising and lowering operators denoted E^. a 
is an r-dimensional vector a = (ai, ...,ar) and r is the rank of the algebra. ^ The lat- 
ter are eigenoperators of the Hi in the adjoint representation belonging to eigenvalue aj: 
[Hi,Ea] = OiEa- For each eigenvalue, or root ai, there is another eigenvalue — and a 
corresponding eigenoperator E^a under the action of Hi. 

Suppose we represent each element of the Lie algebra by an n x n matrix. Then [Hi, Hj] = 
means the matrices Hi can all be diagonalized simultaneously. Their eigenvalues fii are 
given by Hi\fj,) = fii\fi), where the eigenvectors are labelled by the weight vectors /i = 

(/ii, [sni- 

A weight whose first non-zero component is positive is called a positive weight. Also, a 
weight /i is greater than another weight /i' if /i — fi' is positive. Thus we can define the 
highest weight as the one which is greater than all the others. The highest weight is unique 
in any representation. 

The roots a, = a{Hi) of the algebra G are the weights of the adjoint representation. Recall 
that in the adjoint representation, the states on which the generators act are defined by 

^The rank of an algebra is defined through the secular equation (see subsection 16.111 . For a non- 
semisimple algebra, the maximal number of mutually commuting generators can be greater than the rank 
of the algebra. 
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the generators themselves, and the action is defined by 



Xa\X,) = adX,(X,) = [Xa^X,] (2.17) 
The roots are functionals on the Cartan subalgebra satisfying 

a.m{E^) = [H^, Eo] = a{H,)E^ (2.18) 

where Hi is in the Cartan subalgebra. The eigenvectors are called the root vectors. 
These are exactly the raising and lowering operators E±a for the weight vectors /x. There 
are canonical commutation relations defining the system of roots belonging to each simple 
rank r-algebra. These are summarized below: ^ 

[Hi, Hj\ = 0, [Hi, E^] = aiE^, [E^, = aiHi (2.19) 
One can prove the fundamental relation 

2a ■ fi 



-ip-q) (2.20) 



where a is a root, yU is a weight, and p, q are positive integers such that -Ea|yU + pa) = 0, 
— qa) = ^. This relation gives rise to the strict properties of root lattices, and 
permits the complete classification of all the complex (semi)simple algebras. 

^For the reader who wants to understand more about the origin of the structure of Lie algebras, we 
recommend Chapter 7 of Gilmore • 

■^Here the scalar product • can be defined in terms of the metric on the Lie algebra. For the adjoint 
representation, /j, is a root /3 and 

2a-/3 _ 2K{H^,Hp) _ 2/3(g„) 

where K denotes the Killing form (see paragraph ia. 4|l . There is always a unique element Ha in the algebra 
such that K{H, Ha) — a{H) for each H e Hq (see for example Ch. 10). In general for a linear form 
H on the Lie algebra, 

2a^ ^ (2.22) 

a'^ a{Ha) 

Then ^ is a highest weight for some representation if and only if this expression is an integer for each 
positive root a. 
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Eq. ()2.20|) is true for any representation, but has particularly strong implications for the 
adjoint representation. In this case /i is a root. As a consequence of eq. ()2.2Up . the possible 
angle between two root vectors of a simple Lie algebra is limited to a few values: these turn 
out to be multiples of | and | (see e.g. jHSj, Ch. VI). The root lattice is invariant under 
reflections in the hyperplanes orthogonal to the roots (the Weyl group). As we will shortly 
see, this is true not only for the root lattice, but for the weight lattice of any representation. 

Note that the roots a are real-valued linear functionals on the Cartan subalgebra. There- 
fore they are in the space dual to Hq. A subset of the positive roots span the root lattice. 
These are called simple roots. Obviously, since the roots are in the space dual to Hq, the 
number of simple roots is equal to the rank of the algebra. 

The same relation ()2.20|) determines the highest weights of all irreducible representations. 
Setting p = 0, choosing a positive integer g, and letting a run through the simple roots, 
a = a* (i = 1, ...,t), we find the highest weights /x* of all the irreducible representations 
corresponding to the given value of q jSH]- For example, for g = 1 we get the highest weights 
of the r fundamental representations of the group, each corresponding to a simple root a*. 
For higher values of q we get the highest weights of higher-dimensional representations of 
the same group. 

The set of all possible simple root systems are classified by means of Dynkin diagrams, 
each of which correspond to an equivalence class of isomorphic Lie algebras. The classical 
Lie algebras SU(n + 1, C), SO(2n + 1, C), Sp(2n, C) and SO(2n, C) correspond to root 
systems A^, 5„, C„, and D„, respectively. In addition there are five exceptional algebras 
corresponding to root systems E^^ £^7, iJg, F4 and G2- Each of these complex algebras in 
general has several real forms associated with it (see section^)). These real forms correspond 
to the same Dynkin diagram and root system as the complex algebra. Since we will not 
make reference to Dynkin diagrams in the following, we will not discuss them here. The 
interested reader can find sufficient material for example in the book by Georgi [5^ . 

The (semi)simple complex algebra G decomposes into a direct sum of root spaces ^Hj : 

G = Ho©^G„ (2.23) 

OL 

where G^ is generated by \E±a\. This will be evident in the example given below. 

Example: The root system An~\ corresponds to the complex Lie algebra SL(n, C) and all 
its real forms. In a later section we will see how to construct all the real forms associated 
with a given complex Lie algebra. Let's see here explicitly how to construct the root lattice 
of SU(3, C), which is one of the real forms of SL(3, C). 
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The generators are determined by the commutation relations. In physics it is common to 
write the commutation relations in the form 



(2.24) 



(an alternative form is to define the generators as Xj = iTi and write the commutation re- 
lations as [Xj,Xj] = —fijkXk) where fijk are structure constants for the algebra SU(3, C). 

Using the notation g = e**""^" for the group elements (with t°- real and a sum over a implied), 
the generators in the fundamental representation of this group are hermitean^: 




/ 1 


V 1 



To 



( ° 

i 












(2.25) 








— b I 

; 



1 

2^3 



1 
1 






0-2 



In high energy physics the matrices 2Ta are known as Gell-Mann matrices. The generators 
are normalized in such a way that tr(T(jTf,) = \5ah- Note that Ti, T2, T3 form an SU(2, C) 
subalgebra. We take the Cartan subalgebra to be Hq = {TajTg}. The rank of this group 
is r = 2. 



Let's first find the weight vectors of the fundamental representation. To this end we look 
for the eigenvalues Hi of the operators in the abelian subalgebra Hq: 



''Note that we have written an expHcit factor of i in front of the generators in the expression for the 
group elements. This is often done for compact groups; since the KiUing form fsubsection I3.4|l has to 
be negative definite, the coordinates of the algebra spanned by the generators must be purely imaginary. 
Here we use this notation because it is conventional. If we absorb the factor of i into the generators, we 
get antihermitean matrices Xa — iTa] we will do this in the example in subsection 13.11 to comply with 
eq. Of course, the matrices in the algebra are always antihermitean. 
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(2.26) 



therefore the eigenvector (10 0)^ corresponds to the state \n) where 



is distinguished by its eigenvalues under the operators Hi of the Cartan subalgebra. In 
the same way we find that (0 1 0)"^ and (0 1)"^ correspond to the states labelled by weight 
vectors 



respectively, fi, /i', and /i" are the weights of the fundamental representation p = D and 
they form an equilateral triangle in the plane. The highest weight of the representation D 

IS /i = 

There is also another fundamental representation D of the algebra SU(3, C), since it 
generates a group of rank 2. Indeed, from eq. ()2.2()|1 . for p = 0, g = 1, there is one highest 
weight /i*, and one fundamental representation, for each simple root a*. The highest weight 
/i of the representation D is 



The highest weights of the representations corresponding to any positive integer q can be 
obtained as soon as we know the simple roots. Then, by operating with lowering operators 
on this weight, we obtain other weights, on which we can further operate with lowering 
operators until we have obtained all the weights in the representation. For an example of 
this procedure see [22] , Ch. IX. 

Let's see now how to obtain the roots of SU(3, C). Each root vector corresponds to ei- 
ther a raising or a lowering operator: Ea is the eigenvector belonging to the root = a{Hi) 
under the adjoint representation of Hi, like in eq. ()2.H2|1 . Each raising or lowering operator 
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is a linear combination of generators Tj that takes one state of the fundamental represen- 
tation to another state of the same representation: E±a\f^) = N±a,^J.\^J' ± «)• Therefore the 
root vectors a will be differences of weight vectors in the fundamental representation. We 
find the raising and lowering operators E±a to be 

^±(i,o) = ^(Ti±zT2) 

^±(i,^) = 72(^4±^T5) (2.30) 

These generate the subspaces Gq, in eq. ()2.23|) . In the fundamental representation, we find 
using the Gell-Mann matrices that these are matrices with only one non-zero element. For 
example, the raising operator Ea that corresponds to the root a = (1, 0) is 



-^+(1,0) = — 




f2.311 



This operator takes us from the state = | — |, to the state = ||, ^^)- The 
components of the root vectors of SU(3, C) are the eigenvalues ai of these under the 
adjoint representation of the Cartan subalgebra. That is. 



Hi\Ea) = adHi{Ea) = [Hi, Ea] = ai|i?Q 



(2.32) 



This way we easily find the roots: we can either explicitly use the structure constants of 
SU{3) in [Ta,Tb] = ifabcTc = —iC^f,Tc (note the explicit factor of i due to our conven- 
tions regarding the generators) or we can use an explicit representation for Hi, Ea like in 
eqs. ()2.25|) . ()2.30|) . ()2.3H) . to calculate the commutators: 



adi^i(E±(i,o)) = [H,, E±(i,o)] = [Ts, ^(Ti ± iT^)] = -^{iT2 ± T,) = ±E±(i,o) = «f ^±(i,o) 
adi72(^±(i,o)) = [H2, ^±(1,0)] = [Tg, -j^in ± tT2)] = = a^^±(i,o) 

(2.33) 
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The root vector corresponding to the raising operator -E+{i,o) is thus a = {af, a^) = (1, 0) 
and the root vector corresponding to the lowering operator -^-(1,0) is ~« = (aj~,a2) = 
(—1,0). These root vectors are indeed the differences between the weight vectors /i = 
Q, and n' = (— |, of the fundamental representation. 

In the same way we find the other root vectors (^±|,±^^, (^=f|,±-^), and (0,0) (with 
multiplicity 2), by operating with Hi and H2 on the remaining E±aS and on the i^j's. The 
last root with multiplicity 2 has as its components the eigenvalues under Hi, H2 of the 
states \Hi) and 1-^2): Hi\Hj) = [Hi, Hj] = 0; i, j G {1, 2}. The root vectors form a regular 
hexagon in the plane. The positive roots are (1,0), = (^^, and = (^|, — The 
latter two are simple roots. (1, 0) is not simple because it is the sum of the other positive 
roots. There are two simple roots, since the rank of SU{3) is 2 and the root lattice is 
two-dimensional. 

The root lattice of SU{3) is invariant under reflections in the hyperplanes orthogonal to the 
root vectors. This is true of any weight or root lattice; the symmetry group of reflections in 
hyperplanes orthogonal to the roots is called the Weyl group. It is obtained from eq. ()2.20p : 
since for any root a and any weight /i, 2(a ■ fi)/a'^ is the integer q — p. 



^' = ^- ^h-tla (2.34) 
is also a weight. Eq. (|2.34j) is exactly the above mentioned reflection, as can easily be seen. 



2.6 The Weyl chambers 

The roots are linear functionals on the Cartan subalgebra. We may denote the Cartan 
subalgebra by Ho and its dual space by Hg. A Weyl reflection like the one in ()2.34|) can 
be defined not only for the weights or roots n in the space Hq, but for an arbitrary vector 
q e Hq or, in all generality, for a vector q in an arbitrary finite-dimensional vector space: 



Sa{q) = q — a*{q)a (2.35) 

Note that q G Hq is in the space dual to Hq and may denote a root. In ()2.35p the function 
a*{q) is a linear functional on Hq such that a* (a) = 2. We will be concerned only with 
the crystallographic case when a*{q) is integer. We denote the hyperplanes in Hq where 
the function a*{q) vanishes by iJ*^"^: 
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= {g G H* : a*{q) = 0} 



(2.36) 



^(") is orthogonal to the root a, and Sa{q) is a reflection in this hyperplane. 

By identifying the dual spaces Hq and Hg (this is possible since they have the same 
dimension), we can consider hyperplanes like the ones in ()2.36p in the space Hq. The role 
of the linear functional a* (q) is then played by 



2q-a 2q{H^) 
« W = — — = -jTrr 2.37 

where a{Ha) = K{Ha, Ha)- Here K is the Killing form (a metric on the algebra to be 
defined in paragraph 13. 4|) and Ha is the unique element in Hq such that K{H, Ha) = a{H). 
The open subsets of Hq where roots are nonzero are called Weyl chambers. Consequently, 
the walls of the Weyl chambers are the hyperplanes in Hq where the roots q{H) are zero. 



2.7 The simple root systems 



We have just shown by an example, in subsection 12.51 how to obtain a root system of 
type An- In general, for any simple algebra the commutation relations determine the 
Cartan subalgebra and raising and lowering operators, that in turn determine a unique 
root system, and correspond to a given Dynkin diagram. In this way we can classify all 
the simple algebras according to the type of root system it possesses. The root systems 
for the four infinite series of classical non-exceptional Lie groups can be characterized as 
follows jSH] (denote the r-dimensional space spanned by the roots by V and let {ei, ...e„} 
be a canonical basis in R"): 

An-i'. Let V be the hyperplane in R" that passes through the points (1, 0, 0, ...0), (0, 1, 0, 0), 

(0, 0, 0, 1) (the endpoints of the Cj, i = l,...,n). Then the root lattice contains the 
vectors {ci — ej,i ^ j}. 

Bn'- Let V be R"; then the roots are {icj, ±ej ± ej,i ^ j}. 
Cn- Let V be R"; then the roots are {±2ej, ±ej ± e^, i ^ j}. 
Dn- Let V be R"; then the roots are {icj ± Cj, i 7^ j}. 
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The root lattice BCn, that we will discuss in conjunction with restricted root systems, is 
the union of -B„ and C„. It is characterized as follows: 

BCn'- Let V be R"; then the roots are {±ej, ±2ej, ie^ ± ej,i 7^ j}. 

Because this system contains both and 2ej, it is called non-reduced (normally the only 
root collinear with a is —a). However, it is irreducible in the usual sense, which means it 
is not the direct sum of two disjoint root systems and C„. This can be seen from the 
root multiplicities (cf. Table IH). 

The semisimple algebras are direct sums of simple ones. That means the simple constituent 
algebras commute with each other, and the root systems are direct sums of the correspond- 
ing simple root systems. Therefore, knowing the properties of the simple Lie algebras, we 
also know the semisimple ones. 



3 Symmetric spaces 

In the previous section, we have reminded ourselves of some elementary facts concerning 
root spaces and the classification of the complex semisimple algebras. In this section we 
will define and discuss symmetric spaces. 

A symmetric space is associated to an involutive automorphism of a given Lie algebra. 
As we will see, several different involutive automorphisms can act on the same algebra. 
Therefore we normally have several different symmetric spaces deriving from the same Lie 
algebra. The involutive automorphism defines a symmetric subalgebra and a remaining 
complementary subspace of the algebra. Under general conditions, the complementary 
subspace is mapped onto a symmetric space through the exponential map. In the following 
subsections we make these statements more precise. We discuss how the elements of the Lie 
group can act as transformations on the elements of the symmetric space. This naturally 
leads to the definition of two coordinate systems on symmetric spaces: the spherical and 
the horospheric coordinate systems. The radial coordinates associated to each element of 
a symmetric space through its spherical or horospheric decomposition will be of relevance 
when we discuss the radial parts of differential operators on symmetric spaces in section 
ini In the same section we explain why these operators are important in applications to 
physical problems, and in Part II we will discuss some of their uses. 

In all of this paper we will distinguish between compact and non-compact symmetric 
spaces. In order to give a precise notion of compactness, we will define the metric tensor 
on a Lie algebra in terms of the Killing form in subsection 13.41 The latter is defined as a 
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symmetric bilinear trace form on the adjoint representation, and is therefore expressible 
in terms of the structure constants. We will give several examples of Killing forms later, 
as we discuss the various real forms of a Lie algebra. The metric tensor serves to define 
the curvature tensor on a symmetric space (subsection 15. HI . It is also needed in computing 
the Jacobian of the transformation to radial coordinates. This Jacobian is relevant in 
calculating the radial part of the Laplace-Beltrami operator (see paragraph 16. 2j) . 

We will close this section with a discussion of the general algebraic form of coset represen- 
tatives in subsection 13.51 



3.1 Involutive automorphisms 

An automorphism of a Lie algebra G is a mapping from G onto itself such that it preserves 
the algebraic operations on the Lie algebra. For example, if a is an automorphism, it 
preserves multiplication: [a{X),(y(Y)] = a{[X,Y]), for X,YeG. 

Suppose that the linear automorphism a : G ^ G is such that = 1, but a is not the 
identity. That means that a has eigenvalues ±1, and it splits the algebra G into orthogonal 
eigensubspaces corresponding to these eigenvalues. Such a mapping is called an involutive 
automorphism. 

Suppose now that G is a compact simple Lie algebra, a is an involutive automorphism of 
G, and G = K © P where 

a{X) =X for X G K, a{X) = -X for X G P (3.1) 

From the properties of automorphisms mentioned above, it is easy to see that K is a 
subalgebra, but P is not. In fact, the commutation relations 

[K,K]CK, [K,P]CP, [P,P]CK (3.2) 

hold. A subalgebra K satisfying ()3.2j] is called a symmetric subalgebra. If we now multiply 
the elements in P by z (the 'Weyl unitary trick"), we construct a new noncompact algebra 
G* = K © zP. This is called a Cartan decomposition., and K is a maximal compact 
subalgebra of G*. The coset spaces G/K and G* /K are symmetric spaces. 

Example: Suppose G = SU{n,C), the group of unitary complex matrices with determi- 
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nant +1. The algebra of this group then consists of complex antihermitean^ matrices of 
zero trace (this follows by differentiating the identities UW = 1 and detU = 1 with respect 
to t where U{t) is a curve passing through the identity at t = 0); a group element is written 
as g = e*"^" with real. Therefore any matrix X in the Lie algebra of this group can 
be written X = A + iB, where A is real, skew-symmetric, and traceless and B is real, 
symmetric and traceless. This means the algebra can be decomposed as G = K©P, where 
K is the compact connected subalgebra SO(n, R) consisting of real, skew-symmetric and 
traceless matrices, and P is the subspace of matrices of the form iB, where B is real, 
symmetric, and traceless. P is not a subalgebra. 

Referring to the example for SU(3, C) in subsection 12.51 we see, setting Xa = iTa, that 
the {Xa} split into two sets under the involutive automorphism a defined by complex 
conjugation a = K. This splits the compact algebra G into K © P, since P consists of 
imaginary matrices: 




1 
-1 |,i 




K={X2,X5,X7} = 

P = {Xi, X3, X4, Xg, Xg} 

/ 1 

I 



K spans the real subalgebra SO(3,R). Setting X2 
commutation relations for the subalgebra are = 
zHo = {X^,X^} is here entirely in the subspace P. 
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(3.3) 
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3) 



X5 = L2, X^ = Li, the 
The Cartan subalgebra 



Going back to the general case of G = SU(n, C), we obtain from G by the Weyl unitary 
trick the non-compact algebra G* = K © iP. iP is now the subspace of real, symmetric, 
and traceless matrices B. The Lie algebra G* = SL(n, R) is then the set of n x n real 
matrices of zero trace, and generates the linear group of transformations represented by 
real n x n matrices of unit determinant. 



The involutive automorphism that split the algebra G above was defined to be complex 
conjugation a = K. The involutive automorphism that splits G* is defined by cr{g) = 
{g^)~^ for g G G*, as we will now see. On the level of the algebra, (y{g) = (g^)^^ means 

^See the footnote in subsection 12. 51 
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5"(X) = — X^. Suppose now g = e*"^ G G* with X real and traceless and t a real parameter. 
If now X is an element of the subalgebra K, we then have 5"(X) = +X, i.e. —X^ = X 
and X is skew-symmetric. If instead X G iP, we have 5"(X) = — X^ = —X, i.e. X is 
symmetric. The decomposition G* = K © iP is the usual decomposition of a SL(n, R) 
matrix in symmetric and skew-symmetric parts. 

G/K = SU{n,G)/ SO{n, R) is a symmetric space of compact type, and the related sym- 
metric space of non-compact type is G*/K = SL{n, R)/ SO{n, R). 



3.2 The action of the group on the symmetric space 

Let G be a semisimple Lie group and K a compact symmetric subgroup. As we saw in 
the preceding paragraph, the coset spaces G/K and G*/K represent symmetric spaces. 
Just as we have defined a Cartan subalgebra and the rank of a Lie algebra, we can define, 
in an exactly analogous way, a Cartan subalgebra and the rank of a symmetric space. A 
Cartan subalgebra of a symmetric space is a maximal abelian subalgebra of the subspace 
P (see paragraph 15. 2j) . and the rank of a symmetric space is the number of generators in 
this subalgebra. 

If G is connected and G = K © P where K is a compact symmetric subalgebra, then each 
group element can be decomposed as g = kp (right coset decomposition) or g = pk (left 
coset decomposition), with k E K = e^, p E P = e^. P is not a subgroup, unless it is 
abelian and coincides with its Cartan subalgebra. However, if the involutive automorphism 
that splits the algebra is denoted a, one can show ([HI!, Ch. 6) that gpc^g^^) G P. This 
defines G as a transformation group on P. Since cr{k~^) = k~^ for k E K, this means 

p' = kpk-^ G P (3.4) 

if k E K, p E P. Now suppose there are no other elements in G that satisfy a{g) = g than 
those in K. This will happen if the set of elements satisfying a{g) = g is connected. Then 
P is isomorphic to G/K. Also, G acts transitively on P in the manner defined above (cf. 
subsection 12. 3p . The tangent space of G/K at the origin (identity element) is spanned by 
the subspace P of the algebra. 
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3.3 Radial coordinates 



In this paragraph we define two coordinate systems frequently used on symmetric spaces. 
Let G = K © P be a Cartan decomposition of a semisimple algebra and let Ho C P be a 
maximal abelian subalgebra in the subspace P. Define M to be the subgroup of elements 
in K such that 

M = {keK : kHk-^ = H, H e Hq} (3.5) 

This set is called the centralizer of Hq in K. Under conjugation hj k & K, each element 
H of the Cartan subalgebra is preserved. Further, denote 

M' = {keK : kHk~^ = H\ H, H' e Hq} (3.6) 

This is a larger subgroup than M that preserves the Cartan subalgebra as a whole, but 
not necessarily each element separately, and is called the normalizer of Hq in K. If K is 
a compact symmetric subgroup of G, one can show ([SZI, Ch. 6) that every element p of 
P ~ G/K is conjugated with some element h = for some H G Hq by means of the 
adjoint representation^ of the stationary subgroup K: 

p = khk-^ = kha{k-^) (3.8) 

where k G K/M and H is defined up to the elements in the factor group M' /M. This factor 
group coincides with the Weyl group that was defined in eq. ()2.34|) : since the space Hq can 
be identified with its dual space Hq, we can identify M'/M with the Weyl group of the 
restricted root system (see paragraph 15. 2|1 . The effect of the Weyl group is to transform the 
algebra Hq C P into another Cartan subalgebra Hq C P conjugate with the original one. 
This amounts to a permutation of the roots of the restricted root lattice corresponding to 
a Weyl reflection. Equation ()3.8|) means that every element g & G can be decomposed as 
g = pk = k'hk'~^k = k'hk", and this is very much like the Euler angle decomposition of 
SO{n). 

Thus, if Xq is the flxed point of the subgroup K, an arbitrary point x E P can be written 
^Note that 

e^ilc-^ = c-'^^'H = f2 ^^^^ (3-7) 
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X = khk ^xo = khxo 



(3.9) 



The coordinates {k{x), h{x)) are called spherical coordinates. k{x) is the angular coordinate 
and h{x) is the spherical radial coordinate of the point x. Eq. ()3.8|1 defines the so called 
spherical decomposition of the elements in the coset space. Of course, a similar reasoning 
is true for the space P* ~ G*/K. 

This means every matrix p in the coset space G/K can be diagonalized by a similarity 
transformation by the subgroup K, and the radial coordinates are exactly the set of eigen- 
values of the matrix p. These "eigenvalues" are not necessarily real numbers. This is easily 
seen in the example in eq. ()3.3|1 . It can also be seen in the adjoint representation. Suppose 
the algebra G = K©P is compact. From eq. ()2.13p . in the adjoint representation Hi G Hq 
has the form 



/O ... 



Hi 



Oii 



(3.10) 



-Vi J 



where the matrix is determined by the structure constants {[Hi,Hj] = 0, [Hi,E±a] = 
±aiE±a ■■■ and ±ai, ...,±r]i are the roots corresponding to H^). Since the Killing form 
must be negative (see subsection 13. 4j) for a compact algebra, the coordinates of the Cartan 
subalgebra must be purely imaginary and the group elements corresponding to Hq must 
have the form 



/ 1 



Jt H 



(3.11) 



with t = {t^jt"^, ...t^) and real parameters. In particular, if the eigenvalues are real for 
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p G P*, they are complex numbers for p E P. 

Example: In the example we gave in the preceding subsection, the coset space G* /K 
= SL{n, R)/ SO{n) P* = e*^ consists of real positive-definite symmetric matrices. Note 
that G = K © P implies that G can be decomposed as G = PK and G* as G* = 
P*K. The decomposition G* = P*K in this case is the decomposition of a SL{n,R) 
matrix in a positive-definite symmetric matrix and an orthogonal one. Each positive- 
definite symmetric matrix can be further decomposed: it can be diagonalized by an SO{n) 
similarity transformation. This is the content of eq. ()3.8j) for this case, and we know it to be 
true from linear algebra. Similarly, according to eq. ()3.8j) the complex symmetric matrices 
in G/K = SU{n, C)/ SO{n) ~ P = e^ can be diagonalized by the group K = SO{n) to a 
form where the eigenvalues are similar to those in eq. ()3.1H) . 

In terms of the subspace P of the algebra, eq. (|3.8|) amounts to saying that any two Cartan 
subalgebras Hq, Hq of the symmetric space are conjugate under a similarity transformation 
by and we can choose the Cartan subalgebra in any way we please. However, the number 
of elements that we can diagonalize simultaneously will always be equal to the rank of the 
symmetric space. 

There is also another coordinate system valid only for spaces of the type P* ~ G* / K. This 
coordinate system is called horospheric and is based on the so called Iwasawa decomposition 
P7j of the algebra: 

G = N+©Ho©K (3.12) 

Here K, Hq, N+ are three subalgebras of G. K is a maximal compact subalgebra, Hq is 
a Cartan subalgebra, and 

N+ = E (3.13) 

is an algebra of raising operators corresponding to the positive roots a{H) > with respect 
to Ho (G^ is the space generated by Ea). As a consequence, the group elements can be 
decomposed g = nhk, in an obvious notation. This means that if Xq is the fixed point of 
K, any point x G G*/K can be written 

X = nhkxQ = nhxo (3-14) 
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The coordinates h{x)) are called horospheric coordinates and the element h = h{x) 

is called the horospheric projection of the point x or the horospheric radial coordinate. 



3.4 The metric on a Lie algebra 

A metric tensor can be defined on a Lie algebra jHOl ISH IHHl IHZj • For our purposes, it will 
eventually serve to define the curvature of a symmetric space and be useful in computing 
the Jacobian of the transformation to radial coordinates. In sections El and |H1 we will see 
the importance of this Jacobian in physical applications in connection with the radial part 
of the Laplace-Beltrami operator. 

If {Xi} form a basis for the Lie algebra G, the metric tensor is defined by 

g,, = K{X,, X,) = tr(adX,adX,) = C[,q, (3.15) 

The symmetric bilinear form K{Xi,Xj) is called the Killing form. It is intrinsically asso- 
ciated with the Lie algebra, and since the Lie bracket is invariant under automorphisms of 
the algebra, so is the Killing form. 

Example: The generators Xj = Li, X^ = L2, X2 = L3 of 50(3) given in eq. ()3.3p obey 
the commutation relations [Li,Lj] = CfjL^ = ^e^kL^. From eq. ()3.15|) . the metric for this 
algebra is gij = —^Stj. The generators and the structure constants can be normalized so 
that the metric takes the canonical form g^j = —Sij. 

Just like we defined the Killing form K{Xi,Xj) for the algebra G in eq. ()3.15p using the 
adjoint representation, we can define a similar trace form Kp and a metric tensor gp for 
any representation p by 

gp^.j = Kp{X,,X,) = tT{p{X,)p{X,)) (3.16) 

where p{X) is the matrix representative of the Lie algebra element X. If p is an automor- 
phism of G, Kp{Xi,Xj) = K{X^,Xj). 

Suppose the Lie algebra is semisimple (this is true for all the classical Lie algebras except 
the Lie algebras GL{n,C), U{n,C)). According to Cartan's criterion, the Killing form is 
non-degenerate for a semisimple algebra. This means that det^fj^ 7^ 0, so that the inverse 
of gij, denoted by g^^ , exists. Since it is also real and symmetric, it can be reduced to 
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canonical form Qij = diag(— 1, —1, 1, 1) with p — I's and {n — p) +l's, where n is the 
dimension of the algebra. 

p is an invariant of the quadratic form. In fact, for any real form of a complex algebra, the 
trace of the metric, called the character of the particular real form (see below and in |31j ) 
distinguishes the real forms from each other (though it can be degenerate for the classical 
Lie algebras |3T). The character ranges from — n, where n is the dimension of the algebra, 
to +r, where r is its rank. All the real forms of the algebra have a character that lies 
in between these values. In subsection 14.11 we will see several explicit examples of Killing 
forms. 

A famous theorem by Weyl states that a simple Lie group G is compact, if and only if the 
Killing form on G is negative definite. Otherwise it is non-compact. This is actually quite 
intuitive and natural (see [211 , Ch. 9, paragraph 1.2). On a compact algebra, the metric 
can be chosen to be minus the Killing form, if it is required to be positive-definite. 

The metric on the Lie algebra can be extended to the whole coset space P ~ G/K, 
P* ~ G*/K as follows. At the origin of G/K and G*/K, the identity element /, the 
metric is identified with the metric in the algebra, restricted to the respective tangent 
spaces P, iP. Since the group acts transitively on the coset space (cf. paragraph I2.3|l . and 
the orbit of the origin is the entire space, we can use a group transformation to map the 
metric at the origin to any point M in the space. The metric tensor at M will depend on 
the coset representative M. It is given by 



where gij{I) is the metric at the origin (identity element) of the coset space. (j3.17j) follows 
from the invariance of the line element ds"^ = gijdx^dx^ under translations. If {Xi} is a 
basis in the tangent space, and dM = exp{dx'^Xi) is a coset representative infinitesimally 
close to the identity, we need to know how dx^ transforms under translations by the coset 
representative M. We will not discuss that here, but some generalities can be found for 
example in Ch. 9, paragraph V.4. of ref. f^. In general, it is not an easy problem unless 
the coset has rank 1. 

Example: The line element ds"^ on the radius-1 2-sphere S0{3)/ S0{2) in polar coordi- 
nates is ds"^ = dO"^ + s\.v?0 dcj)^. The metric at the point (^, (p) is 




9« = I „ ..... 1 . 9'' = (318) 
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where the rows and columns are labelled in the order 6, 0. 

The distance between points on the symmetric space is defined as follows. The length of 
a vector X = J^it^^i in the tangent space P (this object is well-defined because P is 
endowed with a definite metric) is identified with the length of the geodesic connecting the 
identity element in the coset space with the element M = exp(X) |31j . 



3.5 The algebraic structure of symmetric spaces 

Except for the two algebras SL(n, R) and SU*(2n) (and their dual spaces related by the 
Weyl trick), for which the subspace representatives of K, P and iP consist of square, 
irreducible matrices (for SL(n, R), we saw this in the example in subsection 13.11 and for 
SU(n, C) explicitly in eq. ()3.3j] ). the matrix representatives of the subalgebra K and of the 
subspaces P and iP in the fundamental representation consist of block-diagonal matrices 
X e K, F e P, F' G zP of the form ^ 



in the Cartan decomposition. Here = —A, B'^ = —B and C = iC. In fact, for any 
finite-dimensional representation, the matrix representatives of K and P are antihermitean 
(thus they become antisymmetric if the representation of P is real) and as a consequence, 
those of iP are hermitean (symmetric in case the representation of iP is real) [31]. This is 
true irrespective of whether the matrix representatives are block-diagonal or square. 

The exponential maps of the subspaces P and iP are isomorphic to coset spaces G/K 
and G*/K, respectively (see for example [301 EZ|)- The exponential map of the algebra 
maps the subspaces P and iP into unitary and hermitean matrices, respectively. In the 
fundamental representation, these spaces are mapped onto [HI] 




(3.19) 




(3.20) 
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where X is a spherical and X a hyperbohc function of the submatrix C (C): 



This shows exphcitly that the range of parameters parametrizing the two cosets is bounded 
for the compact coset and unbounded for the non-compact coset, respectively. We already 
saw an explicit example of these formulas in subsection 12.31 



4 Real forms of semisimple algebras 



In this section we will introduce the tools needed to find all the real forms of any (semi) simple 
algebra. The same tools will then be used in the next section to find the real forms of a 
symmetric space. When thinking of a real form, it is convenient to visualize it in terms of 
its metric. As we saw in paragraph 13 . 41 the trace of the metric is called the character of the 
real form and it distinguishes the real forms from each other. In the following subsection 
we discuss various real forms of an algebra and we see how to go from one form to another. 
In each case, we compute the metric and the character explicitly. We also give the simplest 
possible example of this procedure, the rank-1 algebra. In subsection 14.21 we enumerate 
the involutive automorphisms needed to classify all real forms of semisimple algebras and 
again, we illustrate it with two examples. 



4.1 The real forms of a complex algebra 

In general a semisimple complex algebra has several distinct real forms. Recall from sub- 
section 12.51 that a real form of an algebra is obtained by taking linear combinations of its 
elements with real coefficients. The real forms of the complex Lie algebra G 



Y.c'Hi + Y.C'Ea, (c\ complex), (4.1) 

i a 

where Hq = {Hi} is the Cartan subalgebra and {E±ct} are the sets of raising and lowering 
operators, can be classified according to all the involutive automorphisms of G obeying 
(T^ = 1. Two distinctive real forms are the normal real form and the compact real form. 
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The normal real form of the algebra ()4.1|) . which is also the least compact real form, 
consists of the subspace in which the coefficients d, are real. The metric in this case 
with respect to the bases {Hi,E±a} is (with appropriate normalization of the elements of 
the Lie algebra to make the entries of the metric equal to ±1) 



/ 1 



9ij 



1 

1 



1 

1 OJ 



(4.2) 



where the r I's on the diagonal correspond to the elements of the Cartan subalgebra (r is 
obviously the rank of the algebra), and the 2x2 matrices on the diagonal correspond to 
the pairs £"-1-0, of raising and lowering operators. This structure reflects the decomposition 
of the algebra G into a direct sum of the root spaces: G = Hq © J^a^a- This metric 
tensor can be transformed to diagonal form, if we choose the generators to be 



Example: In our example with SU(3, C), K and iP are exactly the subspaces spanned by 
{X2, X5, Xj} and {iXi, iX^, 1X4^, iXq, iXg} (cf. eq. (j3.3|) ). and {Ea—E_a) and —i{Ea+E^a) 
are exactly the Gell-Mann matrices (cf. eq. ()2.30|) ). 

Then gij takes the form 



/ 1 



9ij 



1 
-1 



(4.4) 



1 




-1/ 
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where the entries with a minus sign correspond to the generators of the compact subalgebra 
K, the first r entries equal to +1 correspond to the Cartan subalgebra, and the remaining 
ones to the operators in iP not in the Cartan subalgebra. This is the diagonal metric 
tensor corresponding to the normal real form. The character of the normal real form is 
plus the rank of the algebra. 

The compact real form of G is obtained from the normal real form by the Weyl unitary 
trick: 



The character of the compact real form is minus the dimension of the algebra, and the 
metric tensor is gij = diag(— 1, —1). 

Example: We will use as an example the well-known SU(2, C) algebra with Cartan 
subalgebra Hq = {^3} and raising and lowering operators {J±}- We have chosen the 
normalization such that the non-zero entries of Qij are all equal to 1: 



■^3 = 573^3, J± = l{n±iT2) (4.6) 
where in the defining representation of SU(2, C) 



-3= , . -1= . „ (4.7) 






The normalization is such that 



[Js, J±] = ±j^^J±, [J+,J_] = ^J3 (4.8) 

In equation ()2.16p we constructed the adjoint representation of this algebra, albeit with a 
different normalization. Using the present normalization to set the entries of the metric 
equal to 1, we see that the non-zero structure constants are C^_^ = — C+3 = —C^_ = Cl^ = 
C'l_ = — C^^_|_ = The entries of the metric are given by eq. g^j = K{Ji, Jj) = 

CigCjr with summation over repeated indices, so we see that the metric of the normal real 
form SU(2, R) in this basis is 
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9ij 



/ 1 \ 

1 
Vo 1 0^ 



(4.9) 



where the rows and columns are labelled by 3, +, — respectively. This corresponds to 
eq. (021). 

To pass now to a diagonal metric, we just have to set 



Si 

S2 



S3 — J3 

J++J- _ 1 



V2 



J+-J- 



2y/2 



V2 



2V2 



T2 



(4.10) 



like in eq. ()4.3j) . The commutation relations then become 



[Si,S2] = -^S 



[S2, S3] - -:^Si, [S3, Si] - ^^^S 



V2' 



(4.11) 



These commutation relations characterize the algebra SO(2, 1; R). From here we find the 
structure constants = — = C23 = — C32 = — = Cf^ ~ diagonal 
metric of the normal real form with rows and columns labelled 3,1,2 (in order to comply 
with the notation in eq. ()4.4|) ) is 



9ij 




(4.12) 



which is to be compared with eq. ()4.4p . According to eq. ()4.3|) . the Cartan decomposition 
of G* is G* = K © iP where K = {S2} and iP = {S3,Si}. The Cartan subalgebra 
consists of S3. 



Finally, we arrive at the compact real form by multiplying S3 and Si with i. Setting 
zSi = Si, S2 = S2, iS3 = S3 the commutation relations become those of the special 
orthogonal group: 
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[Si, S2] = -^S3, S3] = -;^Si, [S3, Si] = -^S2 (4.13) 

The last commutation relation in eq. ()4.11|) has changed sign whereas the others are un- 
changed. C|]^, C^3, and consequently (733 and gn change sign and we get the metric for 
SO(3,R): 



9ij 





V 








(4.14) 



This is the compact real form. The subspaces of the compact algebra G = K © P are 
K = {S2} and P = {S3, Si}. Weyl's theorem states that a simple Lie group G is compact, 
if and only if the Killing form on G is negative definite; otherwise it is non-compact. In 
the present example, we see this explicitly. 



4.2 The classification machinery 



To classify all the real forms of any complex Lie algebra, with characters lying between 
the character of the normal real form and the compact real form (the intermediate real 
forms obviously have an indefinite metric), it suffices to enumerate all the involutive au- 
tomorphisms of its compact real form. A detailed and almost complete account of these 
procedures for the non-exceptional groups can be found in [21], Chapter 9, paragraph 3. 
To summarize, if G is the compact real form of a complex semisimple Lie algebra G^, G* 
runs through all its associated non-compact real forms G*, G'*, ... with corresponding 
maximal compact subgroups K, K', ... and complementary subspaces iP, iP\ ... as a 
runs through all the involutive automorphisms of G. 

One such automorphism is complex conjugation ai = K, which is used to split the compact 
real algebra into subspaces K and P in eq. (j4.5p . (To avoid confusion: the generators can 
be complex even though the field of real numbers is used to multiply the generators in a 
real form of an algebra. If the generators are also real, we speak of a real representation. 
However, whether we consider the field to be R and the generators to be complex, or the 
opposite, also depends on our definition of basis; cf. one of the footnotes in subsection [231). 
The involutive automorphisms a satisfy aGa~^ = G, cr^ = 1, which implies that a either 
commutes or anticommutes with the elements of the compact algebra G: if aXa~^ = X', 
then aX'a~^ = X, and we get X' = ±X for X, X' G G (see the example below). One 
can show [2H] (Ch. VII), that it suffices to consider the following three possibilities for cr: 
(Ti = K, a2 = Ip,q and = Jp^p where 
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(4.15) 



and Ip denotes the p x p unit matrix. By operating with one (or two successive ones) of 
these automorphisms on the elements of G, we can construct the subspaces K and P, and 
K and iP of the corresponding non-compact real form G*. A complex algebra and all its 
real forms (the compact and the various non-compact ones) correspond to the same root 
lattice and Dynkin diagram. 

Example: The normal real form of the complex algebra G^ = SL(n, C) is the non- 
compact algebra G* = SL(n, R). As we saw in subsection 13.11 this algebra can be decom- 
posed as K © iP where K is the algebra consisting of real, skew-symmetric and traceless 
n X n matrices and iP is the algebra consisting of real, symmetric and traceless n x n ma- 
trices. Under the Weyl unitary trick we constructed, in a previous example, this algebra 
from the compact real form of G*^, SU(n, C) = G = K © P. 

Starting with the compact real form G, we can construct all the various non-compact 
real forms G*, G'*,... from it, by applying the involutive automorphisms ui, (J2, o"3 to the 
elements of G. All the real forms related to the root system An-i are obtained by applying 
the three involutions to G = SU(n, C): 

(Ji) The involutive automorphism ai = K (complex conjugation) splits G = SU(n, C) 
into K © P (we recall this from the example in paragraph 13. 1|) . The non-compact real 
form obtained this way, by the Weyl unitary trick, is exactly the normal real form G* = 
K©iP = SL(n,R). 

(J2) A general matrix in the Lie algebra SU(n, C) can be written in the form 



where A, C are complex pxp and qxq matrices satisfying = —A, = —C, trA+trC = 
(since the determinant of the group elements must be +1), and B is an arbitrary complex 
px q matrix [p + q = n). In eq. (j4.16p . the matrices A, B and C are all linear combinations 
of submatrices in both subspaces K = {^{E^^ — E_a^)} and P = {iHj^ ^(£'„. -|- E_a^)}. 
The action of the involution (T2 = Ip,q on X is 




(4.16) 




(4.17) 
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Therefore, we see that the subspaces K' and P' are given by the matrices 



Indeed, we see that Jp ^ transforms the Lie algebra elements in K' into themselves, and 
those in P' into minus themselves. The transformation by Ip q mixes the subspaces K and 
P, and splits the algebra in a different way into K' © P'. The matrices 



A iB 
-iB^ C 



GK'ezP' (4.19) 



define the non-compact real form G'*. This algebra is called SU(p, q; C) and its maximal 
compact subalgebra K' is SU(p) ® SU(q) ® U(l). 

(T3) By the involutive automorphism a^ai = Jp^pK one constructs in a similar way (for 
details see [HI]) a third non-compact real form (for even n = 2p) G"* = K"©iP" associated 
to the algebra G = SU(2p, C). G"* is the algebra SU*(2p) and its maximal compact 
subalgebra is USp(2p). 

This procedure, summarized in the formula below, exhausts all the real forms of the simple 
algebras. 



G* = K © iP 



G*^ ^ G = K © P ^ G'* = K' © ?P' (4.22) 



0"2 

\ G"* = K " © iP" 



^The algebra SU*(2p) is represented by complex 2p x 2p matrices of the form 

(4.20) 

where tiA+trA* — 0. USp(2p) denotes the complex 2px 2p matrix algebra of the group with both unitary 
and symplectic symmetry (USp(2p, C) can also be denoted U(p, Q) where Q is the field of quaternions). 
A matrix in the algebra USp(2p, C) can be written as 

where = —A, = B, and the superscript ^ denotes reflection in the minor diagonal. 
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Example: Note that it may not always be possible to apply all the above involutions 
cr2, CTs to the algebra. For example, complex conjugation a\ does not do anything to 
SO(2n + 1, R), because it is represented by real matrices, neither is cxs a symmetry of this 
algebra, since the adjoint representation is odd-dimensional and (T3 has to act on a 2p x 2p 
matrix. The only possibility that remains is 02 = Ip,q- For a second, even more concrete 
example, let's look at the algebra SO(3, R), belonging to the root lattice Bi. This algebra 
is spanned by the generators Li, L2, given in subsection 12.31 A general element of the 
algebra is 



X = t ■ 



L = i 

2 



-t3 



(4.23) 



This splitting of the algebra is caused by the involution l2^i acting on the representation: 




and it splits it into SO(3) = K©P = SO(2) © SO(3)/SO(2). Exponentiating, as we saw 
in subsection 12.31 the coset representative is a point on the 2-sphere 



M 



sm^(ti)2 


+ (t2)2 






sm^(ti)2 


+ (t2)2 




(t2)2 



v 



COS^(tl)2 + (t2)2 ) 



/ 

V 



2 I 2 I 2 

X + y + z 



(4.25) 



By the Weyl unitary trick we now get the non-compact real form G* = K©iP: SO(2, 1) 
SO(2) © SO(2, l)/SO(2). This algebra is represented by 



It' 



it' 
It" 



-it^ -it^ 



-it^ -it^ 



(4.26) 



and after exponentiation of the coset generators 
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M 



,^2 smhV(fl)2 + (f2)2 \ 
,^i sinhv/(fi)2 + (f2)2 

coshy^(ti)2 + (t2)2 y 



(4.27) 



The surface in consisting of points (x, y, z) satisfying this equation is the hyperboloid 
H'^. Similarly, we get the isomorphic space S0{1,2) / S0{2) by applying J12: SO(l,2) = 
K © zP = SO(2) © SO(l, 2)/SO(2) and in terms of the algebra 



X = - 
2 





1 

®2 


/ 


f 











-if 



(4.28) 



5 The classification of symmetric spaces 



In this section we introduce the curvature tensor and the sectional curvature of a symmetric 
space, and we extend the family of symmetric spaces to include also flat or Euclidean-type 
spaces. These are identified with the subspace P of the Lie algebra itself, and the group 
that acts on it is a semidirect product of the subgroup K and the subspace P. As we 
will learn, to each compact subgroup K corresponds a triplet of symmetric spaces with 
positive, zero and negative curvature. The classification of these symmetric spaces is in 
exact correspondence with the new classification of random matrix models to be discussed 
in Part II. These spaces exhaust the Cartan classification and have a definite metric. They 
are listed in Table ^ together with some of their properties. 

In paragraph 15.21 we introduce restricted root systems. In the same way as a Lie algebra 
corresponds to a given root system, the "algebra" (subspace P or iP) of each symmetric 
space corresponds to a restricted root system. These root systems are of primary impor- 
tance in the physical applications to be discussed in Part II. The restricted root system can 
be of an entirely different type from the root system inherited from the complex extension 
algebra, and its rank may be different. We work out a specific example of a restricted 
root system as an illustration. In spite of their importance, we have not been able to find 
any explicit reference in the literature that explains how to obtain the restricted root sys- 
tems. Instead, we found that they are often referred to in tables and in mathematical texts 
without explicitly mentioning that they are restricted, which could easily lead to confusion 
with the inherited root systems. In reference [HI] the root system that is associated to 
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each symmetric space is the one inherited from the complex extension algebra, whereas for 
example in Table Bl of reference jH] and in |3H| the restricted root systems are listed. 

There are also symmetric spaces with an indefinite metric, so called pseudo-Riemannian 
spaces, corresponding to a maximal non-compact subgroup H. For completeness, we will 
briefly discuss how these are obtained as real forms of symmetric spaces corresponding to 
compact symmetric subgroups. This does not require any new tools than the ones we have 
already introduced, namely the involutive automorphisms. 



5.1 The curvature tensor and triplicity 

Suppose that K is a maximal compact subalgebra of the non-compact algebra G* in the 
Cartan decomposition G* = K © iP, where iP is a complementary subspace. K and P 
(alternatively K and iP) satisfy eq. ()3.2p : 



K is called a symmetric subalgebra and the coset spaces exp(P) ^ G/K and exp(zP) ~ 
G* /K are globally symmetric Riemannian spaces. Globally symmetric means that every 
point on the manifold can be moved to any other point by a particular group operation (we 
discussed this in paragraph 12.31 for a rigorous definition of globally symmetric spaces see 
Helgason [^0], paragraph IV. 3). In the same way, the metric can be defined in any point 
of the manifold by moving the metric at the origin to this point, using a group operation 
(cf. eq. ()3.17p in paragraph 13. 4|) . The Killing form restricted to the tangent spaces P and 
iP at any point in the coset manifold has a definite sign. The manifold is then called 
"Riemannian". The metric can be taken to be either plus or minus the Killing form so 
that it is always positive definite (cf. paragraph I3.4|l . 

A curvature tensor with components Rjp.i can be defined on the manifold G/K or G*/K 
in the usual way pm I32j. It is a function of the metric tensor and its derivatives. It was 
proved for instance in jSU], Ch. IV, that the components of the curvature tensor at the 
origin of a globally symmetric coset manifold is given by the expression 



where {Xi} is a basis for the Lie algebra. The sectional curvature at a point p is equal to 



[K,K]CK, [K,P]CP, [P,P]CK 



(5.1) 




(5.2) 
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lC = g{[[X,YlXlY) 



(5.3) 



where g is an arbitrary symmetric and nondegenerate metric (such a metric is also called 
a pseudo-Riemannian structure, or simply a Riemannian structure if it has a definite 
sign) on the tangent space at p, invariant under the action of the group elements. In 
()5.3|) . g{Xi,Xj) = gij and {X, y } is an orthonormal basis for a two-dimensional subspace 
S of the tangent space at the point p (assuming it has dimension > 2). The sectional 
curvature is equal to the gaussian curvature on a 2-dimensional manifold. If the manifold 
has dimension > 2, ()5.3|1 gives the sectional curvature along the section 5*. 

Eqs. ()5.2p and ()5.3|) . together with eq. ()5.1|) show that the curvature of the spaces G/K 
and G* /K has a definite and opposite sign ( 30j, par. V.3). Thus, we see that if G is a 
compact semisimple group, to the same subgroup K there corresponds a positive curvature 
space P ~ G/K and a dual negative curvature space P* ~ G* /K. The reason for this 
is exactly the same as the reason why the sign changes for the components of the metric 
corresponding to the generators in iP as we go to the dual space P. We remind the reader 
that the sign of the metric can be chosen positive or negative for a compact space. The 
issue here is that the sign changes in going from G* /K to G/K. 

Example: We can use the example of SU (2) in paragraph 14.11 to see that the sectional 
curvature is the opposite for the two spaces G/K and G* /K. If we take {X, Y} = {S3, Si} 
as the basis in the space iP and {S3, Si} (Sj = zS^) as the basis in the space P, we see 
by comparing the signs of the entries of the metrics we computed in eqs. ()4.12|1 and ()4.14|1 
that the sectional curvature /C at the origin has the opposite sign for the two spaces 
SO{2,l)/SO{2) and ^0(3)/^0(2). 

Actually, there is also a zero-curvature symmetric space X° = G^/K related to X~^ = G/K 
and X" = G* /K, so that we can speak of a triplet of symmetric spaces related to the same 
symmetric subgroup K. The zero-curvature spaces were discussed in [H] and in Ch. V of 
Helgason's book [301, where they are referred to as "symmetric spaces of the euclidean 
type". That their curvature is zero was proved in Theorem 3.1 of jSUj, Ch. V. 

The fiat symmetric space X° can be identified with the subspace P of the algebra. The 
group G^ is a semidirect product of the subgroup K and the invariant subspace P of the 
algebra, and its elements g = {k, a) act on the elements of in the following way: 

g{x) = kx + a, k e K, x, a G X° (5.4) 

if the x's are vectors, and 



43 



g{x) = kxk + a, k & K, x, a G X' 



(5.5) 



if the x's are matrices. We will see one example of each below. 

The elements of the algebra P now define an abelian additive group, and X° is a vector space 
with euclidean geometry. In the above scenario, the subspace P contains only the operators 
of the Cartan subalgebra and no others: P = Hq, so that P is a subalgebra of G°. The 
algebra G° = K © P belongs to a non-semisimple group G^, since it has an abelian ideal 
P: [K, K] C K, [K, P] C P, [P, P] = 0. Note that K and P still satisfy the commutation 
relations 1)5.11) . In this case the coset space is fiat, since by ()5.1|) . i?^;, = for all the 
elements X e P. Eq. ()5.2|) is valid for any space with a Riemannian structure. Indeed, it 
is easy to see from eqs. ()5.2|) . ()5.3|) that i?^-^ = /C = if the generators are abelian. Even 
though the Killing form on non-semisimple algebras is degenerate, it is trivial to find a 
non-degenerate metric on the symmetric space X" that can be used in ()5.Hj) to find that 
the sectional curvature at any point is zero. For example, as we pass from the sphere to 
the plane, the metric becomes degenerate in the limit as [Li, L2] ~ L3 — [Pi, P2] = (see 
the example below). Obviously, we do not inherit this degenerate metric from the tangent 
space on like in the case of the sphere, but the usual metric for R^, gij = 6ij provides 
the Riemannian structure on the plane. 

Examples: An example of a flat symmetric space is E2/K, where = E2 is the euclidean 
group of motions of the plane R^: g{x) = kx + a, g = {k, a) G G^ where k E K = 5*0(2) 
and a G R^. The generators of this group are translations Pi, P2 E Hq = P and a rotation 
J G K satisfying [Pi, P2] = 0, [J, Pi] = -e'^Pj, [J, J] = 0, in agreement with eq. (fOll 
defining a symmetric subgroup. The abelian algebra of translations J2f=ifPi, G R, is 
isomorphic to the plane R^, and can be identified with it. 

The commutation relations for E2 are a kind of limiting case of the commutation relations 
for SO(3) ~ SU(2) and SO(2, 1). If in the limit of infinite radius of the sphere S"^ we 
identify Si with Pi, E2 with P2, and S3 with J, we see that the commutation relations 
resemble the ones described in eq. ()4.11|) and ()4.13p - we only have to set [Ei,E2] = 0, 
which amounts to setting = —C21 —>■ 0. From here we get the degenerate metric of the 
non-semisimple algebra E2: 



/ -1 



9ij 



(5.6) 



where the only nonzero element is (733. This is to be confronted with eqs. ()4.12|) and ()4.14|) 
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which are the metrics for SO(2, 1) and SO(3). This is an example of contraction of an 
algebra. 

An example of a triplet {X~^ , , X~} corresponding to the same subgroup K = SO{n) 
is: 

1) X^ = SU{n,C)/ SO{n), the set of symmetric unitary matrices with unit determinant; 
it is the space exp(P) where P are real, symmetric and traceless n x n matrices. (Cf. the 
example in subsection 13.11 ) 

2) X^ is the set P of real, symmetric and traceless n xn matrices and the non-semisimple 
group (7° is the group whose action is defined by g{x) = kxk~^ + a, g = {k, a) G where 
k & K = SO{n) and x, a G X". The involutive automorphism maps g = {k,a) G G° into 
g' = (k, -a). 

3) X^ = SL{n, R)/ SO{n) is the set of real, positive, symmetric matrices with unit deter- 
minant; it is the space exp(zP) where P are real, symmetric and traceless n x n matrices. 

We remark that the zero-curvature symmetric spaces correspond to the integration mani- 
folds of many known matrix models with physical applications. 

The pairs of dual symmetric spaces of positive and negative curvature listed in each row 
of Table ^ originate in the same complex extension algebra |S5 with a given root lattice. 
This "inherited" root lattice is listed in the first column of the table. In our example in 
paragraph 14.21 this was the root lattice of the complex algebra G'" = SL(n, C). The same 
root lattice An-i characterizes the real forms of SL(n, C): as we saw in the example these 
are the algebras SU(n, C), SL(n, R), SU(p,q;C) and SU*(2n), and we have seen how 
to construct them using involutive automorphisms. 

However, also listed in Table^is the restricted root system corresponding to each symmetric 
space. This root system may be different from the one inherited from the complex extension 
algebra. Below, we will define the restricted root system and see an explicit example of one 
such system. While the original root lattice characterizes the complex extension algebra 
and its real forms, the restricted root lattice characterizes a particular symmetric space 
originating from one of its real forms. The root lattices of the classical simple algebras 
are the infinite sequences y4„, C^, -Dn, where the index n denotes the rank of the 
corresponding group. The root multiplicities rrio, mi, nig listed in Table ^ (where the 
subscripts refer to ordinary, long and short roots, respectively) are characteristic of the 
restricted root lattices. In general, in the root lattice of a simple algebra (or in the graphical 
representation of any irreducible representation), the roots (weights) may be degenerate 
and thus have a multiplicity greater than 1. This happens if the same weight /i = (/zi, Hr) 
corresponds to different states in the representation. In that case one can arrive at that 
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particular weight using different sets of lowering operators E^a on the highest weight of 
the representation. Indeed, we saw in the example of SU{3,C) in subsection 12.51 that the 
roots can have a multiplicity different from 1. The same is true for the restricted roots. 

The sets of simple roots of the classical root systems (briefly listed in subsection I2.7|l 
have been obtained for example in j^H EE]- In the canonical basis in R", the roots of 
type {icj ± ej,i ^ j} are ordinary while the roots {±2ej} are long and the roots {icj} 
are short. Only a few sets of root multiplicities are compatible with the strict properties 
characterizing root lattices in general. 



5.2 Restricted root systems 

The restricted root systems play an important role in connection with matrix models and 
integrable Calogero-Sutherland models (these models will be introduced in section [7j). We 
will discuss this in detail in Part II. In this subsection we will explain how restricted root 
systems are obtained and how they are related to a given symmetric space. ^ 

As we have repeatedly seen in the examples using the compact algebra SU(n, C) (in 
particular in subsection 14. 2p . the algebra SU(p, q; C) (p + q = n) is a non-compact real 
form of the former. This means they share the same rank-(n — 1) root system An-i- 
However, to the symmetric space SU{p,q;C)/{SU{p) ® SU{q) ® U{1)) one can associate 
another rank-r' root system, where r' = mm{p, q) is the rank of the symmetric space. 
For some symmetric spaces, it is the same as the root system inherited from the complex 
extension algebra (see Tabled for a list of the restricted root systems), but this need not be 
the case. For example, the restricted root system is, in the case of SU{p,q;C)/{SU{p) 
SU{q) ® ?7(1)), BCr'- When it is the same and when it is different, as well as why the rank 
can change, will be obvious from the example we will give below. 

In general the restricted root system will be different from the original, inherited root 
system if the Cartan subalgebra is a subset of K. The procedure to find the restricted root 
system is then to define an alternative Cartan subalgebra that lies partly (or entirely) in P 
(or iP). 

To achieve this, we first look for a different representation of the original Cartan subalgebra, 
that gives the same root lattice as the original one (i.e., for the SU(p, q; C) algebra). 
In general, this root lattice is an automorphism of the original root lattice of the same kind, 
obtained by a permutation of the roots. Unless we find this new representation, we will 

^The authors are indebted to Prof. Simon Salamon for explaining how the restricted root systems are 
obtained. 
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not be able to find a new, alternative Cartan subalgebra that lies partly in the subspace 
P. 



Once this has been done, we take a maximal abelian subalgebra of P (the number of 
generators in it will be equal to the rank r' of the symmetric space G/K or G*/K) and 
find the generators in K that commute with it. These generators will be among the ones 
that are in the new representation of the original Cartan subalgebra. These commuting 
generators now form our new, alternative Cartan subalgebra that lies partly in P, partly 
in K. Let's call it Aq. 

The new root system is defined with respect to the part of the maximal abelian subalgebra 
that lies in P. Therefore its rank is normally smaller than the rank of the root system 
inherited from the complex extension. We can define raising and lowering operators E'^ in 
the whole algebra G that satisfy 



[X:,E',] = a[E', (X^G AoHP) 



(5.7) 



The roots a'^ define the restricted root system. 

Example: Let's now look at a specific example. We will start with the by now familiar 
algebra SU(3, C). As before, we use the convention of regarding the Tj's as the generators, 
without the factor of i (recall that the algebra consists of elements of the form J2a f^^a = 
iJ2a't"'Ta; cf. the footnote in conjuction with eq. ()2.25|) ). In subsection 12.51 we explicitly 
constructed its root lattice A2. Let's write down the generators again: 



/ 1 
1 

V 




/ 
-i 
\0 i 



To 




I 1 


VO 







/ 
1 
1 



(5. 



The splitting of the SU(3, C) algebra in terms of the subspaces K and P was given in 
eq. (Q: 
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K = {iT2, iT^, zTy}, P = {iTi, zTs, iT^, iT^, iTg} 



(5.9) 



The Cartan subalgebra is {iT^jiT^}. The raising and lowering operators were given in 
in terms oi Tf. 



i?±(i,o) = 75(Ti±zT2) 

^±(i,^) = 75m±^T5) (5.10) 

Now let us construct the Cartan decomposition of G'* = K' © iP' = SU(2,1;C). We 
know from paragraph 14.21 that K' and P' are given by matrices of the form 





e K', m n e P' (5.11) 



where A and C are antihermitean and tiA + trC = 0. Combining the generators to form 
this kind of block-structures (or alternatively, using the involution o"2 = /2,i) we need 
to take linear combinations of the Xj's, with real coefficients, and we then see that the 
subspaces K' and iP' are spanned by 



K' 



iP' 



/ 1 \ 




1 




V o; 




\,iT3,iTs} 












\ 1 , 





1 
-1 



' 2 



/ 1 
-1 
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2V3 
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V 



-2 



(5.12) 
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lo 1 



{^4, ^5, Tq^ T^] 



where the block-structure is evidenced by leaving blank the remaining zero entries. K' 
spans the algebra of the symmetric subgroup SU{2) ® U{1) and iP' spans the comple- 



48 



mentary subspace corresponding to the symmetric space SU{2, 1)/{SU{2) ® U{1)). iP' is 
spanned by matrices of the form 



B 
fit 



(5.13) 



We see that the Cartan subalgebra zHq = {iT^,iT^} hes entirely in K'. It is easy to see 
that by using the alternative representation 



T' 



of the Cartan subalgebra (note that this is a valid representation of SU(3, C) generators) 
while the other Tj's are unchanged, we still get the same root lattice A2. The eigenvec- 
tors under the adjoint representation, the -Eq's, are still given by eq. (j5.1(J|) . However, their 
eigenvalues (roots) are permuted under the new adjoint representation of the Cartan subal- 
gebra, so that they no longer correspond to the root subscripts in ()5.10p . This permutation 
is a Weyl reflection; more specifically, it is the reflection in the hyperplane orthogonal to 
the root (-i, ^). 

Now we choose the alternative Cartan subalgebra to consist of the generators T4, Tg: 



Ao = {T4,T^}, [T^X] = 0, iT^ e P', iT'^ E K' (5.15) 

(Note that unless we first take a new representation of the original Cartan subalgebra, we 
are not able to find the alternative Cartan subalgebra that lies partly in P'.) The restricted 
root system is now about to be revealed. We define raising and lowering operators E'^ in 
the whole algebra according to 



E'^, ~ (T5 ± zTs) ~ (Te ± zTa) E'^^ ~ (T7 ± tT,) (5.16) 

2 2 

The ±a subscripts are the eigenvalues of T4 G iP' in the adjoint representation: 



[T,, E'^,] = ±E'^„ [T,, E'^.] = ±^E'^., [T,, E'^r] = ±^E^ (5.17) 
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These roots form a one-dimensional root system of type BCi. We see that the multiphcity 
of the long roots is 1 and the multiplicity of the short roots is 2 = 2{p — q). This result is 
general (cf. Table If we had ordinary roots, their multiplicity would be 2, but for this 
low-dimensional group we can have only 3 pairs of roots. Note that we can rescale the 
lengths of all the roots together by rescaling the operator T4 in fl5.17j) . but their characters 
as long and short roots can not change. The root system BCi is with respect to the part 
of the Cartan subalgebra lying in iP' only, thus it is called restricted. 

According to eq. ()3.8j) . every element p of P c:^ G/K is conjugated with some element 
h = {H G Ho) through p = khk~^, where k G K/M' and H is defined up to the 
elements in the factor group M' /M. Thus, the decomposition p = khk~^ is not unique. 
The factor group M'/M transforms a Cartan subalgebra Hq C P into another Cartan 
subalgebra Hq C P conjugate with the original one. This amounts to a permutation of 
the roots of the restricted root lattice corresponding to Weyl reflections. The factor group 
M'/M then coincides with the Weyl group of the restricted root system. If we fix the Weyl 
chamber of if, H is unique and k is defined up to transformations by the subgroup M. 
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Table 1: Irreducible symmetric spaces of positive and negative curvature originating in 
simple Lie groups. In the fourth and fifth columns the symmetric spaces G/K and G* /K 
are listed for all the entries except the simple Lie groups themselves, for which the sym- 
metric spaces G and G^ /G are listed. Note that there are also zero curvature spaces 
corresponding to non-semisimple groups and isomorphic to the subspace P of the algebra, 
when P is an abelian invariant subalgebra. These are not listed in the table, but can be 
constructed as explained in subsection 15.11 The root multiplicities listed pertain to the 
restricted root systems of the pairs of dual symmetric spaces with positive and negative 
curvature. 



Root 

space 


Restricted 
root space 


Cartan 
class 


G/K (G) 


G*/K {G^/G) 


mo 


nil 




An-1 


An-1 
An-1 
An-1 

BCq {p>q) 
Cq {p=q) 


A 

AI 

All 

AIII 


SU{N) 

SU{N) 
SO{N) 
SU{2N) 
USp{2N) 

SUip+q) 


SL(N,C) 

SU{N) 
SL(N,R) 

SO{N) 
SU*{2N) 
USp{2N) 

SUip,q) 


2 
1 
4 

2 






1 






2(p-g) 


SU{p)xSU{q)xU{l) 


SU{p)xSU{q)xU{l) 


Bn 


Bn 


B 


S0{2N+1) 


S0(2N+1,C) 
SO{2N+l) 


2 





2 


Cn 


Cn 
Cn 

BCq (p>q) 
Cq (p=q) 


C 
CI 

CII 


USp{2N) 

USp{2N) 


Sp(2N,C) 
USp{2N) 
Sp{2N,R) 


2 
1 

4 


2 
1 

3 





4{p - q) 


SU{N)xU(l) 
USp{2p+2q) 


SU{N)xU{l) 
USp{2p,2q) 


USp{2p)xUSp(2q) 


USp{2p)xUSp(2q) 


Dn 


Dn 
Cn 
BCn 


D 

Dlll-even 
Din-odd 


S0{2N) 

so{m) 


SO{2N,C) 
SO(2N) 
SO'{4N) 


2 

4 
4 



1 
1 




4 


SU{2N)xU{l) 
SO{4N+2) 


SU{2N)xU{l) 
SO*(4N+2) 


SU(2N+l)xU(l) 


SU{2N+l)xU{l) 


Bn (p+g=2iv+i) 

Dn (p+q=2N) 


Bq (p>q) 
Dq {p=q) 


BDI 


SO{p+q) 


SO{p,q) 


1 





p-q 


SO(p)xSO(q) 


SO{p)xSO{q) 
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5.3 Real forms of symmetric spaces 



Involutive automorphisms were used to split the algebra G into orthogonal subspaces to 
obtain the real forms G, G*, G'*... of a complex extension algebra G*". By re-applying the 
same involutive automorphisms to the spaces K, P, and iP, these spaces with a definite 
metric tensor can in turn be split into subspaces with eigenvalue +1 and —1 under this 
new involutive automorphism r. Thus, 



a 


: G - 


-^K©P, 


G* 


= K© 


iP 


T ; 


: K - 


->Ki©K2, 


H = 


= Ki© 


' iK2 


T : 


: P - 


^Pl©P2, 


M -- 


= PiE 


)iP2 


T : 


: iP - 


-^iPi ©iPa, 


iM 


= ^Pl 


©P2 



(5.18) 



As we already know, i^' is a compact subgroup, and exp(P) and exp(iP) define symmetric 
spaces with a definite metric (Riemannian spaces). In the same way, if is a non-compact 
subgroup, exp(M) and its dual space exp(iM) define symmetric spaces with an indefinite 
metric. These are pseudo-Riemannian symmetric coset spaces of a non-compact group 
by a maximal non-compact subgroup^. The original algebra G is thereby split into four 

components Ki, K2, Pi, P2, depending on their eigenvalues (++, H — , — h, ) under the 

two successive automorphisms cr, r. By applying all the possible a's and all the possible 
r's, or by replacing either a or r by the involutive automorphism cxr = rex, we obtain all 
the possible real forms of the symmetric spaces associated with the compact algebra G. 

Example: The complex algebra SO(3, C) has a root system of type Its compact 
real form is SO (3, R), and its only non-compact real form is SO(p, q; R) ~ SO(q, p; R) 
(p + q = 3), obtained by applying the involution cr=/p g (/g,p) to SO (3, R). In paragraph 
14.21 we constructed two Riemannian symmetric spaces associated with the algebra SO (3), 
the sphere S0{'i)/S0{2) and the double-sheeted hyperboloid 50(2, l)/S0{2). The Killing 
form has a definite opposite sign for the two spaces. 

The single-sheeted hyperboloid, described by the equation —x^ + + = 1 in R^, 
corresponds to the pseudo-Riemannian symmetric space S0{2, 1)/S'0(1, 1) associated with 
the same algebra. It is obtained by applying two consecutive involutive automorphisms 
<^=-^2,i5 'T=h,2 to the algebra G = SO(3,R). Like in eq. ()4.26p . J2,i and the Weyl unitary 
trick transforms G into G*. Let's now apply Ii^2 to G*: 

^Note that not all the theorems governing symmetric spaces corresponding to maximal compact sub- 
groups apply to the case at hand. A prime example is the decomposition involving radial coordinates in 
subsection 13.31 We will not discuss the symmetric spaces involving maximal non-compact subgroups in 
any detail in this paper. 
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= (KieK2)ei(Pi®P2) 

where in this example, Ki is empty. The spaces Ki, K2, Pi, P2 consist of the generators 
in G with the following combinations of eigenvalues under the two successive involutions, 
err: 



Ki : + + K2 : + - Pi : - + P2 
Thus we see that Ki is empty and the others are spanned by 



(5.20) 
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(5.21) 



The new symmetric space is obtained by doing the Weyl unitary trick on the split spaces 
(Ki©K2) and (PiSPa): 



it-^ 



H = Ki © iK2 = 1 I -it^ 



M = Pi ® iP2 = i 



t' 



(5.22) 



V -it'' -t' 



The second involution r (plus the Weyl trick) gives rise to a non-compact subgroup 
H—S0{1, 1) and to the symmetric space M ~ expM and its dual M* ~ exp(iM). The 
coset M ~ SO{2, l)/SO{l, 1) is represented by 
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expM 



y 



/ • \2 I 2 I 2 

[ix) + y + z 



(5.23) 



—tx —y 



The real forms of the simple Lie groups do not include all the possible Riemannian sym- 
metric coset spaces. For example, the compact Lie group G is itself such a space, and so is 
its dual /G (here the algebra G*^ = G* © iG* is the complex extension of all the real 
forms G*). By starting with a compact algebra G and applying to it all the combinations 
of the two involutive automorphisms a, r, we construct, in the way just described, all 
the remaining pseudo-Riemannian symmetric spaces associated to the corresponding root 
system. A complete list of these spaces can be found in Table 9.7 of reference |31j . 

Note that all the properties of the Lie algebra G (Killing form, rank, and so on) can be 
transferred to the vector subspaces P, iP [Hlj. The only difference is that the subspaces 
are not closed under commutation. 



In this section of the paper we have discussed symmetric spaces of positive, zero and 
negative curvature that will be relevant for matrix models of the circular, gaussian, and 
transfer matrix type, respectively. In Part II we will define and discuss various types of 
random matrix ensembles and their applications to various physical problems, and we will 
associate them to the corresponding symmetric spaces in Table ^ 



6 Operators on symmetric spaces 



The differential operator uniquely determined by the simplest Casimir operator on a sym- 
metric space (and especially its radial part) plays an important role both in mathematics 
and in the physical applications of symmetric spaces. Its eigenf unctions provide a complete 
basis for the expansion of an arbitrary square-integrable function on the symmetric space, 
and are therefore important in their own right. Their importance in the applications to be 
discussed in Part II is evident when considering that the radial part of the Laplace-Beltrami 
operator on an underlying symmetric space determines the dynamics of the transfer matrix 
eigenvalues of the Dorokhov-Mello-Pereyra-Kumar equation (DMPK equation for short) 
in the theoretical description of quantum wires, and maps onto the Hamiltonians of in- 
tegrable Calogero-Sutherland models. Here we will define some concepts related to the 
Laplace-Beltrami operator and discuss its eigenfunctions. 
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6.1 Casimir operators 



Let G be a semisimple rank-r Lie algebra. A Casimir operator (invariant operator) Ck 
[k = l,...,r) associated with the algebra G is a homogeneous polynomial operator that 
satisfies 

[Cfc,X,]=0 (6.1) 

for all Xi G G. The simplest (quadratic) Casimir operator associated to the adjoint 
representation of the algebra G is given by 

C = g^'XiX, (6.2) 

where g'^^ is the inverse of the metric tensor defined in ()3.15p and the generators Xj are in 
the adjoint representation. More generally, it can be defined for any representation p of G 
by 

C, = g^^p{X,)p{X,) (6.3) 

where g'^^ is the inverse of the metric ()3.16|) for the representation p (cf. subsection 13. 4j) . 
The Casimir operators lie in the enveloping algebra obtained by embedding G in the 
associative algebra defined by the relations 

X{YZ) = {XY)Z [X, Y]=XY - YX (6.4) 
(note that in general, XY makes no sense in the algebra G). 

The number of functionally independent Casimir operators is equal to the rank r of the 
group. Other Casimir operators can be formed by taking polynomials of the independent 
Casimir operators Ck {k = l,...,r). Since the Casimir operators commute with all the 
elements in G, they make up the center of the associative algebra ()6.4|) . 

Note that Casimir operators are defined for semisimple algebras, where the metric tensor 
has an inverse. This does not prevent one from finding operators that commute with 
all the generators of non-semisimple algebras. For example, for the euclidean group E3 
of rotations {Ji, J2, J3} and translations {Pi, P2, P3}, = J^PiPi and P- J = J^PiJi 
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commute with all the generators (cf. the comments in the paragraph following eq. ()5.5p ). 
Also the operators that commute with all the generators of a non-semisimple algebra are 
often referred to as Casimir operators. 

All the independent Casimir operators of the algebra G can be obtained as follows. Suppose 
p is an ra-dimensional representation of the rank-r Lie algebra G. The secular equation 
for the algebra G is defined as the eigenvalue equation 



CdimG \ n 

fp{X,) - A4 = Y.i->^r~'Mn = (6.5) 
4=1 / k=0 

where the (pk{t^) are functions of the real coordinates f. In general, they will not all be 
functionally independent (for example, v^o(^*) is a constant). There will be r functionally 
independent coefficients (pk{t^) multiplying the powers of —A jSI]. When writing down 
the secular equation, it is easiest to take a low-dimensional representation. By making 
the substitution f Xi in the functionally independent coefficients, they become the 
functionally independent Casimir operators of the algebra G: 



Mt') — > Q{X^) (6.6) 



Example: The generators Li, L2, of the SO(3) algebra were given explicitly in the 
adjoint representation in equations ()2.5|1 . ()2.6|1 in subsection 12.31 The secular equation for 
this algebra is then 



det (t ■ L - \h 



-X tV2 tV2 
-tV2 -A e/2 
-tV2 -A 



-\r + {-\)-/ = (6.7) 



The equation has one functionally independent coefficient, which is proportional to the 
trace of the matrix (t ■ L)^. It equals v^i(t) = |t^. The rank of S0{3) is 1 and the only 
Casimir operator is 



Ci^L"^ = Ll + Ll + Ll (6.8) 
obtained by the substitution f Li in fiit). The Casimir operator can also be obtained 
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from eq. ()6.2|) by using the metric gij = —^Sij for SO(3) given in the example in subsection 
13.41 We know from elementary quantum mechanics that 



[L^Ll] = [L^L2] = [L^L3]=0 (6.9) 

is an immediate consequence of the commutation relations, so we see that this operator 
indeed commutes with all the generators. Even though the commutation relations are 
not the same in polar coordinates or after a general coordinate transformation, ()6.9p will 
nevertheless be true. 

Example: SU (3) is a rank-2 group and therefore its characteristic equation will have two 
independent coefficients. If we denote a general SU(3) matrix (ajj) we get the characteristic 
equation 



an — A ai2 ctis 
det I 022 - A I = (-A)^ + (-A)^(aii + 022 + 033) 

031 032 (233 ~ A 

(6.10) 

+ (-A) (011022 + 022^33 + 033^11 - 012021 - 023^32 - 031013) 

+ (011(022033 - O23O32) + 012(023031 - O21O33) + 013(032021 - O31O22)) = 



The term proportional to (— A)^ vanishes, because the trace of any matrix in the SU(3) 
algebra is zero. The two independent coefficients are then v?2(oij) and (psicnj)- Substituting 
the values in terms of the coordinates f of the algebra J2i t^Ti for the aij (for example. 
Oil = + 7^^^! 012 = + it'^, etc.), we see that the expression for (p2it^) becomes 



Mtl = E(tT (6.11) 

i=l 

and therefore the substitution ()6.6p gives the first Casimir operator 

Ci = Hi + Hi + Y^iE^E^^ + E.^E^) =j2Tf (6.12) 

a. i=l 

as expected. Making the same substitution in ^p^{f) gives the second Casimir operator for 
S'f/(3), which has a more complicated form. 
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6.2 Laplace operators 



The Casimir operators can be expressed as differential operators in the local coordinates 
on the symmetric space. This is due to the fact that each infinitesimal generator Xa G G 
is a contravariant vector field on the group manifold. An element in the Lie algebra can 
be written 



X = Y.X-{^)d. ^ T.X'-ix)^ (6.13) 

where are local coordinates [HIIIISS] (for example, Li = (r x P)i = x^d^ — x^d2)- That 
the generators transform as lower index objects follows from the commutation relations. 

Example: As an example we take the group 5*0(3). Under a rotation R = R{t^,t'^,t^) = 
exp{J2t''Lk), the vector x = x^e^ G transforms as 



x-^x' = x'^e', (6.14) 
where the transformation laws for the components and the natural basis vectors are 

x'' = R'.x^, e'i = ejR\ (6.15) 
and R^^ = R^ . The one-parameter subgroups of 5*0(3) are rotations 

R{f) = exp(rL„), {n = 1, 2, 3) (6.16) 

(no summation) where L„ are 5*0(3) generators. It is easy to show using the commutation 
relations for L„ (given after eq. ()2.7p ) that under infinitesimal rotations the L„ transform 
like the lower index objects Cj: 

RLiR-^ = LjR\ (6.17) 

Expressed in local coordinates as differential operators, the Casimirs are called Laplace 
operators. In analogy with the Laplacian in R", 
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n Q2 

P^^A^E^ (6.18) 

which is is invariant under the group En of rigid motions (isometries) of R", the Laplace 
operators on (pseudo-)Riemannian manifolds are invariant under the group of isometries 
of the manifold. The isometry group of the symmetric space P ~ G/fi' is G, since G acts 
transitively on this space and preserves the metric, so the Laplace operators are invariant 
under the group operations g E G. 

The number of independent Laplace operators on a Riemannian symmetric coset space is 
equal to the rank of the space. As we defined in paragraph 13.21 the rank of a symmetric 
space is the maximal number of mutually commuting generators Hi in the subspace P (cf. 
also subsection I5.2|l . If Xa, Xp, ... G K and Xj, Xj, ... G P, it is also equal to the number 
of functionally independent solutions to the equation 

CdimP \ n 

J2 t'p{X,) -XlA= E(-A)"-V/(t') = (6.19) 
k=l ) 1=0 

where n is the dimension of the representation p and where now in the determinant we 
sum over all Xk G P. This is equivalent to setting the coordinates P for all the X^ G K 
equal to zero in the secular equation. In the example in the preceding paragraph, the rank 
of the symmetric space S0{3)/ S0{2) (the 2-sphere) is 1, which in this case is also the 
rank of the group 5*0(3). 

The Laplace-Beltrami operator on a symmetric space is the special second order Laplace 
operator defined (when acting on a function (0-form) /) as 



Asf = g'W,D,f = g^^{d,d, - r^/,)/ = -±^l-g^^^\l^f, g ^ detg,, (6.20) 



Here Z?j denotes the covariant derivative on the symmetric space and g^^ are the components 
of the inverse of the metric tensor. (The metric has an inverse because it is non-degenerate 
on a semisimple algebra and can be mapped over the entire symmetric space. For euclidean 
type spaces, we have the usual metric 6ij.) Di is defined in the usual way [SOI IS21 El! , 
example it acts on the components of a contravariant vector field in the following way: 
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(6.21) 



where Tj.- are Christoffel symbols (connection coefficients). The last term represents the 
change in x^ due to the curvature of the space. We remind the reader that on a Rieman- 
nian manifold, the Tl.- are expressible in terms of the metric tensor, hence the formula in 
eq. §M>- 

Example: Let's calculate the Laplace-Beltrami operator on the symmetric space 5*0(3) /S0{2) 
in polar coordinates using ()6.2()j) and the metric at the point [6, cj)) given in the second ex- 
ample of subsection 13.41 

^^^■=(o sin%)' ^^'=(o sin%) ^^-22) 

Substituting in the formula and computing derivatives we obtain the Laplace-Beltrami 
operator on the sphere of radius 1: 

Ab = dl + cot^ de + sin"^^ 9j (6.23) 

Of course this operator is exactly L^. We can check this by computing = yd^ — zdy, 
Ly = zdx — xdz, and = xdy — yd^ in spherical coordinates (setting r = 1) and then 
forming the operator + Ly + L^, remembering that all the operators have to act also on 
anything coming after the expression for each L^. We find that in spherical coordinates, 
expressed as a differential operator, is exactly the Laplace-Beltrami operator. 

In general, a Laplace-Beltrami operator can be split into a radial part A'^ and a transversal 
part. The radial part acts on geodesies orthogonal to some submanifold S, typically a 
sphere centered at the origin [HH] . 

Example: For the usual Laplace-Beltrami operator in expressed in spherical coordi- 
nates, 

Ab = dl + 2r-^dr + [pi + cot^ de + sin'^^ df) (6.24) 

the ffist two terms 



60 



(6.25) 



constitute the radial part with respect to a sphere centered at the origin and the expression 
in parenthesis multiphed by is the transversal part. The transversal part is equal to 
the projection of on the sphere of radius r and equals the Laplace-Beltrami operator 
on the sphere, given for r = 1 in eq. ()6.23p . This is a general result. For any Riemannian 
manifold V and an arbitrary submanifold the projection on S of the Laplace-Beltrami 
operator on V is the Laplace-Beltrami operator on S (see Helgason jSHj, Ch. II, paragraph 
3). 

The radial part of the Laplace-Beltrami operator on a symmetric space has the general 
form 



A's = 4tEt^"^^'^^ (j = 0,-,+) (6.26) 

where r' is the dimension of the maximal abelian subalgebra Hq in the tangent space P (the 
rank of the symmetric space) and J^^'^ is the Jacobian, to be given in equation (j6.3(J|) . of 
the transformation to radial coordinates. The sum goes over the labels of the independent 
radial coordinates defined in subsection 13.31 q = \ogh{x) = (g^, ...,q^ ) where h{x) is the 
exponential map of an element in the Cartan subalgebra and are canonical coordinates 
on Hq (in P they were denoted (g, a) = q ■ a; we will see in a moment that they are 
indeed given by q ■ a where a is a restricted root). 

The adjoint representation of a general element H in the maximal abelian subalgebra Hq 
follows from a form similar to eq. ()3.11|) (with or without a factor of i depending on whether 
we have a compact or non-compact space), but now the roots are in the restricted root 
lattice. For a non-compact space of type P* 



( 



\ / 



log/i = if = q ■ H 



q • a 



-q-r]/ \ 



q-^ ) 



(6.27) 
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Hence = q ■ a and 



/ 1 



h = e 



H 



(6.28) 



Example: For the simple rank-1 algebra corresponding to the compact group SU{2), the 
above formulas take the form (cf. eq. ()2.16j) ) 



/ 



H = 0Hi = 9 



I 1 



h = e 



-1 



(6.29) 



The radial coordinate is g = (g^) =6. 



There is a general theory for the radial parts of Laplace-Beltrami operators [211 • It is of 
interest to consider the radial part of the Laplace-Beltrami operator on a manifold V with 
respect to a submanifold W of V that is transversal to the orbit of an element w & W 
under the action of a subgroup of the isometry group of V. Of special interest to us is the 
case in which the manifold is a symmetric space G/K and the Lie subgroup is K. 

The Jacobian J^^'' = ^J\g\ (where g is the metric tensor at an arbitrary point of the 
symmetric space) of the transformation to radial coordinates takes the form 



^(-n?) = n.6/?+(sinh(g-))-'^ (6.30) 
^(+H9) = n.6/?+(sm(g"))"^" 



for the various types of symmetric spaces with zero, negative and positive curvature, re- 
spectively (see jnn], Ch. 1, par. 5). In these equations the products denoted na6-R+ '^^^^ 
all the positive roots of the restricted root lattice and is the multiplicity of the root a. 
The multiplicities rria were listed in Table CI 
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A remark on equation ()6.30|) is in order here. Strictly speaking, in the euchdean case we 
have not defined any restricted root lattice. The formula for the Jacobian J^^\q) for the 
zero-curvature space is understood as the infinitesimal version of the formula pertaining 
to the negative-curvature space. In the proof of the formula for the zero-curvature 
case we consider the Jacobian for a mapping from K/M x Hq onto P (where M is the 
centralizer in K of the Cartan subalgebra Hq), whereas in the proof of the formula for 
the negative-curvature symmetric space we consider the Jacobian for a mapping from 
K/M X e*^° onto a subset of the symmetric space e*^. can in both cases be interpreted 
as the dimension of the subspace of raising operators corresponding to the root a in the 
same algebra (this is, in the case of X~, the root multiplicity of the restricted root a; see 
|39j . Ch. I, paragraph 5). In the example illustrating the construction of restricted root 
systems in subsection 15.21 this space was spanned by the raising operators E'^ = Tg + iT2 

2 

and E'l = Tj + iTi for the root a = ^, whereas for the root a = 1 it was spanned by only 

one raising operator E[ = + iT^. The above can perhaps also be understood in terms 
of the limiting procedure discussed in paragraph 18.3.21 (cf. also equations ()6.42I6.43|) in 
section . The m^'s will in the following be referred to as "root multiplicities" also for 
the zero curvature spaces, keeping the above in mind. 

Example: On the hyperboloid with metric 





9^:i = n .^.1.2. , 9''=[ a ..-2n (6-31 



equations ()fi.2fipfi.3()|l give the radial part of for H^. The radial coordinate is 6 so we 
get, in agreement with ()6.30|) 

J^'^ = \f\9\= sinh^, A'^ = ii^b ^inh^ dg = {d^ + coth^ dg) (6.32) 

In the same way we can also easily derive the equation ()6.25|) and equation ()6.33|) below. 
In particular, comparing with ()6.32|) we immediately get the radial part of the Laplace- 
Beltrami operator acting on the two-sphere S"^ transversally to a one-sphere S"^ centered 
on the north pole: 



Mb = {dl + coW dg) (6.33) 
which is exactly the radial part appearing in equation (j6.23|) . 
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The eigenfunctions of the radial part of the Laplace-Beltrami operator on a symmetric 
space are called zonal spherical functions. In the applications of symmetric spaces to 
random matrix theory (for example in quantum transport), the properties of these zonal 
spherical functions are of central importance. For this reason we shall devote the following 
subsection (and the appendix at the end of this review) to a detailed discussion of their 
properties. 

The so-called DMPK operator will be discussed in section |H| The differential equation 
involving this operator describes the evolution of the distribution of the set of eigenvalues 
of the transfer matrix of a quantum wire with an increasing length of the wire. One of the 
most interesting applications of symmetric spaces in the random matrix theory of quantum 
transport lies in the identification of the DMPK operator with a simple transformation 
of the Laplace-Beltrami operator on the symmetric space defining the random matrix 
universality class. We will discuss this in more detail in Part II of this review (see subsection 

Em- 



6.3 Zonal spherical functions 

The properties of the so called zonal spherical functions are important for the research 
results to be discussed in Part II. Since there is a natural mapping from the Hamiltonians 
of integrable Calogero-Sutherland systems onto the Laplace-Beltrami operators of the 
underlying symmetric spaces, these eigenfunctions play an important role in the physics 
of integrable systems. But they are also relevant in transport problems in connection with 
the DMPK equation for a quantum wire. The known asymptotic expressions for these 
eigenfunctions allows one to solve this equation in general or in the asymptotic regime, 
because of the simple mapping from the DMPK evolution operator to the radial part of 
the Laplace-Beltrami operator. (For an example of their use see HT | H2].) 

When p is an irreducible representation of an algebra, the associated Casimir operators 
Ck,p are multiples of the identity operator [321 EZI (Schur's lemma). This means that it has 
eigenvalues and eigenfunctions. Since the Casimir operators (and consequently the Laplace 
operators) form a commutative algebra, they have common eigenfunctions. There exists an 
extensive theory regarding invariant differential operators and their eigenfunctions [321 • Of 
particular interest are the differential operators on a group G or on a symmetric space G/K 
that are left-invariant under the group G and right-invariant under a maximal compact 
subgroup K. 

Suppose the smooth complex-valued function (j)x{x) is an eigenfunction of such an invariant 
differential operator D on the symmetric space G/K: 
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D(t)x{x) = 7D(A)0A(a;) 



(6.34) 



Here the eigenfunction is labelled by the parameter A and 7d(A) is the eigenvalue. If in 
addition (j)\{kxk') = (j)\{x) (x G G/K, k & K) and (pxi^) = 1 (e =identity element), the 
function 0a is called spherical. A spherical function satisfies jSH] 

/ (j)x{xky) dk = (j)x{x)(f)x{y) (6.35) 

J K 

where dk is the normalized Haar measure on the subgroup K. We will see examples of this 
formula below. 

The common eigenfunctions of the Laplace operators on the symmetric space G/K are 
invariant under the subgroup K. They are termed zonal spherical functions. Because of 
the bi-invariance under K, these functions depend only on the radial coordinates h: 

0a(x) = Mh) (6.36) 

Example: Let's study for a moment the eigenfunctions of the Laplace operator on G/K = 
S0{3)/ S0{2). We know from quantum mechanics that the eigenfunctions of are the 
associated Legendre polynomials Pi(cos6'), and — /(/ + 1) is the eigenvalue under (our 
definition of L differs by a factor of i from the definition common in quantum mechanics) : 

VPi{cos9) = -l{l + l)Pi(cos^) (6.37) 

where cos6' is the 2;-coordinate of the point P = {x,y,z) on the sphere of radius 1. In 
spherical coordinates, P = (sin6' cos</), sin^sin^, cos^). As we can see, the eigenfunctions 
are functions of the radial coordinate 6 only. The subgroup that keeps the north pole fixed 
is = SO (2) and its algebra contains the operator = d^. Indeed, P;(cos6') is unchanged 
if the point P is rotated around the z-axis. 

Example: In terms of Euler angles, a general S'0(3)-rotation takes the form 
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R{a,P,^) = g{a)mh{^) 



( cosa — sina \ / cos/3 — sin/3 \ / C0S7 — sin7 
1 I I sin/3 cos/3 | 10 

y sm7 C0S7 
(6.38) 



1^ sina cosa ^ 



where g and h are rotations around the y axis by the angles a and 7 respectively, and 
is a rotation around the z axis by the angle (3. Under such a rotation, the north pole 
(0,0,1) goes into (—cosa cos/3 sin7 — sina C0S7, —sin/3 sin7, —sina cos/3 sin7 + cosq;cos7). 
This means that eq. ()6.35|) takes the form 

— / P/(—sinQ; cos/3 sin7 + cosa COS7) (i/3 = P;(cosa;)Pi(cos7) (6.39) 
zvr JO 

(To avoid confusion, note that the eigenvalue I in eq. ()6.37j) is not equal to A. In fact, 
A = / + l/2.) 

Example: For the symmetric space G/K = E2/SO{2) the spherical functions are the 
plane waves: 



V'(r) = e*'^" (6.40) 

where is a complex number, li g, h denote translations in the x-direction by a distance b, a 
respectively, and is a rotation around the origin of magnitude 0, then the transformation 
g{b)k{(f))h{a) moves the point x G by a distance y/a'^ + b'^ + 2abcos(f). Therefore we 
obtain from ()6.35|) 



^ ^ (^Y^a2 + 62 + 2a6cos0) dcp = ip{a)i;{b) (6.41) 
We introduce a parameter a into the the Jacobians ()6.30j) as in reference jH] , 

J^'Kq) = naei?+(«"^sinh(ag"))"^- (6.42) 
^^+^(g)=naeR+(«"'sin(ag-))-'^ 
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The parameter a corresponds to a radius. For example, for the sphere S0{3)/ S0{2) it is 
the radius of the 2-sphere. 

The various spherical functions corresponding to the spaces of positive, negative and zero 
curvature are then related to each other by the simple transformations [Oj 

0f (g) = lim„_o0^Hg) 

(6.43) 

There exist integral representations of spherical functions for the various types of spaces 
G/K ^ EH]- We will list an integral representation only of 0^ \q) below, recalling that 
formulas for the other types of spherical functions can be obtained by ()6.43|1 . If (f)\^ \x) is 
spherical and h is the spherical radial part of x, 

Here A is a complex-valued linear function on the maximal abelian subalgebra Hg of iP 
and p is the function defined below in eq. ()6.46|) . In eq. ()6.44p they act on the unique 
element H{kx) G Hq such that kx = ne^'^^^^k' in the Iwasawa decomposition introduced 
in subsection 13.31 It was shown by Harish-Chandra that two functions (f)^^\x) and 
(j)[^\x) are identical if and only if A = sz/, where s denotes a Weyl reflection. The Weyl 
group is the group of reflections in hyperplanes orthogonal to the roots and was defined in 
subsection 12. 5| eq. ()2.34|1 . and discussed further in subsection 12.61 

Equation ()6.44|) may seem a bit cryptic, but it becomes much more clear if one uses the 
explicit expression for the Iwasawa decomposition. It essentially becomes the integral over 
the product of all the possible lower principal minors (raised to suitable powers) of the 
matrix ke°'^k~^ (where a is the free parameter introduced above). The explicit expression 
can be found in [9J. This integral representation is of great importance, because in some 
cases it can be exactly integrated leading to explicit expressions for the zonal spherical 
functions. This point will be discussed in more detail in the appendix at the end of this 
review. 

The eigenvalues of the radial part of the Laplace-Beltrami operator corresponding to the 
eigenfunctions on zero, negative and positive curvature symmetric spaces are given by the 
following equations (see |^ and jSlI, Ch. IV, par. 5): 
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A'b0^^ = (-S-P^)0^^ (6.45) 



where p is the function defined by 



Example: Take the symmetric space S0{?))/ S0{2). From Tabled we see that this space 
has p — g = 2 — 1 = 1 short root of length 1. Then 




(6.47) 



and setting a = 1, the eigenvalue is — + 1/4 = — /(/ + 1). 

Remarkably enough, for a few classes of symmetric spaces explicit expressions for the zonal 
spherical functions can be obtained. We will make use of this result in section 1^.4. II when 
discussing the exact solution of the DMPK equation in the (3 = 2 case. These solutions 
play an important role in several branches of physics. For this reason we decided to discuss 
them in some detail in the appendix of this review. 



6.4 The analog of Fourier transforms on symmetric spaces 

Much of the material presented in this subsection is taken from the book by Wu-Ki Tung 

A continuous smooth (C°°) spherical function / is said to be elementary if it is an eigen- 
function of any differential operator that is invariant under left translations by G and right 
translations by K. Thus the eigenfunctions of the Laplace operators are elementary. The 
elementary spherical functions are related to irreducible representation functions for the 
group G. The irreducible representation functions are the matrix elements of the group 
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elements g in the matrix representation p. Let's clarify this statement by an example 
well-known from quantum mechanics. 



Example: The angular momentum basis for 5*0(3) is defined by 



L2|/m >= /(/ + l)|/m > 

Lal/m >= m|/m > (6.48) 
L±\lm >= Jl{l + 1) - m{m ± l)\lm > 



where / labels the representation. The irreducible representation functions are, in the 
angular momentum basis |Zm >, the matrix elements Z^'(i?)"^ such that 



R\lm >= \lm! > D\R)1^ (6.49) 

where R = exp(t ■ L) is a general 5*0(3) rotation. It is known that for R G 50(3), if a, f3, 
7 are the Euler angles of the rotation R = R{a, (3, 7), these matrix elements take the form 



D^{R)"i = D\a,i3,^)"'J^ = e-^""^'(i^(/3)"^e-^^"^, d^{f3)"i =< lm'\e-'^^^\lm > (6.50) 



The associated Legendre functions P™(cos^) and the special functions Yi^{9,<p) called 
spherical harmonics are essentially this kind of matrix elements: 



The irreducible representation functions /^^(i?)"^ satisfy orthogonality and completeness 
relations. In fact, they form a complete basis in the space of square integrable functions de- 
fined on the group manifold. This is the Peter- Weyl theorem. Prom here the corresponding 
theorems follow for the special functions of mathematical physics. 

Example: For 50(3) the orthonormality condition reads 

{21 + 1)J drDliRnD^'iRf^, = SfS^'S^ dI{R)Z ^ [D^RTJ* (6-52) 
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where R = R{a,P,'j) is an S0{3) rotation expressed in Euler angles and dr is the in- 
variant group integration measure normahzed to unity, dr = dad{co^(3)d'-) / '^ti'^ . That the 
irreducible representation functions form a complete basis for the square-integrable func- 
tions on the SO (3) group manifold can be expressed as 



f{R) = T.fimD\R)^ (6.53) 

Imn 

where f{R) is square-integrable. Using ()(i.52|l we obtain 

= {21 + 1)J dTDj{R)\f{R) (6.54) 

If R{a, P, 7) = R{(p, 9, 0) we get the special case of the spherical harmonics on the unit 
sphere (setting J An/ {21 + 1)/°^ = fi^): 



f{e,<P) = J2fimYim{0A) (6.55) 

Im 



fim = J f{e, mUO, 4>)d{cose)d<P (6.56) 
and further, for R{a,P,'-y) = R{0,9,0) we get the completeness relation for the associated 



Legendre poljTiomials Pi{cos6) = Yi^^ An / {21 + 1) 

f{e) = Y.m{cose) (6.57) 

{21 + 1) 



fi = / f{e)p;{cose)d{cose) (6.58) 



where // = /;q. These are analogous to Fourier transforms. In the above example, we 
considered a symmetric space with positive curvature. For a space with zero or negative 
curvature, we have an integral instead of a sum in ()6.57|) : 



f{q) = I /(A)0?(g)d/i(A) <xw' [ f{X)f^dX (6.59) 



|c(A)P 
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/(A) = / fiq) J^'\q)dq (6.60) 



where q are canonical radial coordinates. The integration measure dii{X) was determined 
by Harish-Chandra [SI to be well-defined and proportional to w'^\c{X)\~'^dX, where c(A), 
to be discussed in detail below, is a known function whose inverse is analytic (in this 
context, see also EH]) and w is the order of the Weyl group (the number of distinct 
Weyl refiections). J^^\q)dq is the invariant measure on the space of radial coordinates. In 
equation ()6.59j) the arbitrary square-integrable function /(g) is expressed in terms of the 
complete set of basis functions (j)x\q)- 

One can show that the dimension of the space of eigenf unctions of A'^ is less than 

or equal to w. It is a remarkable fact that the eigenf unctions of the radial part of the 
Laplace-Beltrami operator A'^ have the property of being eigenf unctions of the radial part 
of any left-invariant differential operator on the symmetric space as well ([12], Ch. IV). 

In the asymptotic expansions that we will discuss in subsection 19.41 of Part II of this 
review in the context of quantum transport, a crucial role is played by the function c(A) in 
eq. ()6.59p . It encodes all the information relating the transport problem to an underlying 
symmetric space. The explicit form of c(A) is 



c(A) = n C"(A) (6.61) 

oe-R+ 



with 



r(aV2) 

where F denotes the Euler gamma function, a is a generic root belonging to the restricted 
root lattice of the symmetric space, rria denotes its multiplicity and the product is restricted 
to the sublattice of positive restricted roots only. We shall see an explicit example of 
these functions in equations ()9.23p and ()9.24|) for the three symmetric spaces which are 
relevant for the random matrix description of quantum wires. 

An important feature of the zonal spherical functions is that they satisfy, for large values 
of \h\, the following asymptotic expression for spaces of zero and negative curvature PHH]: 
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(6.63) 



where h = e is a spherical coordinate, H is an element of the maximal abelian subalgebra, 
A is a complex-valued linear function on the maximal abelian subalgebra, and the function 
p was defined in eq. ()6.46j) . 

This expansion will play a major role later in this review. We will refer to it in the following 
as the "Harish-Chandra asymptotic expansion" . 



7 Integrable models related to root systems 

As mentioned in the introduction, an important role in our analysis is played by the class 
of integrable models known as Calogero-Sutherland (CS) models, which turn out to be 
deeply related to the theory of symmetric spaces. These models describe n particles in one 
dimension, identified by their coordinates g^, and interacting (at least in the simplest 
version of the models) through a pair potential f (g* — q^). The Hamiltonian of such a 
system is given by 



where the coordinate q is q = (g^,...,g"), are the particle momenta, and the 

particle mass is set to unity. In eq. ()7.1|1 is the subsystem of positive roots of the 
root system R = {a^, ....a"} related to a specific simple Lie algebra or symmetric space, 
and n is the dimension of the maximal abelian subalgebra Hq and of its dual space Hq. 
The components of the positive root a = a'^ G -R"*" are a^, a^. The number of positive 
roots is Z//2, where v is the total number of roots. In general, the coupling constants Qa 
are the same for equivalent roots, namely those that are connected with each other by 
transformations of the Weyl group W of the root system (see subsections 12 .51 IT^ and . 

Several realizations of the potential f (g") have been studied in the literature (for a review 
see ref. P): 



Pi = 




(7.1) 
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viiiO 
viiiiO 

vviiO 



sinh-2(0 

no 



(7.2) 



Here V{q) denotes the Weierstrass P-function, to be defined in eq. ()10.2U|) . We will 
mainly be interested in the first three realizations. The potential expressed in terms of the 
Weierstrass P-function will be discussed in subsection 110. 3[ while we will not deal with 
the last two cases which we reported here only for completeness. The reader is referred to 
[Oj for a discussion of these two potentials. 

The most relevant features of the Calogero-Sutherland models are: 

• Under rather general conditions (see |9j for a detailed discussion) they are completely 
integrable, in the sense that they possess n commuting integrals of motion. 

• For particular values of the coupling constants (i.e. for those related to the root mul- 
tiplicities by eq. ()7.6p which we will discuss below) the Calogero-Sutherland Hamil- 
tonians can be mapped onto the radial parts of the Laplace-Beltrami operators on 
suitably chosen symmetric spaces. As we will see, these spaces have negative curva- 
ture for the sinh-type models and positive curvature for the sin-type models. 



7.1 The root lattice structure of the CS models 

In the original formulation of the CS model, the interaction among the particles was simply 
pairwise |H1. Only later it was realized that this particular choice was the signature of an 
underlying structure, namely the root lattices of Lie algebras of type An (see the example 
below). Also, the model could be extended to any root lattice canonically associated to 
a simple Lie algebra or symmetric space, keeping its relevant properties: complete inte- 
grability and mapping to the radial part of a Laplace-Beltrami operator for special values 
of the couplings 9J. This corresponds exactly to the choice of Hamiltonian introduced in 



eq. (ED). 
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Let us look at two examples which may clarify this construction. In these examples we 
use the root lattices An and C„ and the potential of type II. Remember that the coupling 
constants ga are the same for equivalent roots, thus we expect a single coupling constant 
in the first example and only two copuling constants in the second example. In order to 
fix the notation, let us denote with {ei, Cn} a canonical basis in the space R". 



An'. As mentioned in subsection 12.71 the An root system is contained in the hyperplane 
in R""*"^ with equation xi + X2 + ... + Xn+i = 1- The root system R is given hj R = 
{a^, a"} = {ei — ej, i ^ j}. In this case VT, the Weyl group, is the permutation 
group of the set {cj}. The corresponding CS Hamiltonian is in this case 

^ = "2g%^ + §sinh^(g^-g.) ^^'^^ 

The arguments of the sinh-function are 5'^ = ■ = (q'^, (t^ ■ (cj — Cj) = — 
[i < j) where a'^ is a positive root of the root lattice R. As anticipated we have only 
pairwise interaction and a single coupling constant appears in the model. This is the 
model originally considered in [Sj. 

Cn'. The Cn root system is R = {±2ej, ie^ ± ej,i ^ j}. The Weyl group is the product 
of the permutation group and the group of transformations that changes the sign of 
the vectors {cj}. The corresponding Hamiltonian is: 




^ = E 7^ + E -ifcrv + E ..^.J: + J ■ , (7.4) 



As anticipated, we have different coupling constants for the ordinary roots (corre- 
sponding to pairwise interaction) and for the long roots. In the following sections 
we will come back to this particular choice of Hamiltonian which turns out to be of 
great importance in the random matrix description of quantum wires. 



7.2 Mapping to symmetric spaces 



The most interesting property of these Hamiltonians, which is a direct consequence of the 
underlying structure of the symmetric space, is that for the potentials of type I, II and 
III there exists an exact transformation of the Hamiltonian Ti. into the radial part of the 
Laplace-Beltrami operator on the symmetric space corresponding to the root lattice R of 
Ti. In particular the target spaces have negative curvature for the sinh-type models and 
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positive curvature for the sin-type models. Denote the radial part of the Laplace-Beltrami 
operator by A'^; the transformation is then given by (see appendix D of ref. jH] for a proof) 



7^ = e(g)^(A'B±p')r'(g) (+ for II, - for III, p = for I) (7.5) 
if and only if the coupling constants Qa in Ti take the following root values 



2 _ ma{ma + 2m2a - 2)|ap 

9a g I ' •'-'J 

In equation ()7.6p is the multiplicity of the root a and |a| its length, and p in eq. ()7.5|) 
is the vector defined in (16.4611 



For the ordinary roots eq. (|7.6|) simplifies to 



|2 



2 _ ma{ma - 2)\a\ 

9a- g U- 



The function ^{q) in ()7.5|) is related to the Jacobian of the transformation to radial coor- 
dinates (cf. equation ()6.30p ). It is given by 



m = naei?+[sinh(g")]""/2 II (7.9) 
n„6ij+[sin(g-)]™'^/2 III 

for quantum systems with potentials of type I, II, and III, respectively. The radial part of 
the Laplace-Beltrami operator, A'^, on the symmetric space can be expressed in terms of 
^ (g) as follows (see equation ()6.26p in subsection IHZ 
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and the symmetric space has zero, negative, and positive curvature for the cases I, II, and 
III, respectively. 

At this point a number of results can be obtained for the corresponding quantum systems 
merely by using the theory of symmetric spaces. A detailed collection of results pertaining 
to spectra, wave functions, and integral representations of wave functions can be found in 
the original article 0, and we will not further elaborate on them here. The crucial point 
in this relation for our analysis is another one. From eq. ()7.8p we see that whenever the 
ordinary root multiplicity is = 2, the interaction between particles in the Hamiltonian 
vanishes. The degrees of freedom decouple and the problem can be solved explicitly. Such 
a decoupling could hardly have been inferred by looking at the original form of the radial 
part of the Laplace-Beltrami operator, equation ()6.26|) . This decoupling is also the reason 
why in these cases explicit expressions for the zonal spherical functions exist (see paragraph 
16.31 and, for a more detailed discussion, the appendix at the end of this review). 
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Part II 



In the second part of this review we will use the background material presented in Part I 
to make evident the close relationships between symmetric spaces and random matrix 
ensembles. To this end, we discuss in detail the identifications that have to be made 
between matrix ensembles and symmetric space characteristics in section |H1 We show 
that the integration manifolds of random matrix theories with physical applications in the 
description of complex nuclei, gauge field theories and mesoscopic systems, can be exactly 
identified with the irreducible symmetric spaces based on non-exceptional groups that 
were classified by Cartan, and their euclidean counterparts. We identify the random matrix 
eigenvalues with the spherical radial coordinates frequently used in the theory of Lie groups, 
and show that the Jacobians determining the probability distribution functions of random 
matrix eigenvalues are exactly the Jacobians appearing in equation ()6.42|) determined by 
the metric g on the symmetric space: J = \f\g\- The structure of the restricted root 
lattice associated to the symmetric space determines this Jacobian completely through the 
multiplicities of its roots. In connection with the discussion of transfer matrix ensembles, 
we define the Dorokhov-Mello-Pereyra-Kumar (DMPK) equation, and in subsection 18.51 
we include a general discussion of Fokker-Planck equations associated to random matrix 
ensembles. We also discuss the so-called Coulomb gas analogy used in random matrix 
theory. We show that the DMPK operator connected to any random matrix ensemble is 
simply related to the radial part of the Laplace-Beltrami operator (defined in section Ej) 
on the appropriate symmetric space. We conclude by summarizing the results in Table |21 

We devote section El to some consequences of the identifications between random matrix 
theories and symmetric spaces. In subsection 19.11 we discuss the new classification of 
disordered systems that the equivalence gives rise to, as a natural consequence of the 
Cartan classification of symmetric spaces. In subsections I9.2H9.3I we will study how the 
properties of the symmetric spaces are refiected in the orthogonal polynomials associated 
to a random matrix partition function, and see how the refiection or translation symmetries 
of the Jacobians are direct consequences of the correponding properties of the restricted 
root lattices. 

An introduction to the relationship between the restricted root lattices of symmetric spaces 
and integrable Calogero-Sutherland models describing interacting many-particle systems 
in one dimension was given in section [7| Olshanetsky and Perelomov |^ showed that 
the dynamics of these systems are related to free diffusion on a symmetric space. This 
relationship is due to the fact that the Hamiltonians of integrable Calogero-Sutherland 
models map onto the radial parts of the Laplace-Beltrami operators of the underlying 
symmetric spaces, as explained in section [7| of Part I. 

Since the Dorokhov-Mello-Pereyra-Kumar equation for a disordered conductor can be 
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mapped onto a Schrodinger-like equation in imaginary time (021; see also jH]) featuring 
a Calogero-Sutherland type Hamiltonian with root values of the coupling constants, in- 
formation on the solutions of the DMPK equation jlHl 1^ can be extracted from known 
properties of the zonal spherical functions of the underlying symmetric space. As we have 
seen in section IHl these are the eigenfunctions of the radial part of the Laplace-Beltrami 
operator. Some results relating to this will be described in detail in subsection 19.41 as 
well as an application of the theory of symmetric spaces in a quantum transport problem 
involving the magnetoconductance. 

Finally in section El we discuss some possible extensions of the present results and some 
non-standard applications of the symmetric space formalism. We go beyond the Cartan 
classification and discuss non-Cartan parametrization of symmetric spaces in subsection 
110.11 In subsection 110.21 we discuss generalizations of the DMPK equation and the al- 
ternative clustered solutions to the DMPK equation that are a consequence of the exact 
integrability of Calogero-Sutherland models. Finally, in the last subsect ion II . 31 we discuss 
the Weierstrass "P-function and show that it describes three types of Calogero-Sutherland 
potential in its various limits, corresponding to symmetric spaces of different curvature. 
These limiting potentials correspond in the Calogero-Sutherland model to particles inter- 
acting on a circle, hyperbola, or line, respectively, reflecting the corresponding triplet of 
symmetric spaces underlying the respective Calogero-Sutherland models. 



8 Random matrix theories and symmetric spaces 

8.1 Introduction to the theory of random matrices 
8.1.1 What is random matrix theory? 

Random matrix theory has evolved into a rich and versatile field with applications in 
several branches of physics and mathematics. In the theory of random matrices one studies 
the statistical properties of the eigenvalues of large matrices with randomly distributed 
elements. Historically, large random matrices were first employed by Wigner j37j and Dyson 
jini to describe the energy levels in complex nuclei, where they modelled the Hamiltonian 
of the system. Disregarding the continuous part of the spectrum, one may represent the 
Hamiltonian of such a system quantum mechanically by a large hermitean matrix acting 
on a finite-dimensional Hilbert space. In principle, the eigenvalues of the Hamiltonian 
then give us the energy levels of the nucleus, and the eigenvectors give us the eigenstates. 

However, for a system involving hundreds of nucleons (or, in a small metal sample, hundreds 
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of electrons) we do not know the form of the Hamiltonian, and even if we did, the huge 
number of degrees of freedom involved would prevent us from solving the Schrodinger 
equation exactly. Choosing instead to represent the Hamiltonians of an ensemble of such 
systems by an ensemble of large hermitean matrices, whose elements are random variables 
with some given distribution, we obtain a statistical theory for its eigenvalues Ai,...,AAr, 
where is a large number equal to the size of the random matrices If the random 
matrices in the chosen ensemble have the same global symmetry properties as the actual 
Hamiltonian, this statistical theory describes also the eigenvalues of the Hamiltonian, to 
the extent that they depend only on symmetry, i.e. they are universal. Even systems 
with just a few degrees of freedom lend themselves to such a description. In general, a 
description in terms of random matrices is possible whenever we are dealing with a chaotic 
or disordered system (see |2H] for a discussion). 

Such random matrices were studied by Dyson in a series of papers He showed that 
they fall into one of three classes named the orthogonal, unitary and symplectic ensemble, 
depending on whether or not the Hamiltonian possesses time-reversal invariance and ro- 
tational invariance. This classification was reviewed in Ch. 2 of the book by Mehta [TT] . 
The corresponding random matrices are real symmetric, complex hermitean, or self-dual 
quaternion, respectively (see below). 

8.1.2 Some of the applications of random matrix theory 

In typical applications, large matrices with given symmetry properties and randomly dis- 
tributed elements substitute a physical quantum operator. This operator could be a Hamil- 
tonian like in the example of nuclear energy levels or in quantum chaotic scattering (in 
the latter case an effective Hamiltonian appears in the scattering matrix), or it could be 
for instance a scattering matrix, transfer matrix, or Dirac operator. Then one studies the 
statistics of its eigenvalue spectrum and extracts the universal behavior. 

An early application of the Wigner-Dyson ensembles was in the theory of scattering in 
chaotic quantum systems. With the help of random matrices one is able to study the 
statistics of resonance poles and scattering phase shifts in nuclei, atoms, and molecules, as 
well as microwave cavities and ballistic systems (for a self-contained review see !;48.5). 

A major area where random matrices are employed is in the description of the infrared limit 
of gauge theories, notably QCD. In these applications an integration over an appropriate 
random matrix ensemble replaces the integration over gauge field configurations in the 
partition function. This is achieved by substituting a suitable random matrix for the 
Dirac operator appearing in the fermion determinant. Because the magnitude of the quark 

^°For symplectic ensembles the eigenvalues are two-fold degenerate. 
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condensate is proportional to the spectral density of the Dirac operator at the origin of the 
spectrum, random matrices are useful in this context in obtaining insights concerning the 
spontaneous breaking of chiral symmetry. This is one of the most fundamental phenomena 
in QCD, since it determines the hadronic mass spectrum. As it is a phenomenon that takes 
place at very low energy, it is inaccessible to perturbation theory and therefore one of the 
most difficult to study from a theoretical point of view. The random matrix approach 
brings considerable simplification and has indeed contributed to the understanding of this 
phenomenon. The corresponding random matrix ensembles are called chiral ensembles, 
because of the chiral symmetry of the Dirac operator. We will discuss chiral random 
matrix theories in more detail in subsection 18. 3. HI 

Random matrix theories are also largely and successfully applied in the theoretical descrip- 
tion of mesoscopic systems, where they may be used to model the properties of scattering 
and transfer matrices. A mesoscopic conductor is a micrometer-sized metal grain or wire 
at low temperature with randomly distributed impurities that act as scattering centers for 
electrons originating from current leads (alternatively, random scattering takes place at 
the edges and depends on the shape of the conductor). Examples include disordered wires, 
quantum dots, and normal metal-superconductor heterostructures. 

The theory of quantum transport is concerned with the statistics of the transmission eigen- 
values in such systems. The transmission eigenvalues are related to the transfer matrix 
in a way which will be discussed in detail in paragraph 18.3.41 These eigenvalues directly 
determine the physical conductance. We will discuss quantum wires in some detail below. 
We will also briefly discuss applications in normal metal-superconductor heterostructures 
as well as in ballistic quantum dots. The ensembles used in these applications are related 
to yet other symmetric spaces than those pertinent to the ensembles used in the study of 
the usual mesoscopic systems. The reason is, of course, that the symmetries of the physical 
operator whose eigenvalues we study are different. 

The same chiral ensembles used in the description of chiral gauge theories are also realized 
in the Hamiltonians of random flux and random hopping problems |49, ,50J in the theory 
of quantum transport. 

Some important applications of hermitean matrix models which will not be treated here, 
are in the description of random surfaces in the field of quantum gravity (for an excellent 
introduction see ^H]); where a N~'^ expansion amounts to a genus expansion of the random 
surface, and in string theory. 

The wide range of applications is one of the most fascinating aspects of random matrix 
theory. Notably, the classical random matrix theories have no adjustable parameters. 
The reason for their success in accurately characterizing such a wide range of systems 
(nuclei, disordered metals, chaotic systems) lies in the universal behavior of the eigenvalue 
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correlators in an appropriate double scaling limit. This is the subject of the following 
paragraph. 



8.1.3 Why are random matrix models successful? 

The statistical properties of a sequence of apparently random rumbers like nuclear energy 
levels may be either of the local or of the global type. The local and global spectral 
characteristics are completely disconnected. Systems with identical global characteristics 
may have different local characteristics, and vice versa. 

A typical example of a global, or macroscopic, property is the mean level (eigenvalue) 
density, or one-point correlation function, in the limit of large N . It is often denoted p(A). 
For the gaussian Wigner-Dyson ensembles the mean level density approaches a semicircle in 
the large limit. In general its shape depends on the random matrix potential determining 
the probability distribution of individual matrix elements. A global characteristic spectral 
property does not change appreciably on the scale of a few level spacings. An example of 
a local spectral property is the spacing between two successive levels. It fluctuates from 
level to level. 

One amazing feature of random matrix models is that some local spectral characteristics 
are independent of the distribution of individual matrix elements in the limit of large 
A^, if we simultaneously rescale the eigenvalues in an appropriate way. This is called a 
double scaling limit. In this limit local spectral characteristics are determined solely by 
the global symmetries of the random matrices. This so called microscopic universality 
is manifest only in the appropriate double scaling limit, where the average level spacing 
becomes much smaller than the scale on which it varies. 

The n-point correlation function, or joint probability density for n eigenvalues p„(Ai, A„) 
{1 < n < N) is central in the theory of random matrices and is the most important example 
of such a microscopically universal quantity. It is defined as follows: Pn{^i, \n)d\i...d\n 
is the probability of finding one eigenvalue in each of the intervals [Aj, Aj + dXi] (1 < 
i < n). This means that if one rescales the eigenvalues by the average level spacing^^ 
and simultaneously takes the limit N ^ oo, the n-point function p„ for each ensemble 
becomes a universal function characteristic of the ensemble. There are also other spectral 
characteristics that depend on symmetry only. The number variance (i.e., the variance 
related to the number of eigenvalues in a certain interval), the so-called spectral rigidity A3 
(this quantity was defined for example in fHI^), the distribution of the first eigenvalue, 
and the distribution of spacings between adjacent levels are examples of such quantities. 

^^In general one can consider also other scales. 
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In a hermitean matrix model the eigenvalues are real. The eigenvalue density (one-point 
function) of such a model is typically zero outside some interval or intervals around the 
origin. As we learnt above, its macroscopic shape depends on the probability distribution, 
while locally its microscopic shape is universal. Spectrum edge universality, i.e. universality 
of scaled correlation functions near the edge of the eigenvalue spectrum, is known to differ 
from universality in the bulk of the spectrum. One normally distinguishes between three 
microscopic scaling regimes: the bulk, the hard edge (vicinity of the origin) and the soft 
edge (tail) of the spectrum. The rescaling of the eigenvalues must be chosen appropriately 
in each regime, and in practice amounts to a rescaling by an appropriate integer or fractal 
power of N [221 ISHl ISH IHS] , which is proportional to the inverse local density of eigenvalues 
and thus proportional to the average local level spacing. This procedure in turn amounts 
to magnifying a local region in such a way that the behavior of the eigenvalue correlators 
is universal. 

In case of a (multi) critical hermitean matrix model, the macroscopic spectral density de- 
velops extra (multiple) zeros on the limit between a one-cut (single-interval) support and 
a multi-band support. This is achieved by fine-tuning couplings in the random matrix 
potential determining the probability ditribution. In this case the scaling is more subtle 
and the calculations more involved due to subleading terms in the large- expansion [HH] . 

Several attempts at theoretically demonstrate universality can be found in references 
E31 IMl mi EHI- It is also supported by a wealth of numerical jSZI evidence. In 
the universality of correlation functions near the origin of the eigenvalue spectrum 
of chiral unitary and unitary ensembles of random matrices in the microscopic limit was 
investigated. Universality was shown to follow from the fact that the recursion relation 
for the orthogonal polynomials determining the kernel from which all correlators are de- 
rived, reduces to a universal differential equation (in this case for Bessel functions) in the 
appropriate microscopic limit. 



8.2 The basics of matrix models 

In this subsection we will give a few basic definitions and results regarding the simplest 
kinds of hermitean matrix models. Similar results apply for all the matrix models we will 
study in this paper, and the purpose here is to introduce the reader not at all familiar with 
matrix models into the subject. 

A general hermitean matrix model is defined by a partition function 
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where is a square N x N hermitean matrix with randomly distributed elements Sij 
and dS is a Haar measure. For example, for the gaussian unitary ensemble consisting 
of complex hermitean matrices (this ensemble describes a system with no time-reversal 
invariance) the matrix elements are complex numbers, Sij = S^j + iS^j (except the ones on 
the diagonal that must be real). dS is then given by 

dS = l[dSl^l[dSl, (8.2) 

i<j k<l 

and the probability P{S)dS that a system described by the gaussian unitary ensemble will 
belong to the volume element dS is invariant under the automorphism 

S U-^SU (8.3) 

of the ensemble to itself, where U is any unitary N x N matrix^^. This amounts to a 
change of basis. The invariance of P{S)dS under such automorphisms restricts P{S) to 
depend on the traces of the first powers of S. For the gaussian matrix models reviewed 
for example in the book by Mehta, there was another requirement, namely that the matrix 
elements should all be statistically independent. This excludes everything except the traces 
of the first two powers, and these may occur only in an exponential ^T]. However, this 
requirement is not physically motivated and was subsequently relaxed as other types of 
matrix models were introduced. A probability distribution of the form 

P{S) oc e-'^*'^^^^) (8.4) 

is often used. Here c is a constant and V{S) is a matrix potential, typically a polynomial 
with a finite number of terms. 

By doing an appropriate similarity transformation ()8.3|1 on the ensemble of random matri- 
ces, the Haar measure dS and potential V{S) can be expressed in terms of the eigenvectors 
and eigenvalues of the matrix S. This amounts to (block) diagonalizing S. At the same 
time the Haar measure is factorized into a part that depends only on the eigenvectors (the 
"angular" part) and a part that is a function of the eigenvalues {Aj} of S only. Assume 
for simplicity that the random matrix potential is given by 

^^If S is time-reversal invariant U will be either a real orthogonal matrix or a symplectic matrix, 
depending on whether rotational symmetry is also present. In this case we are dealing with the orthogonal 
or symplectic ensemble. 
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cV{S) = (8.5) 

This kind of potential is called gaussian, for obvious reasons. The Jacobian of the similarity 
transformation depends on {Aj} only and is given by [TT] 



^({M)~n(^^-^.)' (8-6) 

i<j 

for the gaussian unitary ensemble. For the other gaussian ensembles one has a similar 
Jacobian (see eq. ()8.29j) below). The part of P{S)dS that depends on the eigenvalues then 
takes the form 



. ^^.v,j;^.v,...^.vjv - -j; - ^=1^ c/Ai...c/A^ (8.7) 

i<j 

whereas the "angular" degrees of freedom do not appear in the integrand at all and can be 
integrated out to give a constant in front of the integral. This model is easily solvable. This 
is the strength of the random matrix description of disordered systems. The eigenvalue 
correlators (fc-point functions) are defined as 



P/c(Ai,...,A,) = ^^^^— n^^^.'^({^i'-'^^}) (8-8) 

and can be calculated exactly by rewriting the Jacobian as a product of Vandermonde 
determinants of a set of polynomials orthogonal with respect to the measure e"'^^'-'**-*. This 
procedure was reviewed in [TT]. For the gaussian unitary ensemble these polynomials are 
Hermite polynomials Hn{X) and satisfy 



KMHn{\)e-^'d\ = hJmn (8.9) 

where is a normalization. In practice one defines a kernel K]\i{Xi,Xj) that turns out 
to be universal and whose determinant gives all the correlation functions pjj. Defining 
so-called oscillator wave functions that for the gaussian unitary ensemble take the form 
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yp,(A) = ^Hj{X)e-^ (8.10) 



this kernel is given by 

Kn{K, Xj) = E MX^)MXJ) (8.11) 

fc=0 

It can be evaluated explicitly using the Christoffel-Darboux formula for the sum of products 
of orthogonal polynomials. The /c-point function is then 

Pfc(Ai, Afc) = det[i^'jv(Ai, Xj)]i<i,j<k (8.12) 

Different authors use different conventions for the normalization constants involved in these 
formulas. 

So far we have not said anything about the integration manifold in equation ()8.H) . As we 
have seen, we start out with some given class G of matrices S, and the integral is a priori 
over this class. However, by doing the similarity transformation in eq. ()8.3|) . we can perform 
the integration over the subgroup K to which the matrix U belongs without further ado. 
We then end up with an integral over a smaller manifold G/K, the manifold of the random 
matrix eigenvalues, and as we will see in the following subsection, this manifold turns out 
to be a symmetric space for all the commonly used hermitean random matrix ensembles. 
This gives rise to the deep and useful connection between random matrix theories and 
symmetric spaces that is the subject of this review. 

Let us now look more in detail at some explicit examples of random matrix theories. In 
the following subsections we will identify in turn their integration manifolds as symmetric 
spaces, their eigenvalues as radial coordinates on the respective symmetric spaces, and 
their Dyson and boundary indices (to be defined) as determined by the multiplicities of 
the corresponding restricted root lattices. As we will see, the Jacobian of the similarity 
transformation to the space of random matrix eigenvalues and eigenvectors is determined 
explicitly by these root multiplicities and the curvature of the symmetric space manifold. 
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8.3 Identification of the random matrix integration manifolds 

In this subsection we introduce and explicitly define some random matrix ensembles. We 
will see how a description in terms of symmetric spaces emerges naturally from the cor- 
responding integration manifolds. We discuss this connection in detail for the gaussian 
and circular ensembles, which are the simplest and most frequently used random matrix 
theories. 

Historically, the gaussian ensembles were studied first, and only subsequently the circular 
ensembles, in which the requirement of statistical independence of the matrix elements is 
relaxed. We will, however, start by addressing the circular ensembles which turn out to 
be simpler to discuss from a symmetric space point of view. The circular ensembles are 
related to symmetric coset spaces of positive curvature, while the gaussian ensembles are 
related to euclidean type symmetric spaces of vanishing curvature, which are identified 
with a subspace of a Lie algebra. Therefore the identification of the symmetric space is 
slightly more tricky for the gaussian ensembles. As we will see, the identification with 
euclidean spaces is due to the translational invariance of the ensembles. In the following 
subsections we then discuss, in less detail but in the same spirit, a few further examples 
of random matrix ensembles which we have selected for their widespread applications in 
various branches of physics. 



8.3.1 Circular ensembles 

[1] Physical applications of the circular ensembles 

The circular ensembles are used in the description of physical systems which can be char- 
acterized by a unitary matrix S, typically a scattering matrix, each of whose elements 
give the transition probability from one state to another. For instance, these ensembles 
describe the statistics of scattering phase shifts (pi in the theory of mesoscopic systems (for 
a review and references see These systems are modeled by a phase-coherent disor- 

dered scattering region connected to ideal leads. Scattering of electrons against randomly 
distributed impurities takes place in the disordered region. 

The scattering matrix S relates the incoming and the outgoing wave amplitudes of the 
electrons. If we have propagating modes (channels) at the Fermi level, we can describe 
them by a vector of length 2N of incident modes /, /' and a similar vector of outgoing 
modes O, O' in each lead, where unprimed letters denote the modes in the left lead and 
primed letters the modes in the right lead. Then the scattering matrix is defined by 
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(8.13) 



Following a standard notation |2HJ EH] the scattering matrix has the following block struc- 
ture 



where r, r', t, t' are N x N reflection and transmission matrices. The unitarity of S, 
warranted by flux conservation (|/p + |J'p = |0p + |0'p), implies that the four matrices 
tt^, t't'^, 1 — rr^, 1 — rV'^ have the same set of real eigenvalues Ti, T/v called transmission 
eigenvalues (0 < Tj < 1). These determine the physical conductance of the scattering 
region and will be discussed further in connection with transfer matrix ensembles. 

[2] Definiton of the circular ensembles 

The circular ensembles were introduced by Dyson in the early sixties. They are the 
simplest possible matrix models, since the probability distribution is a constant. By prob- 
ability distribution we mean, in general, a distribution P{S) such that P{S)dS is the 
probability that a matrix from the ensemble takes a value between S and S + dS (for an 
introduction see ^T], Ch. 9). The non-trivial features of the model only appear in the 
probability distribution of the eigenvalues. 

Mimicking the gaussian case (which we will discuss in the next subsection) Dyson proposed 
three ensembles labelled by an index j3, which takes one of the values (3 =1, 2, and 4 and 
counts the number of degrees of freedom characterizing each matrix element. The index 
f3 is called the Dyson index of the matrix model. For the gaussian matrix models, the 
three cases correspond exactly to systems with time-reversal invariance and rotational 
invariance (/5 = 1), no time-reversal invariance {(3 = 2), and time-reversal invariance and 
no rotational invariance (/? = 4) jllj . 

• In the circular orthogonal ensemble labelled by /5 = 1, S" is a unitary symmetric 
N X N random matrix. 

• In the circular unitary ensemble {(3 = 2) S is a. unitary N x N random matrix. 

• Finally in the case of the circular symplectic ensemble {(3 = 4) S is a unitary self- 
dual random matrix with N x N quaternionic matrix elements. Generally speaking. 




(8.14) 
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a quaternion can be expressed in terms of 2 x 2 matrices as Q = aoO"o + ia ■ a, where 
the Oj are real or complex numbers and (uo, c?) = (1, (Xi, (J2, (T3) are the 2x2 unit and 
Pauli matrices. (For an introduction to quaternion algebra, see or jSHI, Ch. 8.) 
By self-dual we mean = S where {S^)ij = <J2{S'^)ija2 (see [TH ED] for further 
details). 

The terminology for the circular ensembles may seem odd. It stems from the gaussian 
ensembles, where the stability subgroups of invariance are orthogonal, unitary, and 
symplectic, respectively, for the values 1, 2 and 4 of the Dyson index. 

[3] Properties of the circular ensembles 

Below we recollect some general properties of the circular ensembles. As we will see, these 
properties have analogous realizations in all the ensembles that we will study. 

Invariance. It is easy to see that P{S)dS for the three circular ensembles is invariant 
under the similarity transformations 

S^USU-^ (/5 = 1,4) S-^VSW {(3 = 2) (8.15) 

where U is unitary {(3 = 1) or unitary with quaternion elements (/? = 4) and V, W are 
unitary. 

Diagonalization. In each case the matrix 5* can be diagonalized by a similarity transfor- 
mation like the ones in eq. ()8.15p . The eigenvalues lie on a unit circle and take the form 
e**^' (thereof the name "circular ensembles"). They are doubly degenerate in the symplectic 
case P = 4. 

Probability distribution of the eigenvalues. A simple calculation JT] shows that the 
Jacobian of the transformation from the space of unitary matrices to eigenvalue space is 
given by: 

J;3(0,)ocnie*^'-e*'^^f (8.16) 

i<j 

It gives rise to correlations between eigenvalues. Since the original distribution function 
P{S) of the random matrices is constant, the Jacobian coincides with the probability 
distribution function of the eigenvalues: 



88 



(8.17) 



[4] The symmetric spaces associated to the circuleir ensembles 

In this paragraph our main interest is identifying the integration manifold of the three 
circular ensembles. 

Orthogonal ensemble. Every symmetric unitary matrix S can be written as 



where C/ is a generic unitary matrix. However, this mapping is not one-to-one. If we 
assume that S — U'^U — V'^V, then it is easy to see that the matrix relating the two 
expressions, R = VU~^, is unitary and satisfies R^R = 1. Hence R must be real and 
orthogonal. Thus we see that the manifold of the unitary symmetric matrices is actually 
the coset U{N)/0{N), due to the above mentioned degeneracy. From the point of view 
of the physical properties of the ensemble nothing changes if we perform the restriction to 
an irreducible symmetric space. Then the manifold becomes SU{N)/SO{N). We will 
systematically perform the reduction to irreducible symmetric spaces in the following. 

Unitary ensemble. Finding the integration manifold is trivial in the j3 = 2 case, where we 
simply have unitary matrices without any further constraint and the manifold is simply 
the group U{N). Like for /3 = 1, if one is interested in irreducible symmetric spaces a C/(l) 
factor must be taken out, and the manifold becomes SU{N). 

Symplectic ensemble. This case strictly resembles the /? = 1 case discussed above. Any 
self-dual unitary quaternion matrix can be written as 



where U this time is a 2N x 2N unitary matrix. By following the same reasoning as in 
the case of the orthogonal ensemble, we see that the same matrix S can be obtained using 

^•^In the partition function 



extracting such a U{1) factor from the integration manifold just amounts to redefining Z by a constant. 




s = u''u 



(8.20) 




(8.19) 
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a new unitary matrix V obtained from the previous one by the transformation V = BU 
where B is constrained by B^B = BB^ = 1. By definition of the duahty operator ^ 
this means that B is a. symplectic matrix. Thus, extracting a. U{1) factor, the manifold 
coincides with SU{2N)/Sp{2N). 

By comparing with Table ^ in subsection 15.21 in Part I of the present review, we see that 
the integration manifolds of the three circular ensembles are exactly the first three coset 
spaces (described in the Cartan notation as A, AI and All) of positive curvature in the list 
of possible irreducible symmetric spaces. This is the main result of this paragraph. 



8.3.2 Gaussian ensembles 

[1] Physical applications of the gaussian ensembles 

Some of the physical applications of the gaussian ensembles were discussed already in the 
introduction to this section, in paragraph 18.11 In general gaussian ensembles describe 
hamiltonian ensembles. This is because they correspond to flat symmetric spaces, as we 
will see below. The corresponding positive curvature ensembles are the ensembles of the 
scattering matrices of the same systems, since in general the scattering matrix is expressed 
as 5" = e*^. As we have seen, applications include spectra resulting from complicated 
many-body interactions like those taking place in neutron resonances in atomic nuclei, the 
electronic energy levels inside tiny metal grains at low temperature (that depend on the 
shape of the grain), or generally the spectra of electrons moving in a random potential 
with no further symmetries present. 

The energy spectrum of classically chaotic systems is another example of apparently gaus- 
sian random behavior. A realization of such a system is a chaotic billiard. The motion 
of the billiard ball is determined by the shape of the billiard table (which is chosen to 
be irregular in some way) and it is drastically different for small variations in the initial 
trajectory. Corresponding quantum billards may also be considered, represented by a free 
quantum particle confined to a finite part of space. Its discrete energy spectrum will be 
determined by the Laplacian on the space, and may possess varying degrees of randomness. 
In all these systems the statistical properties of the energy spectrum can be described by 
a gaussian random matrix ensemble. 

[2] Definition of the gaussian ensembles 

Historically the first type of ensembles to be studied, the gaussian ensembles introduced 
by Wigner |17j and Dyson [TJ HO] are defined by the probability distribution 
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Pp{H) ~ e-^*'-^^^) (8.21) 

where H is a. hermitean N x N matrix and V{H) is a quadratic potential. We can define 
a partition function 

Z = J dHPp{H) (8.22) 

wfiere dH is an invariant Haar measure. It can be sliown ^I] tliat tlie form of PpiH) is 
automatically restricted to the form 

Pf3{H) = exp(-a tiH^ + btiH + c) (8.23) 

(a > 0) if one postulates statistical independence of the matrix elements Hij. Note that 
PfsiH) can be cast in the form 

Pp{H) ~ e-"*^^' (8.24) 
by simply completing the square in the exponent. 

Depending on the nature of H we can distinguish three cases labelled by /3 =1, 2, and 4: 

• In the orthogonal ensemble {/3 = 1) H is a.n N x N hermitean symmetric random 
matrix. It immediately follows that the probability distribution Pi{H) and the in- 
tegration measure dH are invariant under all real orthogonal transformations of H. 
This ensemble is used when the random matrix H models the time-reversal and 
rotation invariant Hamiltonian of a system with integral spin ([IHl, CI], Ch. 2). 

• In the unitary ensemble {f3 = 2) H is an N x N hermitean random matrix and 
P2{H)dH is invariant under all unitary transformations of H. A system without 
time-reversal invariance has this type of Hamiltonian lllj. 

• In case of the symplectic ensemble {f3 = 4) H is an N x N hermitean self-dual 
random matrix with quaternionic elements. It can be shown that the entries of such 
a matrix are all real quaternions which can be expressed in terms of 2 x 2 matrices 
as Q = aoCTo + id ■ a, where the are real numbers and (cro,a) = (1, cri, (T2, era) 
are the 2x2 unit and Pauli matrices. Pi{H)dH is invariant under all symplectic 
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transformations of H. The Hamiltonian of a system with time-reversal invariance, 
no rotational invariance and half-odd integral spin is of this type ITT] . 

As we can see, the symmetry classes are distinguished by their behavior under time reversal 
(TR) and spin rotation (SR). 

[3] Properties of the gaussian ensembles 

As we did for the circular ensembles, let us review some general properties of the gaussian 
ensembles. 

Statistical independence of the entries. With the choice of potential in eq. ()8.24|) 

the three ensembles have independently distributed elements, since the potential tiH^ = 

'I2ij \Hij I . 

Invariance. As implied above, Pp{H) and the integration measure dH are separately 
invariant under the transformation 

H UHU-\ (8.25) 

where U is an orthogonal, unitary or symplectic N x N matrix depending on the value of 
f3. However, it is important to notice that the symmetry group of P/siH) dH is larger and 
consists of rotations by the matrix U like in eq. ()8.25|1 . and addition by square hermitean 
matrices: 

H UHU-^ + H' (8.26) 

This will play an important role in the identification of these ensembles with a symmetric 
space. 

Diagonalization. For each matrix H there is a matrix U that maps it onto its eigenvalues: 

H = UAU-\ A = diag(Ai,...,Ajv) (8.27) 

where Aj {i = 1, A^) are real eigenvalues (if /3 = 4, they are twofold degenerate.) A simple 
calculation shows that the Jacobian of the transformation from the space of hermitean 
matrices H to eigenvalue space is given by 
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J/3({Ai}) oc J]^ |Ai - Xj 



(8.28) 



i<j 



where j3 =1, 2, and 4 in the orthogonal, unitary and symplectic case, respectively (see ^I] 
and 61J, Ch. 3). 

Probability distribution of the eigenvalues. Under the transformation ()8.27j) . the 

integration measure dH factorizes into an integral over the symmetry subgroup K and an 
integral over the eigenvalues with integration measure </s({Aj}) Hi d^i- The integral over 
the subgroup can be performed immediately and gives a constant equal to the volume of 
the subgroup. With this in mind and the random matrix potential given by V{H) = 
the joint probability density for the eigenvalues becomes 



The proportionality constant is not relevant if one normalizes the integral to unity (however, 
if needed, exact normalization constants for many integrals like the one in (|8.29|) were given 
by Hua in [HT]). 

The spectral correlations are due to the Jacobian. Eq. ()8.29|1 can be rewritten in the 
general form 



The logarithmic pair potential leads to repulsion between the eigenvalues. This is the 
so called spectral rigidity which is one of the most important features of random matrix 
theory. It is not present in uncorrelated Poisson distributions. The interpretation of ()8.30|) 
is that the probability of finding the i'th eigenvalue in the interval between Aj and Aj + dAj 
is proportional to -?/?({ Aj}) Hi dAj. 

[4] The symmetric spaces associated to the gaussian ensembles 

Let us now rephrase the above results in a group theoretical language. This will allow 
us to make contact with the analysis in section El Our goal here is to show that the 

""^"^This form of the potential gave rise to the name "gaussian ensembles". Of course, a gaussian random 
matrix potential can be used in any context, as it often is, because it is the simplest one to deal with. 



^/3({A.})ocniA.-A,|V^^.-i^ 



(8.29) 




(8.30) 
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integration manifolds of the above ensembles are symmetric spaces of zero curvature. This 
is essentially due to the translational invariance of eq. ()8.26p . 

Let us start with the (3 = 1 case. As a first step let us limit ourselves to traceless matrices 
only. This is the equivalent in the present context of the choice that we made in the previous 
subsection when we eliminated the U{1) factor from the coset. Also in the present case the 
reason is that we want to obtain a description in terms of irreducible symmetric spaces. 
The remaining degree of freedom acts trivially and is usually neglected. 

We now try to describe the set of all real, symmetric, and traceless matrices as a symmetric 
space X = G/K. By comparing with the example in subsection 13. II in Part I of the review, 
this can easily be done by identifying it with the algebra subspace SL(N, R)/SO(N). This 
simply amounts to taking a matrix from the set of generic, traceless, real N x N matrices, 
write it as a sum of a symmetric and an antisymmetric matrix, and then eliminate the 
antisymmetric part. The resulting space is exactly the set of all real, symmetric and 
traceless matrices that we were looking for. By taking only traceless matrices, we have 
eliminated an factor from the algebra GL(N, R), corresponding to the absolute value 
of the determinant, and we end up with the algebra SL(N,R), of which we then choose 
only the symmetric generators. 

Following the discussion of flat symmetric spaces in paragraph 15.11 the set of all real, 
symmetric and traceless matrices (let us call it P) can be realized as a zero-curvature 
symmetric space by the following two steps: 

(1) The group of affine transformations of P = SL(N, R)/SO(N) into itself is given 

by the semidirect product of the subgroup K = SO{N) and the invariant algebra 
subspace P. The action of G° on P is given by g{p) = kpk~^ + a, where g = 
{k,a) G G^, k & K, p, a G = G^/K = P. But this affine transformation is 
exactly of the same type as the symmetry transformation of eq. ()8.26|) that leaves 
the gaussian ensembles invariant (in fact, for the algebra subspace under discussion 
it is this transformation for (3 = 1). 

(2) The algebra subspace P now forms an abelian subalgebra (identical to the Cartan 
subalgebra) of the algebra G°, because it is additive: [P, P] = 0. Therefore G° = 
K © P has an abelian ideal, and is non-semisimple. The curvature tensor of these 
spaces is zero by equations (j5.1|) and (j5.2|) of subsection 15.11 Also, the commutation 
relations (j3.2j) defining a symmetric space are satisfied for the algebra elements in K 
and P. Hence X^ = G^/i^ = P is a fiat symmetric space with a euclidean geometry. 

By performing a similar analysis for f3 =2, 4 we obtain the following general result: The 
gaussian ensembles labelled by (3=1, 2, and 4 consist of hermitean square matrices belong- 
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ing to algebra subspaces SL(N, R)/SO(N), SL(N, C)/SU(N), and SU*(2N)/USp(2N), 
respectively. The algebra SL(N, C) is the algebra of hermitean, complex and traceless 
N X N matrices, and the subalgebra SU(N) defines the subgroup of invariance of the 
unitary ensemble. The algebras SU*(2N) and USp(2N) were formally defined in a foot- 
note in subsection 14.21 Again, they are exactly the algebra of hermitean, self-dual N x N 
quaternion matrices and the subalgebra of unitary symplectic matrices. From Table ^ we 
see that these three symmetric spaces correspond to algebra subspaces in Cartan classes A, 
AI and AIL The integration manifolds of the circular ensembles are the positive curvature 
coset spaces corresponding to the same Cartan classes. 

[5] Remarks concerning fiat symmetric spaces 

The fiat symmetric spaces may be seen as limiting cases of their curved counterparts. 
Below we make a few observations that illustrate this. 

As we discussed insubsection 15.11 the symmetric spaces are naturally organized in triplets 
of negative, zero and positive curvature. The simplest way to look at the zero curvature 
symmetric spaces is to see them as limiting cases of one of the other two symmetric spaces 
(with nonzero curvature) in the same triplet. This limiting procedure is related to group 
contraction. As the simplest example of this, we can take the group S0{3). As we saw 
in subsections 14.11 and 15.11 the algebras SO(3,R), SO(2,l;R) and E2 are all related to 
each other. The coset S0{3)/ S0{2) can be identified with the unit 2-sphere S"^, a positive 
curvature space (cf. subsection 12. 3|) . By performing the Weyl unitary trick on the algebra 
SO(3) we obtain the non-compact algebra SO(2, 1; R). As we showed in subsection 14. 2t 
the coset space S0{2, 1; R)/ 5*0(2) corresponds to the hyperboloid H^, a negative curvature 
symmetric space. The commutation relations for the non-semisimple euclidean group E2 
are a limiting case of the commutation relations for SO(3) and SO(2, 1) in the limit of 
infinite radius (cf. subsection 15.11) . The symmetric space E2/ S0{2) is thus a zero curvature 
symmetric space, isomorphic to the euclidean plane and a limiting case of the corresponding 
positive and negative curvature spaces in the same triplet. 

A simple and important example of a zero curvature symmetric space, which has several 
analogies with the case at hand, is the fiat euclidean space in four dimensions. It is well 
known that this space can be realized as the coset of the euclidean Poincare' group P 
with respect to S'0(4): P ~ P/S0{4). The translations of the Poincare' group play the 
role of P, they are isomorphic to euclidean space and have all the characteristics of a 
symmetric space of vanishing curvature. The above observation that the zero curvature 
spaces can be obtained as limits of positive curvature spaces can be exemplified as follows. 
We can realize the euclidean Poincare' group as a suitable limit of the 5*0(5) group. In 
this limit the coset 50(5)750(4) (which is nothing else than the four dimensional unit 
sphere) exactly becomes the euclidean four-dimensional space. 
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We can rephrase the above result in still another (more geometric) way. The zero curvature 
symmetric spaces can be seen as tangent spaces of their curved partners. This is clear both 
from the above example and from the algebraic structure of the space. 



8.3.3 Chiral ensembles 



A natural generalization of the previous results is the extension to rectangular random 
matrices. The study of this class of ensembles was initiated in the mid-1980 's by Cicuta 
et al. [HS] and further discussed in reference jHSj where they were named Laguerre en- 
sembles (due to the fact that the orthogonal polynomials associated to them are Laguerre 
polynomials). An updated review on this subject can be found for instance in |27j . 

[1] Physical applications of the chiral ensembles 

One of the most relevant applications of the chiral ensembles is in the study of the infrared 
limit of gauge theories. In particular, the universal properties of Dirac spectra that depend 
only on the global symmetry of the euclidean Dirac operator i]/), can be reproduced by 
substituting the integral over gauge fields in the euclidean partition function by an integral 
over an appropriate random matrix ensemble. The work by J. Verbaarschot et al. [HI] 
related to this subject dates back to the 1990's. It is in this framework (for reasons which 
will soon be clear) that these ensembles came to be named "chiral" . In this paragraph we 
discuss some of their characteristics with emphasis on their connection with the theory of 
symmetric spaces. 

Chiral ensembles are also relevant in connection with random flux and random hopping 
problems. In this subsection we will discuss two papers by Mudry, Brouwer, Simons and 
Altland ^49i^ iHfl] where chiral ensembles are realized. 

Chiral random matrix theory and QCD. The interest in the low-lying spectrum of 
the Dirac operator is due to the Banks-Casher formula 



im) = ^ (8.31) 

relating the magnitude of the quark condensate in QCD in a finite volume V to the spectral 
density of the Dirac operator at the origin. The part of the Dirac spectrum near the origin is 
therefore of interest in studying the mechanism of spontaneous breaking of chiral symmetry, 
since the quark condensate is an order parameter for the chiral transition. Since we are 
interested in the connection with symmetric spaces in this review, we will not further report 
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the results of these important apphcations. We only give a sketchy description below and 
refer to the literature cited for more detailed information. 

Due to the presence of chiral symmetry ({75, il/)} = 0) in four space-time dimensions, the 
Dirac operator has the block structure 

where W is rectangular, say p x q [p > q), thus u = p — q corrsponds to the number 
of fermionic zero modes. In QCD, these may describe the zero modes in the field of an 
instanton. Then u equals the winding number (topological charge). Zero modes may also 
originate in the bulk of the spectrum due to repulsion between eigenvalues. The total size 
N = p + q of the matrix in ()8.32|) corresponds to the finite space-time volume V in the 
gauge theory. In the calculations to be sketched below, the rectangular submatrices W in 
the euclidean Dirac operator i p will be replaced by random matrices T chosen from an 
appropriate ensemble. The symmetries of the chiral ensembles of random matrices (i.e. 
the ones having a block-structure like the one in eq. ()8.32|) ) are chosen on the basis of 
the fermion representation (fundamental, adjoint) and the number of quark colors. The 
various possibilities lead to the presence (or absence) of an anti-unitary symmetry 

[iP,Q]=0 (8.33) 

where Q is an anti-unitary operator. The presence or absence of the symmetry ()8.33|1 
and the explicit form of Q determines whether the Dirac operator can be represented as 
a matrix with real [jS = 1), complex [jS = 2) or quaternion real {f3 = 4) elements [HE]- 
The corresponding random matrices are in Cartan class BDI, AIII, and CII. In the gauge 
theory partition function in euclidean space 

Z = f DA e~^[^] n det(^ + rrif) (8.34) 

(where the integral is over the gauge field configurations, S[A] is a gauge field action, the 
product of fermion determinants is over the flavor degree of freedom and is the fermion 
mass for flavor /) the integral over gauge fields is then substituted with a gaussian average 
over a random matrix T. Since the Dirac operator has the block-structure in equation 
(|8.32|) . we can substitute W with T to get a much simpler model 
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with the same symmetries as the original partition function [DT is a Haar measure and 
the size of T is p x g). From this random matrix theory we obtain the eigenvalue density 
and correlators, in particular their universal microscopic limit. 

Both the random matrix theory and the field theory map onto the same low-energy effective 
partition function in the mesoscopic regime describing static Goldstone modes in a finite 
volume. This partition function is expressed as an integral over a coset manifold G'/K' 
related to spontaneous symmetry breaking and the emergence of Goldstone modes. Here 
G' is the original symmetry group of the Lagrangian and K' is the unbroken subgroup. It 
is worth noting that these Goldstone manifolds obtained in jEEl lEZl lEH] for /5 = 1, 2, and 
4 for QCD in both three (non-chiral ensembles) and four (chiral ensembles) space-time 
dimensions are saddle-point manifolds that appear in the large-A^ limit and should not 
be confused with the symmetric spaces G/K identified with the original ensembles. These 
Goldstone manifolds are the fermionic parts Mp of the symmetric supermanifolds Mb x Mp 
discussed by Zirnbauer in 

In addition to the finite volume partition function of static Goldstone modes expressing the 
quark mass dependence, one can obtain the flavor symmetry (or parity in odd space-time 
dimensions) breaking pattern. One can also derive sum rules constraining the eigenvalues 
of the Dirac operator in a finite volume jSHl 168j. 

Chiral ensembles in the context of random flux and random hopping problems. 

It is of interest in the theory of quantum transport to study how the conductance of a 
system behaves analytically as one or more of the dimensions of the system go to infinity. 
This approach is called the scaling theory of localization, and deals with the conditions 
under which a system is metallic, insulating, or undergoes a transition between the two 
states. Such a transition can be disorder-induced, i.e. it can depend on impurities in the 
metal. Another important factor is dimensionality. We postpone a further discussion of 
transport problems to the subsection on transfer matrix ensembles. In this paragraph we 
will only mention some results relating to chiral ensembles. 

In the random magnetic flux problem one studies the localization properties of a spinless 
electron moving in a plane perpendicular to a static magnetic field of random amplitude 
and vanishing mean. In two dimensions, it is difficult to reach a conclusion about the 
localization properties of the states based on numerical data. Therefore it is easier to study 
a quasi-one dimensional configuration. This was done by Mudry et al. in reference 
where an M x lattice (M >> N » 1) corresponding to a thick quantum wire with 
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weak disorder was studied. Although the study of such wires is a vast field, we will limit 
ourselves here to reporting only one of the results of this interesting paper, so as not to go 
beyond our scope, which is to give examples of physical manifestations of chiral random 
matrix symmetry classes. 

In the case at hand, the system is governed by a Hamiltonian Ti (whose detailed form 
the interested reader can look up in the original paper) through a Schrodinger equation 
Till) = ei/j. Away from the band center e = 0, the localization properties of the particle are 
those of the standard gaussian unitary ensemble. In other words, the Hamiltonian 7i of 
the system is in Cartan class A. Exactly at the band center e = 0, an additional symmetry 
of the transfer matrix, called chiral or particle-hole symmetry, changes the symmetry class 
of the Hamiltonian into the chiral unitary symmetry class AIII. (Note that the presence 
of the magnetic field breaks time-reversal symmetry, so the chiral orthogonal and the 
chiral symplectic class can not be realized in this system.) Thus the random magnetic fiux 
problem provides a physical realization of the chiral unitary ensemble. We will refer to the 
above results again in connection with transfer matrices, when we enumerate some of the 
physical manifestations of symmetric spaces of negative curvature. 

Related to the random fiux problem is the random hopping problem, in which a particle 
hops on a lattice with random hopping amplitudes. Also in this problem one finds that 
the point at the center of the band is special. A delocalization transiton in such a system 
in one dimension goes back to work by Dyson (see the discussion and references in jSHI)- A 
quasi-one dimensional hopping model with weak staggering was investigated in [SU]. The 
lattice consisted of coupled chains of length L, where L >> N. In the vicinity of the 
Fermi energy the lattice model was approximated by a continuum model. For /3 = 1,2,4 
respectively, the Hamiltonian has the symmetries of a chiral ensemble and belongs to 
Cartan classes BDI, AIII, and CII, respectively. For our analysis the important point here 
is that the Hamiltonian in each case is a realization of a chiral symmetry class, just like 
the Dirac operator in four-dimensional QCD and the Hamiltonian at the band center in 
the random fiux problem. For details and other results of the quasi-ld random hopping 
problem we refer to the original paper. 

[2] Some properties of the chiral ensembles 

Invariance and diagonalization. Just like in the circular and gaussian ensembles, the 
rectangular random matrix T above can be diagonalized through a transformation with 
unitary matrices U and V: 

T = UAV-^ (8.36) 
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The matrix A is a diagonal real matrix with positive elements corresponding to the non-zero 
eigenvalues of T. The partition function in ()8.35|) is invariant under this transformation. 

Probability distribution of the eigenvalues. It can be shown that the Jacobian for 
the transformation from matrix to eigenvalue space for the chiral ensembles is given by 

^,^({M)ocniA?-A,YnAf^'^"' (8-37) 

i<j k 

where {Aj} are the eigenvalues of the rectangular matrix T. In terms of the new variables 
Xi = Xf ()8.H7|1 takes the well-known form 

;3(^ + l)-2 

JpiM)^E\x^-^/Ex^ ' (8.38) 

i<j k 

Here u is the number of zero modes of the chiral matrices. 

[3] The symmetric spaces associated to the chiral ensembles 

The chiral block-structure of the random matrices already gives a hint that they must 
belong to the subspace P of some algebra (see subsection 13. 5|) . Depending on the number 
of degrees of freedom of the matrix elements (i.e. the value of f3: 1, 2, or 4), it is not hard 
to guess that this subspace is identified respectively with SO(p, q)/(SO(p) (S> SO(q)), 
SU(p,q)/(SU(p) ® SU(q)), or USp(p, q)/(USp(p) ® USp(q)) (in this case p, q have 
to be even^^). Obviously, these algebra subspaces are symmetric spaces of the euclidean 
(zero-curvature) type. They correspond to Cartan classes BDI, AIII and CII in Table ^ 
Above we have enumerated several physical realizations of these symmetry classes. 



8.3.4 Transfer matrix ensembles 

The transfer matrix ensembles appear in the theory of quantum transport, in the random 
matrix theory description of so called quantum wires. In these pages we will only discuss 
the part of the theory which is relevant for our purpose, the study of the mapping between 
random matrix theory and symmetric spaces. For further information on the experimental 
and theoretical issues we refer the reader to the excellent introductory review on random 
matrix theory and quantum transport by Beenakker ^B] . 

i^Our notation is such that USp{2p) = USp{2p, C) = U{p, Q). 
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[1] Physical context of the use of transfer matrix ensembles 



The natural theoretical framework for describing mesoscopic systems is the Landauer the- 
ory Within this approach Fisher and Lee proposed the following expression for the 
conductance in a two-probe geometry (a finite disordered section of wire to which current 
is supplied by two semi- infinite ordered leads): 



where t is the N xN transmission matrix of the conductor (see eq. ()8.14j) ). is the number 
of scattering channels at the Fermi level and Ti, T2 ■ ■ ■ T/v are the eigenvalues of the matrix 
ttl The Ti's are usually referred to as transmission eigenvalues. The constants e and h 
denote the electronic charge and Planck's constant, respectively. 

The transmission eigenvalues Tj are related to the scattering matrix and to the transfer 
matrix, which define equivalent descriptions of the impurity scattering process. In the next 
paragraph we will define the relationship between the physical degrees of freedom Tj and 
the transfer matrix. 

[2] Definition of the transfer matrix 

The scattering matrix of a mesoscopic conductor was defined in eq. ()8.13|) in the context 
of circular ensembles. Like in subsection 18.3.11 if we have N propagating modes at the 
Fermi level, we can describe them by a vector of length 2N of incoming modes /, /' and 
a similar vector of outgoing modes O, O' in each lead. Let again the unprimed letters 
denote the modes in the left lead and the primed letters the modes in the right lead. While 
the scattering matrix S relates the incoming wave amplitudes /, /' to the outgoing wave 
amplitudes O, O' (see eq. fl8.13|) ). the transfer matrix M relates the wave amplitudes in 
the left lead to those in the right lead: 



Following the notation used in subsection 18.3.11 and in [SHlEni, the scattering matrix has 
the block structure 



G = GoTTitt^) = GoY.Tr 



n 




(8.39) 




(8.40) 




(8.41) 
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where r, r', t, t' are N x N reflection and transmission matrices. As we saw in subsection 
18.3. H the unitarity of S imphes that the four matrices tt^ , t't'\ 1 — rr\ 1 — rV"!" have 
the same set of eigenvalues Ti,...,T/v called transmission eigenvalues (0 < Tj < 1). The 
eigenvalues Aj of the matrix 



Q = h^M^M + {M^M)-^ -2) (8.42) 



are related to the transmission eigenvalues by 



A^ = (8.43) 



The Aj are non-negative. In terms of these, M can be parametrized as 




M=\l •] X^\[":,\^UTV (8,44) 



where m, m', f , v' are unitary N x N matrices (related by complex conjugation: u' = u*, 
v' = v* ii M e Sp{2N, R) or M e SO*{AN), see below) and A = diag(Ai, Ajv). In case 
spin-rotation symmetry is broken, the number of degrees of freedom in ()8.44j) is doubled 
and the matrix elements are real quaternions. 

[3] The symmetric spaces associated to the transfer matrix ensembles 

Transfer matrices are strongly constrained by various physical requirements. As a result 
they belong to one of three different symmetric spaces. We discuss in some detail below 
how these symmetric spaces are obtained. 

The physical requirements of flux conservation, presence or absence of time-reversal sym- 
metry, and presence or absence of spin-rotation symmetry lead to conditions on the transfer 
matrix. These conditions determine the group G to which M belongs. Flux conservation 
leads to the following condition on M ^\ 



i.e. M preserves the (2A^ x 2A^) metric S^. This means that M belongs to the pseudo- 
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unitary group SU{N, N) (M has to be continuously connected to the unit matrix so we take 
the connected component of U{N, N)). If the Hamiltonian is invariant under time- reversal, 
the condition on the transfer matrix is [HHl 



M* = S^MS^ 



( 



1 

1 



) 



(8.46) 



It is easy to check that together with the condition of flux conservation, this implies 



But J is the skew-symmetric form invariant under the non-compact symplectic group 
Sp{2N, R) (see e.g. [21 E3 ED]), thus M belongs to this group in case time-reversal sym- 
metry is present. 

If the Hamiltonian contains a spin-orbit interaction (i.e. in case of magnetic impurities in 
the conductor), spin-rotation symmetry is broken. In this case the presence of an extra 
spin degree of freedom doubles the number of components in the vectors /, I', O, O' of 
incoming/outgoing wave amplitudes (the components become spinors). Formally, the same 
conditions ()8.45|1 . ()8.46j) on M are valid for flux conservation and time-reversal symmetry, 
but M is now a matrix of N x N real quaternion elements [H] . If time-reversal symmetry 
is broken, we get M e SU{2N, 2N) like before. If it is conserved, a condition analogous 
to ()8.46|) is valid, with the only difference that the matrices now act on an iV-dimensional 
vector space of real quaternions. In this case M belongs to the group SO* (AN). This 
is the connected component of the group of linear transformations that preserve a skew- 
symmetric bilinear form on a quaternionic vector space (see. e.g. [72], paragraph 7.2; cf. 
also Hiiffmann |2]). 

Invariance of the transfer matrix ensembles. As we can check using the parametriza- 
tion eq. (jH311), rotating M by a matrix W G U{N), SU{N) x SU{N) x f/(l), or U{2N) 
(if M G Sp{2N,R), SU{N,N), or SO* {AN), respectively), gives a new transfer matrix 
M' = WMW~^ = U'TV with the same physical degrees of freedom {Ai, Aat}, since 
the matrix F is unchanged. This means that in each case F belongs to a coset space 
G/K, where M & G and W & K. These three ensembles of transfer matrices, correspond- 
ing to different physical symmetries, are usually named (in analogy with the nomencla- 
ture of the standarad gaussian ensembles) the "orthogonal" {Sp{2N, R)/U{N)), "unitary" 



{SU{N,N)/SU{N)xSU{N)xU{l)) and "symplectic" {SO* (AN) /U {2 N)) ensemble. This 



nomenclature is a bit unfortunate, since the stability subgroups do not correspond to it. 
These coset spaces are irreducible and non-compact. Thus they are symmetric spaces of 





(8.47) 
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negative curvature, as is also evident from Table [TJ They correspond to Cartan classes CI, 
AIII, and CII, respectively. 

New transfer matrix ensembles in the context of random tight— binding models. 

In connection with the chiral ensembles we discussed the Hamiltonian at the band center 
in a random magnetic flux problem. It is of the chiral symmetry class AIII This 
means that the transfer matrix for the same physical system belongs to Cartan class A, 
i.e. it is a physical realization of the coset space SL{N,C)/ SU{N). We also discussed 
the random hopping problem [50^ in the same context and saw that the Hamiltonian for 
f3 = 1,2,4 has chiral symmetry. This means the transfer matrices of the systems belong 
to coset spaces SL{N, R)/SO{N) (Cartan class AI), SL{N,C)/SU{N) (Cartan class A), 
and SU*{2N)/USp{2N) (Cartan class All). These are new realizations of transfer matrix 
ensembles not found in the applications discussed above. In the same type of problems, 
also other transfer matrix ensembles are manifest. We refer to [HI ESI EDI ESI EEl EZ| 
for more details. Of course, to each hamiltonian ensemble corresponds a transfer matrix 
ensemble and a scattering ensemble, so in a certain sense, all the boxes in Table |21 should be 
filled with "physical" ensembles. As we have already noted, the correspondence between 
hamiltonian ensembles and the transfer matrix ensembles for the same systems was given 
in Table I of reference jHj. 

By inspection of the parametrization of eq. ()8.44|1 . we realize that the (generalized) eigen- 
values Aj are not exactly the radial coordinates in the three symmetric spaces. By a change 
of coordinates (see eq. ()8.51|) below), we will be able to "disentangle" the coordinates in 
such a way that this parametrization, and the DMPK equation to be discussed in the next 
paragraph, become much more tractable. 



8.3.5 The DMPK equation 

What is needed at this point to complete the transfer matrix description of the quantum 
wire, is the explicit expression for the probability distribution of the {Aj} as a function of 
an external parameter, the length L of the the quantum wire. The standard approach to 
this rather non-trivial problem has been to find some dynamical principle so as to obtain 
an "equation of motion" for the probability distribution of the {Aj} as a function of L, and 
then (hopefully) to solve the equation and obtain the probability distribution. 

The construction of the equation for the probability distribution was completed during 
the eighties, at least in the case of quasi-one-dimensional wires {L >> W where W 
is the thickness of the wire), by Dorokhov [74J, and independently by Mello, Pereyra, 
and Kumar 59^ (for /5 = 1) by studying the infinitesimal transfer matrix describing the 
addition of a thin slice to the wire. The resulting evolution equation for the eigenvalue 
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distribution P({Aj},s), where s is the dimensionless length of the wire, is usually known 
as the Dorokhov-Mello-Pereyra-Kumar (DMPK) equation. The only assumptions which 
are needed to obtain this equation are: 1) that the conductor is weakly disordered, so 
that the scattering in the thin slice can be treated by using perturbation theory, and 2) 
that the flux incident in one scattering channel is, on average, equally distributed among 
all outgoing channels. It is exactly this second assumption which restricts the DMPK 
equation to the quasi-lD regime, where the finite time scale for transverse diffusion can 
be neglected. The results of j^S] were subsequently generalized to /5 = 2,4 in refs. [7^ 171]. 
The DMPK equation can be written as 



dP 

^ = DP (8.48) 

where s is the wire length measured in units of the mean free path I: s = L/l, and D can 
be written in terms of the {Aj} as follows: 



D =-Y: t^A,(1 + X,)J{X)^J{Xr\ (8.49) 
7 ~i oXi dXi 

with 7 = (3N + 2 — (3. (3 is the symmetry index of the ensemble of scattering matrices, in 
analogy with the well-known Wigner-Dyson classification, and J{X) = ^({A„}) is given 
by 



JW = ]l\^j-M^ (8.50) 

i<j 

The solution of this equation will be discussed in the forthcoming subsections. For the 
moment, let us only point out that J(A), which is the Jacobian of the transformation from 
random matrices to eigenvalues and eigenangles in the gaussian ensembles, has nothing 
to do with the Jacobian between the space of transfer matrices and the space of radial 
coordinates in the present case. Its appearance in the DMPK equation has historical 
reasons, as the authors tried to mimic the well-known gaussian ensembles. This choice of 
coordinates makes the DMPK equation asymmetric and ultimately hard to solve. 

The relationship between the Jacobian and the radial part of the Laplace— 
Beltrami operator. The DMPK equation can be rewritten in a more tractable form by 
making the change of coordinates 
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\i = sinh^ Xi 



(8.51) 



This introduces the "right" Jacobian for the transformation to eigenvalue space in the 
transfer matrix ensemble. This Jacobian turns out to be the function ^^{{xi}) = ^"^{x) 
given by 



J{{xi}) =^^(x) = nisinh^a;^ - sinh.'^ Xi\'^Y[\smh2xk\ (8.52) 

i<j k 

Let us introduce a new operator B defined as 



B = Y.—e{x)j-r\x) , (8.53) 

dxk dxk 

It is easy to see, by direct substitution, that B is related to the DMPK operator by: 



D = ^B (8.54) 
27 



Thus we can rewrite the DMPK equation as 



— = —BP (8.55) 

By comparing eq. ()8.53p with eq. ()6.26p of section 16.21 we see that the operator B is 
related to the radial part of the Laplace-Beltrami operator A'^ on the negative curvature 
symmetric spaces associated to the transfer matrix ensembles by A'^ = J~^BJ. This is 
the first indication of a general, important relation which we will discuss in more detail in 
section 18.51 below. 

Similar Fokker-Planck equations can be derived for all the symmetry classes. In [HH! the 
equation was derived in conjunction with a tight-binding model for the cases in which the 
transfer matrix is in the standard ensembles and the Hamiltonian in the chiral ensembles. 
In [76 it was derived for systems with Hamiltonians of the Bogoliubov-de Gennes (BdG) 
type (see the next subsection). The latter apply to quasiparticle transport at the Fermi 
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level in a disordered superconducting wire. In this case charge is not conserved, but one 
can study transport of heat and spin. 

In standard systems, localization occurs when the wave function of an electron undergoes a 
transition from a Bloch wave extending throughout the sample to a localized exponentially 
decaying form ?/'(r) ~ e"*"/^, where ^ is the localization length. As a result, in the localized 
regime the conductivity is dominated by the lowest eigenvalue and decreases exponentially 
with the length of the wire. In the context of disordered wires, the seven universality 
classes pertaining to chiral and BdG Hamiltonians are more interesting than the standard 
ones in that they can display a departure from exponential localization 

8.3.6 BdG and p— wave ensembles 

[1] Physical applications of the BdG ensembles 

The so called NS ensembles are examples of Bogoliubov-de Gennes, or BdG, ensembles. 
They get their name from the physical structures they describe The symmetry classes of 
these new ensembles are realized in normal metal-superconductor (NS) heterostructures. 
These are mesoscopic systems composed by a normal conductor in conjunction with a 
superconductor. Because of Andreev reflection at the NS interface, this system is different 
from the conventional one. Just like the additional chiral symmetry of the Dirac operator 
and the symmetries of the transfer matrix of a normal mesoscopic conductor give rise to new 
ensembles, so the symmetries of the so called Bogoliubov-de Gennes (BdG) Hamiltonian 
give rise to four new symmetry classes depending on whether time-reversal (TR) and/or 
spin-rotation (SR) symmetry is present. The BdG Hamiltonian is a first-quantized version 
of the BCS Hamiltonian in a mean-field approximation. It has an additional particle-hole 
grading which is absent in the Hamiltonian for a normal metal. This gives rise to a discrete 
particle-hole symmetry. 

[2] The symmetric spaces associated to the BdG ensembles 

In jS] the authors show that the BdG Hamiltonian Ti. (actually z7i) belongs to one of four 
symmetry classes depending on which symmetries are present. Since the presentation in 
jH] is excellent, we will just summarize the results here. The space to which iTi belongs is 
an algebra or an algebra subspace and is either SO(4N) (no TR and no SR), USp(2N) 
(only SR), SO(4N)/U(2N) (only TR), or USp(2N)/U(N) (both TR and SR). All these 
algebra subspaces can be considered to be symmetric spaces of zero curvature, by the 
construction we have repeatedly discussed. By comparing with Table ^ we see that they 
can be identified with Cartan classes D, C, Dili and CI respectively. The Jacobian for the 
transformation by the stability subgroups to radial coordinates is [3] 



107 



JrA{Q^})^I[\l'-l'\'I[\l'^\' (8-56) 

i<j k 

where the pair (r, s) takes the values (2, 0), (2, 2), (4, 1), and (1, 1) for SO(4N), USp(2N), 
SO(4N)/U(2N), and USp(2N)/U(N), respectively. For a review of NS junctions and 
the BdG equation, see jSEj- 

Scattering matrices of the BdG type. The scattering matrix for an NS-type het- 
erostructure is obtained by exponentiation of the Hamiltonian: S = e*^. Since ili. is in the 
algebras or tangent spaces listed above for the respective symmetry classes, the scattering 
matrix is in the corresponding symmetric spaces of positive curvature. Further information 
on these ensembles can be found in |2f)j . 

[3] The symmetric spaces associated to the p— wave ensembles 

Recently, Ivanov [3| found realizations of Cartan classes B and Dlll-odd in the algebra 
SO(2N + 1) and the algebra subspace SO(4N + 2)/U(2N + 1) in p-wave superconduc- 
tors. The corresponding ensembles are called p-wave ensembles and are characterized by 
the presence of a zero mode. We recall that u = \p — q\ zero modes are also present for 
the chiral ensembles for p ^ q (for notation see Table P). These two classes were obtained 
for disordered vortices in p-wave superconductors with or without time-reversal symme- 
try, using the same Bogoliubov-de Gennes Hamiltonian as in [3 . As in the case of BdG 
ensembles, the identification of the p-wave ensembles with algebra (sub) spaces maps them 
onto flat symmetric spaces. For more comments on these ensembles see 0. 

[4] Vicious walks and BdG ensembles 

Recently a new and completely independent realization of the gaussian BdG ensembles has 
been proposed in the context of the two-matrix model description of the vicious walk in 
presence of a wall (see and references therein). In this case the peculiar symmetries 
of the BdG ensemble are a direct consequence of the boundary conditions imposed by the 
presence of the wall in the model. 



8.3.7 S— matrix ensembles 

[1] Physical applications of the S— matrix ensembles 

The important class of matrix models referred to as S-matrix ensembles was introduced a 
few years ago in [7HI EH] to describe the behavior of ballistic chaotic quantum dots. These 
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are microstructures in which impurity scattering can be neglected, and only boundary 
scattering is considered. 



[2] Definiton of the S— matrix ensembles 

One assumes that the scattering matrix, given in eq. ()8.41|) . belongs to the circular ensem- 
ble. Like in subsection 18.3.41 the transmission eigenvalues T„ are the eigenvalues of the 
submatrix combination tt\ and these are not simply related to the scattering phase shifts. 
(Recall that in the circular ensembles, the probability distribution ()8.17j) in the original 
formulation by Dyson is expressed in terms of scattering phase shifts.) By expressing Pjj 
in terms of the variables {Aj} of eq. ()8.43|) . 



one obtains a representation suitable for studying the transport properties through a quan- 
tum dot. 

The resulting distribution P^({Aj}) takes the form of a Gibbs disribution. A direct calcu- 
lation (TS] shows that 



[3] The symmetric space associated to the S— matrix ensemble for /3 = 2 

It appears that the S-matrix ensemble for (3 = 2 corresponds to the compact symmetric 



space SU{2N)/SU{N) x SU{N) x U{1) (for p = g in SU{p + q)/SU{p) x SU{q) x f/(l). 



which corresponds to the multiplicity of the short roots = 2(p — g) = 0, i.e. a root 
lattice of type BCn that reduces to C„ for p = q). Interestingly, as we will see in subsection 
18.4. H the two other ensembles (for (3 = 1, 4) do not find correspondence in the Cartan 
classification This issue will be discussed further in section ITUl 



(8.57) 




(8.58) 
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8.4 Identification of the random matrix eigenvalues and univer- 
sality indices 

Just like the integration manifolds of random matrix theories have been identified with 
symmetric spaces in the previous subsection, we will now identify the random matrix 
eigenvalues with the spherical radial coordinates, the Jacobian from matrix to eigenvalue 
space with the Jacobian of the transformation to spherical coordinates, and the random 
matrix symmetry indices with the multiplicities of the restricted roots of the pertinent 
symmetric spaces. 

Radial coordinates and the Weyl group. In subsection IH.HI of this paper we discussed 
the important role played by the reference frame of spherical coordinates in the study of 
symmetric spaces. In particular the radial spherical coordinates turned out to be given by 
the Cartan subalgebra of the subspace P. We saw that P = could be naturally divided 
into conjugation classes by the action of the stability subgroup K of the symmetric space. 
More precisely, we showed that any element p G P is conjugated with some element of the 
Cartan subalgebra by means of the adjoint representation of K: p = khk~^ where k & K, 
h = and H is in the Cartan subalgebra. The radial coordinates of p coincide with the 
set of eigenvalues of the matrix H. 

In the case of a gaussian ensemble, from eq. ()8.24|1 we see that the probability density 
depends only on these radial coordinates. This is not a coincidence. It happens in all the 
random matrix models that we will study, and can be considered one of their character- 
istic features. This does not mean that random matrix models in which also the angular 
coordinates play a role are not important. They appear in several interesting contexts, but 
they are in general much more difficult to study and require methods which are outside 
the scope of this review. 

We observe that the Weyl group is the symmetry group of the root system of the symmetric 
space and of the radial coordinates at the same time. The Weyl group can easily be found. 
It is simply given by all the possible permutations of the eigenvalues. The Weyl chamber 
(defined in section 12. 6|) is fixed once we choose a particular ordering for the eigenvalues. 
We may for example order them according to increasing magnitude. 

Jacobians and the root lattice. In the previous subsection we have discussed a number 
of different random matrix ensembles. For each of these ensembles we identified the integra- 
tion manifold with a suitable symmetric space of positive, negative or zero curvature. For 
each ensemble we now identify the Jacobian of the transformation to eigenvalue-eigenvector 
space with the Jacobian ()6.42|) f paragraph 16. 3|) of the transformation to spherical coordi- 
nates on the symmetric space identified with the integration manifold. In random matrix 
theory this Jacobian gives rise to the familiar geometric interactions between eigenvalues. 
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and determines the form of the joint probabihty distribution function of the latter. As 
we will see, the detailed form of the Jacobian fits exactly into the common and elegant 
framework provided by the classification of symmetric spaces. The curvature of these 
spaces, together with the restricted long, ordinary and short roots and their respective 
multiplicities, exactly determine its form. 

A few important remarks. Before discussing the various ensembles, let us make two 
observations which will hopefully help in clarifying the pattern which is emerging. In 
general, symmetric spaces of positive curvature correspond to ensembles of circular type, 
those of zero curvature to ensembles of gaussian type (i.e. hamiltonian ensembles), and 
those of negative curvature to transfer matrix ensembles. Algebraically, the curvature of 
the spaces is reflected in the explicit form of the Jacobian of the transformation to spherical 
coordinates. This can be understood by comparing with eq. ()6.42|) which we reproduce 
here for comparison: 

J^~Kq) = naefl+(a~^sinh(ag"))'"- (8.59) 
^^+^(g) =naei?+(a-^sin(ag°))-'^ 

Note that the root multiplicities are the same for the spaces of positive and negative 
curvature in the triplet corresponding to the same symmetric subgroup. Even though 
we have not defined a restricted root system for the zero-curvature spaces arising from 
non-semisimple groups, the same multiplicities characterize also the Jacobian pertaining 
to the zero-curvature member of the triplet. This should be understood as explained in 
the remark following equation ()6.30|) in subsection 16.21 

8.4.1 Discussion of the Jacobians of various types of matrix ensembles 

Let us now come back to the classification scheme. The root multiplicities are listed in 
Table ^ to which we systematically refer in the following. This means that a few integers 
and the sign of the curvature (+, — , 0) are enough to completely characterize a matrix 
ensemble. 

1) The gaussian ensembles 

As we saw in subsection 18.3.21 the standard gaussian ensembles are labelled by the Dyson 
index (3. It is noteworthy that this index is exactly equal to the multiplicity of ordi- 



111 



nary roots of the restricted root lattice characterizing a triplet of symmetric spaces (see 
comments in the previous subsection concerning the "root multiplicities" characterizing a 
non-semisimple algebra G°). The Jacobian in the gaussian case is given by eq. ()8.28|1 : 

J^'\{^^}) - Ul^^' - (8-60) 

i<j 

This case corresponds to the An-i root lattices with only ordinary roots. The Jacobian in 
()8.60|) has exactly the form given in eq. ()8.59|) for zero curvature spaces with a root lattice 
of type An-i {mi = nis = 0). The positive roots in this case are a = — ej {i > j) and the 
expressions entering in the Jacobian are the radial coordinates Aq, = A ■ a = Aj — Aj (see 
paragraph 16. 2j) . The absolute value in eq. ()8.60p corresponds to a choice of Weyl chamber 
(cf. paragraph 12. 6p Ai > A2 > ... > Xn- 

2) The circular ensembles 

The Jacobian of the circular ensembles of unitary scattering matrices was given in ()8.17|1 : 

>//3({0.})«nie*-e^'^^f (8-61) 
It is not hard to see that this Jacobian can be rewritten in the form 

mo 

(8.62) 

which is exactly the form given by eq. ()8.59|) for the Jacobian pertaining to the transfor- 
mation to spherical coordinates on a symmetric space of type An-i with only ordinary 
roots and of positive curvature, if we choose the radius a = 1/2. This allows us to identify 
these ensembles with the positive curvature symmetric spaces of Cartan classes A, AI and 
All, in agreement with what we had already concluded from the corresponding integration 
manifolds. 

3) The chiral ensembles 

The Jacobian given for the chiral ensembles in eq. ()8.37j] 



n 

i<i 



2 sin 
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WM)«niA^-A?i'niA.r^^^-^ (8.63) 

i<j k 

corresponds in the most general case to a root lattice with all types of roots. The chiral 
ensembles are algebra subspaces, and like for the gaussian ensembles we associate to it the 
restricted root lattice of the curved symmetric spaces in the same triplet. The restricted 
root lattices for the chiral ensembles are of type Bg or Dg for (3 = 1 and BCq or Cg for 
/9 = 2, 4. In case the root lattice is not of type A„_i, not only is Cj — ej a positive root but 
also Cj + Cj. Thus the positive roots are as follows: for Bg {cj, Cj ± Cj}, for Dg {e, ± Cj}, 
for BCg {cj, Cj ± Cj, 2ei}, for Cg {cj ± Cj, 2ei} {i ^ j always). Using the root multiphcities 
rrio = P, mi = P — 1, rris = (3\p — g| = /3z/ we see that the Jacobian is of type (j8.59|) . It can 
be rewritten 

■^^"HlAiD-m^- -A^Yni^^r, l3 = mo, a = ms + mi (8.64) 

i<j k 

From the Jacobian it is evident that in addition to the usual repulsion between different 
eigenvalues, Aj also repels its mirror image — Aj, and the eigenvalues are no longer trans- 
lationally invariant. This kind of ensembles are therefore called boundary random matrix 
theories. The boundary random matrix theories include chiral and NS ensembles. 

4) The transfer matrix ensembles 

The "orthogonal" and "symplectic" transfer matrix ensembles correspond to the same 
classes of symmetric spaces as the BdG ensembles (see below), but with negative curvature. 
The "unitary" transfer matrix ensemble also corresponds to a symmetric space of negative 
curvature, but it is of the chiral, not BdG type (see Table Hj). The relevant root systems are 
all of type C„ with long and ordinary roots. The multiplicities of these can be read from 
Table lU mi is always equal to 1 and mo = (3. Therefore the Jacobian given in eq. ()8.52|) . 

J^({A,}) oc n Isinh^Ai - si\i\i^Xjf\{smh{2\k) (8.65) 

i<j k 

is of the general form J^^\{\i]) if we note that the positive roots are {cj ± ej,2ei} and 
use the following identity for hyperbolic functions: 

sinh(Aj — \j) sinh(Aj + \j) = sinh^Aj — sinh^Aj (8.66) 
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5) The BdG ensembles 



In exactly the same way as above, using Tableland the trigonometric identity similar to 
equation ()8.(j(ij) 

sin(A.j — Xj) sin(Aj + Xj) = sin^Aj — sin^Aj (8.67) 
we check that the Jacobian for the BdG ensembles of scattering matrices 

^r,«({Ai}) oc n \sm^Xi - sin^Ajf nsin'(2Afc) (8.68) 

i<j k 

is the one corresponding to positive curvature spaces {J^^\{Xi])) and determined by the 
corresponding restricted root systems of type D„ (only ordinary roots) or C„ (long and 
ordinary roots) according to equation ()8.59|) . This Jacobian corresponds to scattering 
matrices of the BdG type, 5* = e*'^. For the corresponding algebra subspaces to which ili, 
belongs, the Jacobian was given in eq. ()8.56|) . Obviously, its structure is the one given 
by the same restricted root lattices according to eq. (|8.59|) . but for zero curvature spaces 

6) The S— matrix ensembles 

In terms of the transmission eigenvalues T„ the probability distribution of equation ()8.58|) 
(which is the same as the Jacobian for circular ensembles) can be rewritten as: 

P{{Ti})dTi = J{{Ti})dT,^ ll\Ti-T/ll\n\^dn (8.69) 

i<j k 

Changing variables and setting Tj = sin^ 9i we find 

J{{ei})de, ~ Y[\sm'^ei-smHj\'^Y[\smekf~'^\sm2ek\dek (8.70) 

i<j k 

As we have by now understood, this is the typical Jacobian for a positive curvature space, 
described by a root lattices of the BCn type with root multiplicities rrio = P, nig = 
/3 — 2, mi = 1, where the Dyson index as usual can take the three values (3 = 1, 2, 4. 
However, in this case we see a rather unexpected feature: it seems that only one of the three 
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ensembles can be mapped onto a symmetric space. It is the unitary ensemble {(3 = 2), which 
is described by a C„ root lattice, since for /? = 2, = 0. The other two cases correspond 
to choices of the root multiplicities which lie outside the classification of Table H We 
will come back to this problem in subsection where we will see that the ensembles 
characterized by Dyson index (3 = 1 and (3 = A can be understood in terms of a non-Cartan 
parametrization of standard symmetric spaces. 



8.5 Fokker— Planck equation and the Coulomb gas analogy 

An important tool in the theory of random matrices is the so called "Coulomb gas analogy" 
due to Dyson (see fTP, Ch. 8 for a review). The content of this analogy is the following. 
The probability distribution 



of random matrix eigenvalues is identical to the probability density of the positions {Aj} 
of unit charges in one dimension in a stationary potential that is given by 



if the random matrix potential is of the gaussian type. The temperature of this system is 
kT = (3~^. The force exerted by the potential W is of the Coulomb type. If in addition to 
this force one considers the Brownian motion of the unit charges in the presence of a time- 
dependent random fluctuating force and a frictional force, one can derive a Fokker-Planck 
type equation describing the evolution of the Coulomb gas with time t 

This analogy allows us to write the probability distribution of eigenvalues in the gaussian 
and circular ensembles as the asymptotic limit of a Brownian motion process. Similar 
equations can be written also for all the other ensembles of positive and zero curvature 
listed in the previous section. The relevant feature of these equations for the purpose of 
the present review is that the differential operator which describes the Brownian motion in 
the case of the positive curvature ensembles can be exactly related to the radial part of the 
Laplace-Beltrami operator on the corresponding symmetric space. Remarkably enough, 
if one performs a similar mapping for the negative curvature symmetric spaces (i.e. the 
transfer matrix ensembles) one exactly finds the DMPK equation. Below we will discuss 
all this in more detail, but let us first summarize the results: 





m{M) = ;.E^'-EHA.-A, 




i=l i<j 
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The DMPK equation is the exact analog for the negative curvature symmetric spaces 
of the Coulomb gas equation introduced by Dyson for the positive and zero curvature 
ensembles. 



• All these equations describe the free diffusion - this is the ultimate meaning of 
the Laplace-Beltrami operator - of the eigenvalues on the corresponding symmetric 



• For the positive curvature ensembles, in the large time asymptotic limit the eigen- 
values reach the distribution dictated by the Jacobian. In fact, if we had the whole 
Laplace-Beltrami operator (not just its radial part) we would expect a uniform dis- 
tribution on the symmetric space, which is exactly the probability distribution of the 
random matrices in the circular ensembles. Since the Dyson operator is the radial 
projection of the Laplace-Beltrami operator, the final result is that the eigenvalues 
are distributed according to the Jacobian of the transformation to radial coordinates. 

• For the negative curvature ensembles, in the asymptotic limit the eigenvalues depart 
from each other exponentially. This behavior has a deep physical meaning. Indeed, it 
is well known that for quasi-one-dimensional wires, as the length of the wire (which 
corresponds to the time in the Brownian motion analogy) increases, the wire eventu- 
ally reaches the insulating regime. This perfectly agrees with the Brownian motion 
picture. As the time increases the eigenvalues depart arbitrarily far from each other. 
As a result the conductivity is dominated by the lowest eigenvalue and decreases ex- 
ponentially, in agreement with the fact that we have reached the insulating regime. 

Let us now discuss these results in detail. 



8.5.1 The Coulomb gas analogy 

Let us concentrate on the gaussian case for definiteness. As we mentioned above, the prob- 
ability density Pj3{{Xi}) of eq. ()8.29p and ()8.7H) can also be interpreted as the probability 
density for the positions of unit charges constrained to move on a line and interacting 
through the potential ()8.72j) . The integral 



for the probability can thus be interpreted as a standard thermodynamic partition function 
with f3 playing the role of inverse temperature. A nice feature of this representation is that 



spaces. 
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we may obtain the distribution Pfs{{Xi}) as the asymptotic hmit of a Brownian motion 
type process. To this end the following steps are needed: 

First we must add a fictitious time dependence in Pf^. Then we must assume that the 
out-of-equilibrium dynamics of the Coulomb gas is of the Brownian motion type (free 
diffusion). As we implied above, this involves introducing, in addition to the Coulomb 
force, a time-dependent rapidly fiuctuating force giving rise to the Brownian motion, and 
a frictional force. As a consequence the time dependence of Pf3{{Xi},t) must be described 
by the following Fokker-Planck equation jTTj 

dP ^ ( 1 d^P d 1 
f%-^\~-i-m,)P]] (8.74) 

where / is the friction coefficient setting the time scale for the diffusion process, and 



^(A,).-™^-A,+i:^ (8.75) 

is the electric Coulomb-type force experienced by the unit charge at Aj. Standard manip- 
ulations show that the Fokker-Planck equation can be equivalently rewritten as 



BP ^ B ~ 3 ~ 

where J exactly coincides with the joint probability density for the eigenvalues of eq. (|8.29|) 



J({A,}) = Pp{{\}) = Cn |A. - A,| V^Sr^i^'/' (8.77) 
with C an undetermined normalization constant. We have 



limP^({A,},t)=P^({Ai}) (8.78) 

t — >CXD 

One can write down similar Fokker-Planck equations for the all the zero and positive 
curvature spaces. All these Fokker-Planck equations describe the approach to equilibrium 
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of the eigenvalue distribution. As mentioned above, the equation for the negative curvature 
spaces coincides with the DMPK equation. 



8.5.2 Connection with the Laplace— Beltrami operator 

Comparing the Fokker-Planck equation ()8.7(j|l in the case of positive curvature ensembles 
(in this case the potential is absent and JUK}) = J^^\{\i}), that is, the joint probability 
distribution of the eigenvalues is identified with the Jacobian of the transformation to 
radial coordinates) with equation ()6.26|) we see that the Fokker-Planck operator 



^ d ~ 8 ~ 



is related to the radial part of the Laplace-Beltrami operator. 



A'^ = V ^JW-^ (8.80) 

where g° are the radial coordinates identified with the eigenvalues Xj. More precisely, the 
two operators are related by 



A'^ = J-^FJ (8.81) 

Similarly, in case of negative curvature spaces the operator B of equations ()8.531 18.54|) . 
that is proportional to the DMPK operator 



^ = E 7^^'(^)7^r'(x) = 2^D (7 ^ /?iV + 2 - (3) (8.82) 
is mapped to A'^ by the following relation 

A'^ = r\x) B eix) (8.83) 



118 



after identifying ^^(x) with the Jacobian J({xj}) = JIkj I sinh^ Xj — sinh^ Hi I sinh 2xj| 
of the transformation to radial coordinates (eq. ()8.52|) ) and the radial coordinates them- 
selves with the eigenvalues Xi related to Aj by ()8.51|) : Aj = sinh^Xj . 

Equation (jS.SHj) . that relates the DMPK operator to the free diffusion on a symmetric 
space, is one of the main results of the mapping between random matrix ensembles and 
symmetric spaces. As we shall see it also represents the starting point for the exact solution 
of the DMPK equation for the transfer matrix ensembles. 



8.5.3 Random matrix theory description of parametric correlations 

The original motivation behind the Coulomb gas approach to random matrix theory was to 
obtain the equilibrium probability distribution of the eigenvalues in a different and more 
physical way. In this respect the only interesting regime of the corresponding Fokker- 
Planck equations was the asymptotic t —>■ oo limit. However, it was later realized that this 
equation contains interesting physical information for all values of t. In particular it could 
be used for describing the adiabatic response to an external perturbation of the energy 
spectrum of a mesoscopic system ^U|. In this approach the role of the fictitious time is 
played by the perturbation parameter. 

An interesting application appears if the role of the external perturbation is played by an 
external magnetic field. In a random matrix description of a disordered conductor, there is 
a smooth transition between ensembles characterized by various values of /3 as a function 
of the magnetic field B around 5 ~ (for a review see j26j). When such a transition is 
completed, the level distribution becomes independent of the magnetic field. The random 
fluctuation of individual energy levels is still present, however, as a function of B. Such 
random fluctuations are described by Fokker-Planck type equations similar to the DMPK 
equation, but containing the Laplace-Beltrami operator on symmetric spaces of positive 
curvature 



8.6 A dictionary between random matrix ensembles and sym- 
metric spaces 

Probably the simplest way to summarize what we have said in this section is to write down 
the equivalences between random matrix and symmetric space concepts in a table. This is 
what we have done in Table El 
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Table 2: The correspondence between random matrix ensembles and symmetric spaces 



Random Matrix Theories (RMT) 


Symmetric Spaces (SS) 


circular or scattering ensembles 


positive curvature spaces 


gaussian or hamiltonian ensembles 


zero curvature spaces 


transfer matrix ensembles 


negative curvature spaces 


random matrix eigenvalues 


radial coordinates 


probability distribution of eigenvalues 


Jacobian of transformation to radial coordinates 


Fokker-Planck equation 


radial Laplace-Beltrami equation 


Coulomb gas analogy 


Brownian motion on the symmetric space 


ensemble indices 


root multiplicities 


Dyson index (3 


multiplicity of ordinary roots {(3 = rrio) 


boundary index a = (3{v + 1) — 1 


multiplicity of short and long roots {a = rris + mi) 


translationally invariant ensembles 


SS with root lattice of type An 


boundary matrix ensembles 


SS with root lattices of type i?„, C„, Dn or BCn 


pair interaction between eigenvalues 


ordinary roots 



9 On the use of symmetric spaces in random matrix 
theory 

In this section we discuss some of the applications of the mapping between random matrix 
ensembles and irreducible symmetric spaces outlined in section |H1 The main application, 
which is a natural consequence of the Cartan classification of symmetric spaces, is a ten- 
tative classification of the random matrix ensembles. We will discuss this important issue 
in subsection I9.H while in subsection 19.21 we discuss how the symmetries of the spaces are 
reflected in the random matrix ensembles. A second natural application is related to the 
orthogonal polynomial approach to the construction of eigenvalue correlation functions of 
random matrix ensembles JT]. As it turns out, the orthogonal polynomials associated to 
random matrix ensembles whose integration manifolds are symmetric spaces, are all of the 
classical type and can be directly constructed from the knowledge of the curvature and 
root multiplicities of the underlying symmetric space. We will discuss this issue in para- 
graph 19.31 Finally in paragraph 19.41 we discuss the applications of symmetric spaces to the 
transfer matrix ensembles which appear in the description of quantum transport. This is 
probably the field in which the knowledge of the mathematical structure of the underlying 
symmetric spaces is most helpful. 
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9.1 Towards a classification of random matrix ensembles 



The most interesting consequence of the identification of random matrix ensembles with 
irreducible symmetric spaces is that the Cartan classification of the latter naturally induces 
a classification of the matrix ensembles. It is important to observe, however, that we 
should look at this Cartan classification more as a useful framework for organizing the 
various matrix ensembles than as a rigid classification of all possible ensembles. It is not 
hard to construct ensembles that fall outside of the Cartan grid. This can be done by 
suitably constraining the random matrices. There is nothing wrong with the ensembles 
that go beyond the Cartan classification. However, a consequence of such constructions is 
that for these ensembles, the beautiful and powerful properties of symmetric spaces are no 
longer applicable, and it becomes much more difficult to extract meaningful information 
or predictions for the random matrix theory^^. 

At the same time it appears that most of the matrix ensembles which have physically 
interesting realizations belong to the Cartan framework. It is possible that this is merely 
a consequence of the fact that, in trying to describe a physical problem, one prefers to 
use ensembles that are simple to deal with, even though this may imply stronger approx- 
imations. Let's mention an example which illustrates the issue at hand. Recently various 
generalizations of the DMPK transfer matrix ensembles discussed in this review have been 
proposed 81 as an attempt to avoid the "quasi-one-dimensional" approximation involved 
in the standard DMPK equation. The resulting equations cannot be mapped to a sym- 
metric space, and thus cannot be solved exactly or even asymptotically. In spite of this, 
important information on the expected behavior of the eigenvalues can all the same be 
obtained by means of suitable perturbative expansions ^Tj . 

In Table 0] we have included the random matrix ensembles discussed in section |H1 in the 
Cartan classification of irreducible symmetric spaces. This classification was presented in 
Table ^ The scattering matrix for an NS-type heterostructure is obtained by exponen- 
tiation of the Hamiltonian S = e*^. Since ili. is in the algebras or tangent spaces for 
the respective symmetry classes, the scattering matrix is in the corresponding symmetric 
spaces of positive curvature. The scattering matrix ensembles of NS systems have been 
listed in Table El with the notation -8^0,™,,^^ of the same type as the circular and 

S-matrix ensembles. We list the ensembles using this kind of notation, which is more 
consistent than some of the traditional names given to ensembles in the past with abuse 
of language (for example, an ensemble was called "unitary" if /3 = 2, even though the 
stability subgroup was not unitary). 

The notation used here is as follows. An ensemble is labelled by a letter indicating the 
type of ensemble and alluding to the traditional name. Let C stand for circular ensembles, 

^^See the discussion in section HO. II below for a remarkable exception to this fact. 
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G for gaussian ensembles, P for the p-wave ensembles of Ivanov (these are of the BdG 
type, but have a zero mode), B for Bogoliubov-de Gennes ensembles, T for transfer matrix 
ensembles, S for S-matrix ensembles, and x for chiral ensembles), a superscript taking 
the values +, 0, — indicating the curvature of the space, and three subscripts mo,mi,ms. 
Since the pair {(3, a} can be the same for some pairs of distinct ensembles, it is better to 
keep all three root multiplicities in the label. For example, B^q o indicates the BdG (NS) 
ensemble corresponding to a symmetric space of positive curvature with root multiplicities 
{rrio, mi, rris} = {2, 0, 0}. In this way each ensemble is uniquely labelled. 

Let us stress that the empty spaces of Table El do not mean that a corresponding random 
matrix ensemble doesn't exist. Following the previous discussion, it is easy to see that 
for each symmetric space (of arbitrary curvature) one can construct a perfectly consistent 
matrix ensemble. The empty boxes simply mean that such an ensemble still has not found 
a relevant physical application or realization (or, possibly, that we are not aware of such 
a realization). It is likely that with time all the empty boxes in Table El will be occupied 
with physically interesting applications. Note also that the type of restricted root system 
changes within the same Cartan class for the ensembles labelled by two integers {p,q}, 
depending on whether p > q or p = q. The S-matrix and transfer matrix ensembles have 
p = q {u = 0) and have been written on the corresponding line, while the chiral ensembles 
may have p > q or p = q. 

It is important to note that the symmetric space associated to a random matrix theory 
can be given either for an ensemble of random Hamiltonians 7i, for an ensemble of random 
transfer matrices Ai, or for an ensemble of random scattering matrices S. Thus, the 
symmetric space associated to the transfer matrix group of a given system (for example a 
quantum wire) is different from the symmetric space associated to the Hamiltonian Ti or 
the scattering matrix 5* of the same physical system. A table of correspondences between 
the Ai and Ti. descriptions was given in [6|. 



9.2 Symmetries of random matrix ensembles 

Some known symmetries of the random matrix ensembles can be understood in terms of 
the symmetries of the associated restricted root lattice. In particular, ensembles of An 
type are characterized by translational invariance of the eigenvalues. This translational 
symmetry is seen to originate in the root lattice: all the restricted roots of the An lattice 
are of the form (cj — ej). The Wigner-Dyson (circular and gaussian) ensembles are of the 
translation invariant type. 

For all the other types of restricted root lattices (-B„, C„, D„ and BCn) this invariance is 
broken and substituted by a new Z2 symmetry giving rise to the reflection symmetry of 
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Table 3: Irreducible symmetric spaces and some of their random matrix theory realizations. 
The random matrix ensembles with known physical applications are listed in the columns 
labelled and X~ and correspond to symmetric spaces of positive, zero and negative 

curvature, respectively. Extending the notation used in the applications of chiral random 
matrices in QCD, where v is the winding number, we set z/ = p — q. The notation is C 
for circular, G for gaussian, x for chiral, B for Bogoliubov-dc Gcnncs, P for p-wave, T for 
transfer matrix and S for S-matrix ensembles. The upper indices indicate the curvature, 
while the lower indices correspond to the multiplicities of the restricted roots characterizing 
the spaces with non-zero curvature. To the euclidean type spaces ~ /K, where 
the non-semisimple group is the semidirect product X (g) P, we associate the root 
multiplicities of the algebra G = K © P. 



Restricted 
root space 


Gartan 
class 


G/K {G) 


G*/K (G^'/G) 


mo 


m 




X+ 




X- 


An-1 
An-1 

An-1 


A 

AI 

All 


SU{N) 

SU{N) 
SO{N) 
SU{2N) 
USp{2N) 


SU{N) 
SL{N,R) 

SOiN) 
SU*l2N) 
USp{2N) 


2 

1 

4 












r+ 


^2,0,0 

pO 

1,0,0 

pO 

^^4,0,0 


^2,0,0 

Ti,o,o 

T4,0,0 


BGq {p>q) 


AIII 


SU{p+q) 


SU{p,q) 


2 


1 


2u 




X2,1,2j/ 




Cq ip=q) 


SU{p)xSU{q)xU{l) 


SU{p)xSU{q)xU{l) 


+ 

^2,1,0 


^2,1,0 


Bn 


B 


S0{2N + 1) 


S0{2N+1,C) 
SO(2N+l) 


2 





2 




pO 

2,0,2 




Gn 
Gn 


C 
CI 


USp{2N) 

VSp{2N) 


Sp{2N,C) 
USp(2N) 
Sp{2N,R) 


2 
1 


2 
1 






1^2^2,0 


bo 

-^2,2,0 

BO 

^1,1,0 


^2,2,0 

Ti,i,o 


SU{N)xU{l) 


SU{N)xU{l) 


BGq (p>q) 


CII 


USp(2p+2q) 


USp{2p.2q) 


1 


3 








-V 1.3.1,/ 




Cq (p=<f) 


/-.SV(2/,)x(-,SV(2r/) 


rSp{2p)yA-Sp(2q) 




1,3,0 


Dn 
Gn 

BC^ 


D 

Dili 

Dili 


SO{2N) 

so{m) 


SO(2N,C) 
SO{2N) 
SO* (AN) 


2 

4 
1 



1 
1 




1 


R+ 

^2,0,0 

R+ 


B° 

^2,0,0 

^4,1,0 

pO 

^ 1,1.1 


T2,0,0 
'^4,1,0 


SU{2N)xU{l) 

SO ( IN +2 ) 


SU{2N)xU{l) 

SO* (IN +2) 


SlJ(2N + l)xU(l) 


SlJ(2N + l)xlJ(l) 






Bq ip>q) 


BDI 


SOip+q) 


SO(p,q) 


1 





V 




Xi,o,i/ 




Dq ip=q) 


SO(j>}xSO{q) 


SO(j))xSO{q) 




Ti,o,o 
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the eigenvalues that we discussed in the context of the Jacobians. Since these ensembles 
are characterized by the presence of a boundary (not always immediately evident) with 
respect to which they are reflection invariant, they are called boundary random matrix 
theories (BRMT in the following). They include all the remaining ensembles that are not 
of the circular or gaussian type. 



9.3 Orthogonal polynomials 

An important role in the study of matrix ensembles is played by the set of polynomials or- 
thogonal with respect to the random matrix theory integration measure. These polynomials 
come into the picture when rewriting the Jacobian for the transformation to eigenvalue 
space in terms of a product of Vandermonde determinants. By adding linear combinations 
of the rows, the determinant can be written as a determinant of monic polynomials (a 
polynomial Pn{x) is called monic if Pn{x) = x" + (9(x"~^)), for example 



N 

TT(xi-Xo)~ det x*""^ = det Pj_i(xo) (9.1) 

J--^/ l<i,i<7V J l<i,j<N * ^ ^ 

If these polynomials are then chosen orthogonal with respect to the measure^^ w{x)dx = 
g-Ary(x)^^^ where V{x) is the random matrix potential 

y^e-^^(^)p„(x)P„(x)tix = hj^n (9.2) 

(here I is some interval on the real axis and /i„ is a normalization factor), the eigenvalue 
correlation functions can be expressed in terms of a Christoffel-Darboux kernel An 
arbitrary /c-point correlation function is defined as 



p{xi,...,Xk) = n dxj P{xi,...,xn) (9.3) 

j=fc+i 

where P{xi, ...,xn) is the joint eigenvalue distribution of the random matrix model. The 
quantity p{xi, ...,Xk)dxi...dxk equals the probability of finding one eigenvalue in each of 

^^The factor N in the exponent is common in QCD related applications. If the weight function is odd, 
so called pseudo-orthogonal polynomials may be used. We will not discuss them here, since our goal is 
just to remind the reader of the general mechanisms that make orthogonal polynomials useful. 
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the intervals between Xj and Xj + dxj {j = 1, ...,k). 

The general formula for the fc-point function turns out to be 

p{xi,...,Xk) = —^^ ^det^^KN{xi,Xj) (9.4) 

In particular, the spectral density is simply 

pix) = KN{x,x) (9.5) 

In eqs. fl9.4|l and ()9.5j) . Kiy( the Christoffel-Darboux kernel defined in terms of 

orthogonal polynomials 

K^ix,,x,) = iV-^ef (^(^?)+^(^?» Y: K'P,ix,)P,ix^) (9.6) 

A;=0 



It turns out that for all the matrix ensembles related to symmetric spaces, the associated 
orthogonal polynomials belong to the set of so called classical orthogonal polynomials. 

With the term classical one usually denotes three families of orthogonal polynomials: the 
Jacobi, Laguerre and Hermite polynomials, whose unifying feature is the so called Ro- 
driguez formula^^ which allows to construct the polynomials once the weight function p{x) 
and the function X{x) (which specifies the domain of support of the polynomials) are 
given: 



Here An is a normalization constant which can be obtained explicitly, but is irrelevant for 
our purposes. 

The domain-specifying function X{x) is a polynomial in x of degree < 2. The three 
possibilities are: 

^^Actually this formula is a generalization due to Tricomi of the original formula obtained by Rodriguez 
for the case of Legendre polynomials. 
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• if X{x) is simply a constant, then the polynomials are defined on the whole real line, 
the weight function must be p{x) = , and we find the Hermite polynomials. 

• if X{x) is a polynomial of first degree, one can always shift the origin so as to 
obtain X{x) = x. In this case the orthogonal polynomials are defined on the positive 
real axis, the weight function must be p{x) = x^e~^, and we find the Laguerre 
polynomials. 

• if X{x) is a polynomial of second degree, then one can always normalize it so as 
to obtain X{x) = {1 — x'^). In this case the orthogonal polynomials are defined on 
the interval [—1, 1], the weight function must be p{x) = (1 — + xY, and we 
find the Jacobi polynomials. The well known Gegenbauer, Chebyshev and Legendre 
polynomials are only special cases of Jacobi polynomials. 

Each one of these polynomial families is in one-to-one correspondence with a particular 
random matrix ensemble of table El The mapping is complete, i.e. all the ensembles of 
table El are covered. 

In particular the Hermite polynomials are related to the gaussian Wigner-Dyson ensem- 
bles (i.e., those corresponding to symmetric spaces defined by an An root lattice). Simple 
changes of variables allow one to show that the Laguerre polynomials are related to the 
gaussian BRMT's (these are the chiral, BdG andp-wave ensembles of zero curvature spaces) 
and the Jacobi polynomials to the circular BRMT's (i.e., ensembles related to symmetric 
spaces of positive curvature not of the type, like for instance the BdG scattering en- 
sembles). A thorough discussion of the correlation functions for these ensembles using the 
orthogonal polynomials listed above can be found for instance in jSEj (see also |M]). 

These results are well-known (for instance, in the early papers on chiral matrix ensembles, 
these ensembles were named Laguerre ensembles) and do not require any particular refer- 
ence to the symmetric space description of random matrix theories. What is interesting 
in our framework is that the parameters which define the polynomials can be explicitly 
related to the multiplicities of short and long roots of the underlying symmetric space and 
thus, by the identification in Tabled with the boundary universality indices of the BRMT. 
The relation is the following: 

Laguerre polynomials: 




— A _x 



X = 



rus + mi - 1 




2 
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Jacobi polynomials: 




(-1)" (1 - x)"" rf" / (1 + xY+P \ 
2"n! {l + x)P dx"" \{l - x)-"^-" J 



(-1 < x < 1) 



TUg -\- mi — 1 



nil - I 



(9.9) 



2 



a = 



2 



We see that A and p have the same expression in terms of and m;. Thus the BRMT's 
corresponding to Laguerre and Jacobi polynomials with the same X = p indices belong to 
the same triplet in the classification of Table El They are respectively the zero curvature 
(Laguerre) and positive curvature (Jacobi) elements of the triplet. This explains the so 
called "weak universality" which was observed a few years ago for the boundary critical 
indices of these ensembles i.e. the fact that the form (near the boundary) of scaled k- 
level correlators is the same for Laguerre and Jacobi ensembles if the symmetry parameter 
P is the same (see jSH] for a discussion of this point). The weak universality turns out to 
be simply a consequence of the organization in triplets of the symmetric spaces! 



9.4 Use of symmetric spaces in quantum transport 

One of the most interesting applications of symmetric spaces is in the transformation of the 
Fokker-Planck equations in random matrix theory into Schrodinger equations in imaginary 
time (where, as before, the time coordinate in transfer matrix ensembles is identified with 
the dimensionless length s of a quantum wire). As a consequence of this transformation, 
in the case (3 = 2 the degrees of freedom in the Schrodinger equation decouple and it can 
be solved exactly. 

This result traces back to the original work by Dyson jH3] and was later extended to 
boundary random matrix theories by various authors. We mentioned this important result 
already at the end of section [7| In paragraph 19.4. II we will study this mapping in detail in 
a case which is particularly relevant from a physical point of view, namely for the DMPK 
equation of transfer matrix ensembles. 

Another important consequence of the mapping is that in the "interacting cases" (for (3 = 1 
and 4) one can use the results discussed in subsections 16.31 and 16.41 on the zonal spherical 
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functions to obtain important information on the asymptotic behavior of the solutions to 
the DMPK equation. We will discuss this issue in paragraph 19.4.21 In subsection 19.4.31 we 
will see an example of a change of symmetry class (hence also of the underlying symmetric 
space description) induced in a quantum wire by switching on an external magnetic field. 
Finally in subsection l9.4.4l we shall discuss a general scheme (based on the DMPK equation) 
for constructing the scaling equations for the density of states of a quantum wire that covers 
all the Cartan symmetry classes. 

9.4.1 Exact solvability of the DMPK equation in the 13 = 2 case 

The exact solution of the DMPK equation in the (3 = 2 case was first obtained in a 
remarkable paper by Beenakker and Rejaei. Here we review their derivation in a 
slightly different language, trying to stress the symmetric space origin of their result. 

The starting point is the mapping discussed in paragraph 18.3.51 which we briefly recall 
here. By setting A„ = sinh^x„ (cf. eq. fl8.5H) ) the DMPK equation can be rewritten as 
(see equations (jH3^ and dHSSl)) 

BP 1 

-q; = i^mf ^'b P (9.10) 

where 

i{{xi\) = I sinh^a;^ — sinh^a;j|2 JJ^ | sinh2a;fc|2 (9-11) 

i<i k 

and A'^ is the radial part of the Laplace-Beltrami operator on the underlying symmetric 
space. 

At this point one can follow two equivalent ways. In the flrst one (which was the one 
followed in |13]) one makes use of the results discussed in section [71 in particular eq. ()7.5p . 
to map the DMPK equation into a Schrodinger equation. By comparing with eq. ()7.5|) we 
see that this simply requires the substitution 

P{{Xn}, S) = ^{{Xn})^{{Xn}, s) (9.12) 

A straightforward calculation shows that the DMPK equation then takes the form of a 
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Schrodinger equation in imaginary time: 



-^ = in-U)^ (9.13) 



where f/ is a constant and 7i is a Hamiltonian of the form 



^ 2j^[dxr ^ 27 ^(cosh(2a;0-cosh(2x,))2 ^ ^ 

At this point the main goal has already been reached: it is easy to see that if /3 = 2 the 
equation decouples and an exact solution can be obtained Before going into the details 
of this solution, let's remark that the above equation can be recast in a slightly different 
form, thus completing the chain of identifications DMPK equation — radial part of the 
Laplace-Beltrami operator — Calogero-Sutherland model. By using simple identities for 
hyperbolic functions, this Hamiltonian becomes 



/ 1 ^2 gf \ ( gl gl \ 



where gf = —1/2, g"^ = (3{(3 — 2)/4 and c is an irrelevant constant. This Hamiltonian, 
apart from an overall factor I/7 and the constant c, exactly coincides with the Calogero- 
Sutherland Hamiltonian ()7.4|) discussed in section [71 corresponding to a root lattice R = 
{±2xi, ±Xi ± Xj,i 7^ j} of type Cn with root multiplicities rrio = P, rrii = 1. The values of 
the coupling constants go, gi are exactly the root values given in eq. (|7.8|) 



2 ma{ma — 2)\a 



2 

(9.16) 



of section [71 for which the transformation from Ti into A'^ is possible. 

Let us now come back to the exact solution of the DMPK equation following Beenakker 
and Rejaei. 

As we have seen, for /5 = 2 7i is reduced to a sum of single-particle Hamiltonians TYq, 
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At this point, to solve the DMPK equation one simply has to construct the Green's function 
Go of the single-particle Hamiltonian Ho- This requires solving the eigenvalue equation 



Hoi^ix) = eij{x), (9.18) 

A standard analysis shows that the spectrum of Tio is continuous, with positive eigenvalues 
e = jk? /N and that the eigenfunctions ijjk{x) are real functions given by 



Tik tanh I — I sinh 2x 



1/2 



Pi(i,_i)(cosh2a;). (9.19) 



where Pu{z) denotes the Legendre functions of the first kind. From this we obtain the 
spectral representation of the single-particle Green's function Go 



Go{x,s\y) = {2n)-^ J ^ dk expi-—j iJk{x)iJk{y) 



1 I k'^s\ ( 7Tk\ 

- (sinh 2x sinh 2yY''^ j dk exp f — j k tanh f — j 

X Pwi,_i)(cosh2x)Pi(i,_i)(cosh2y) (9.20) 



The A^-particle Green's function G is related to the single-particle Green's function Gq 
through a Slater determinant, and the probability distribution of eigenvalues is related to 
G through a similarity transformation by the antisymmetrized eigenstate ^1^0(3;) = ^{x) 
{j3 = 2) of the iV-fermion Hamiltonian (for details see |13]). From the expression ()9.20p . 
imposing so called ballistic initial conditions (which essentially amount to requiring that all 
the eigenvalues are concentrated at the origin for s = 0), one finally obtains the probability 
distribution P({x„},s) for the eigenvalues 



Pi{xn},s) = G(s) ]^(sinh^ Xj — sinh^ Xj) ]^(sinh 2a;fc) 

i<j k 
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X Det 



dk exp 




]^)(cosh 2Xr 



(9.21) 



This is the exact solution of the DMPK equation for (3 = 2^^. 

The second approach rehes more heavily on the underlying symmetric space structure. The 
starting point is again the identification made in eq. (j9.1(J|) between the DMPK operator 
and the radial part of the Laplace-Beltrami operator on the underlying symmetric space. 
As a consequence of this identification, if $fc(x) (x = {xi, ■ ■ ■ , x^}, k = {ki, ■ ■ ■ , kjy}) is 
an eigenfunction of A'^ with eigenvalue k"^, then ^(x)^$fc(x) will be an eigenfunction of the 
DMPK operator with eigenvalue k'^/{2'y). The eigenfunctions of the A'^ operator (which 
are known in the literature as zonal spherical functions) have been widely discussed in 
subsections 16.31 and 16.41 As we have seen, by means of the zonal spherical functions one 
can define the analog of the Fourier transform on symmetric spaces: 



fix) = I mMx)j^, (9.22) 

(here we have neglected an irrelevant multiplicative constant). In particular, for the three 
symmetric spaces which are of interest for us one finds: 



\c{k)\' = \mm 



r(i + 4) 



(9.23) 



with 



\m\' 



n 

m<j 



r (i^^^) r (i 



r(f + z^)r(f + 2^ 



(9.24) 



where F denotes the Euler gamma function. This is a completely general result which we 
will use also in the next section. The problem is that in general the explicit form of the 



^^We remark that in [HHI, using the same technique as in 45j, the Fokker-Planck equation for the 
probabihty distribution of eigenvalues in systems with a chiral Hamiltonian was solved exactly in the case 
(3 = 2, and in [THI, the equation corresponding to a system with BdG Hamiltonian was solved in the 
presence of time-reversal symmetry (for two of the four BdG symmetry classes) . 
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zonal spherical functions involved in eq. ()9.22|) is not known. A remarkable exception to 
this situation is represented exactly by the (3 = 2 case for which we have P ES] (see eq. 
flA.4l in the appendix): 



det 



Qi 



Yli<:j[{kf — A;|)(sinh^ Xi — sinh^Xj)] 



(9.25) 



where the matrix elements of Q are: 



= F Q(l + ikm), ^(1 - ^^m), 1; -sinh^a;; 



(9.26) 



and F(a, b, c; z) is the hypergeometric function. 



Equations 



allow us to write the s-evolution of P{{xn}, s) from given initial 



conditions (described by the function fo{k)) as follows: 



dk 



\c{kW 



(9.27) 



By inserting the explicit expression of |c(/i;)p and by using the identity: 



r(.|) 



k Ilk 
— tanh — 
2 2 



(9.28) 



we end up with the following general expression for P{{xn}, s) with ballistic initial condi- 
tions (which, due to the normalization of $fc(x), simply amount to choosing fo{k) equal to 
a constant): 



Pi{xn},s) = [ax)]' dke 



27' 



$fc(x) 



,'n'ki 



|^^^^|,rifc,tanh(Y) 



(9.29) 



Inserting the explicit expression for ^k{x) from equations ()9.25|9.26p into ()9.29|) and using 
the identity 
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P,(z) = F(-z/,z/ + l,l;(l-^)/2) 



(9.30) 



we exactly obtain, as expected, the solution 1)9.211) found by Beenakker and Rejaei. This 
is a remarkable and non trivial consistency check of the correctness of this solution. 

9.4.2 Asymptotic solutions in the /? = 1,4 cases 

The power of the description in terms of symmetric spaces becomes evident in the /3 = 1 
and (3 = 4 cases, in which the interaction between the eigenvalues does not vanish and the 
first approach discussed in the previous subsection does not apply. On the contrary, the 
description in terms of zonal spherical functions (i.e., eq. ()9.27)l ) also holds in these two 
cases. Even though for /3 7^ 2 one does not know the explcit form of the zonal spherical 
functions, one can use the powerful asymptotic expansion ()fi.fi3|l discussed at the end of 
subsection 16.41 to get asymptotic solutions. In our context this expansion reads: 



where rk is the vector obtained acting with r & W on k {W denotes the Weyl group of the 
symmetric space). The important feature of eq. 1)9.31)) is that it is valid for all values of 
k, thus it can be used both in the metallic {k 3> 1) and in the insulating {k -C 1) regimes. 
This leads to expressions for the probability distribution of the eigenvalues and for the 
conductance, which can then be compared to other theoretical results using numerical 
simulations or experiments (see jlO] for a detailed discussion). 



9.4.3 Magnetic dependence of the conductance 

In [SH], the theory of symmetric spaces was applied to give a possible explanation for the 
discrepancy between the random matrix theory and non-linear sigma model analysis of 
the magnetoconductance in the weakly insulating, localized regime close to the Anderson 
transition of a disordered wire. The Anderson transition is a disorder-induced transition 
from the conducting to the insulating regime. For a review of the sigma model approach 
see Hg. 

Magnetoconductance is the change in the conductance of the wire due to the presence of a 
magnetic field (for a review see |26j). More precisely, it is the suppression of weak localiza- 
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tion (a quantum effect due to time-reversal symmetry) which appears when a magnetic field 
destroys time-reversal invariance. The disagreement between the random matrix theory 
approach and the non-linear sigma model approach was evident in the prediction for the 
magnetoconductance in the presence of strong spin-orbit scattering: the sigma model ap- 
proach gave a zero magnetoconductance, while the random matrix approach gave a strong 
negative value due to suppression of the localization length. A negative magnetoconduc- 
tance has been observed in experiments. However, as we will review below, also within 
the random matrix approach a zero magnetoconductance can be expected, depending on 
whether Kramers degeneracy of the eigenvalues is conserved or not. 

The conductance in the insulating regime is related to the localization length ^ by G = 
Goexp(— 2L/^), where L is the length of the sample. In the standard random matrix 
approach, ^ is proportional to [3. Thus the transition (due to the switching on of a magnetic 
field) between the ensemble characterized by /3 = 4 and the ensemble characterized by 
(3 = 2 means a negative contribution to G. 

The key observation made in [5^] was that the matrix ensemble to which the transfer matrix 
belongs corresponds, for /? = 4, to a symmetric space SO*{4N)/U{2N) with a root lattice 
of type Cn- This root lattice is characterized by two types of roots, long and ordinary, and 
therefore by two indices (3 = nio = 4 and rj = mi = 1. In the ensemble labelled by /3 = 4 
the eigenvalues are twofold degenerate. Therefore, if there are N degenerate scattering 
channels for (3 = 4 (so called Kramers degeneracy), for (3 = 2 there are 2N channels. 

Using the mapping of the DMPK equation onto a Calogero-Sutherland model (see sec- 
tion [7j), a generalized DMPK equation was derived to take into account the new in- 
dex 7]. As a consequence, it was shown that the localization length was not affected if 
one performs the simultaneous change /3 = 4 (3 = 2 and — > 2N, while keep- 

ing rj fixed. This corresponds to the transition between ensembles SO*{4:N)/U{2N) —>■ 
SU{2N, 2N)/{SU{2N) x SU{2N) x f/(l)) (the latter ensemble is characterized by /? = 2, 
r] = 1) and gives a zero magnetoconductance. Effectively, the level statistics depends on (3 
and A^, but the localization length depends more generally on (3, rj, and A^. 

Assuming instead that Kramers degeneracy is conserved {N N) we have a transition 
S0*{AN)/U{2N) USp{2N,2N)/{USp{2N) x USp{2N)) (the latter ensemble is char- 
acterized by /5 = 4, 7] = 3) with the result that ^ is smaller in the presence of a magnetic 
field. In this case we have a negative magnetoconductance. Thus we see that both possi- 
bilities, zero and negative magnetoconductance, are feasible also within the framework of 
random matrix theory. 
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9.4.4 Density of states in disordered quantum wires. 



Another interesting application of tlie DMPK equation is tlie construction of the scaling 
equation for the density of states of a disordered quantum wire which was recently discussed 
by Titov, Brouwer, Furusaki and Mudry in jHlEZl- Computing the density of states p(e) in 
the thermodynamic limit — »• oo in a quantum wire of infinite length at energy e turns out 
to be intimately related to the solution of the DMPK equation in presence of absorption 
ini EI! • Absorption is described by the addition of a spatially uniform imaginary potential 
iuj, uj > to the Hamiltonian. 

A DMPK-like equation can be obtained also in this case. It is very similar to the standard 
one (see eq.s ()8.53p and ()8.55p ). The only change is an additional term proportional to the 
absorption constant. 



ds 



— B -\ > 7- — smh 2xj 

27 Vp ~i ox^ 



Pi. 



(9.32) 



j=l ^-^3 

where / is the mean free path and vp the Fermi velocity. 



It is not difficult to see that this equation has the following stationary (i.e. L independent) 
solution: 

|J| ^ 



Zia] 



(9.33) 



where a = '-flu/vp is the adimensional absorption constant and Z{a) is a normalization 
constant which is needed to ensure that Pj^ is normalized to one. The key (and non- 
trivial) point is that, once Z{a) is known it is possible to show (by making an analytical 
continuation from iuj to e) that the density of states p(e) is obtained from Z{a) as follows: 



Pie) 



d 



-Re 



lim 

-i-yle/vp Oa 



a-^\ogZ{a) 



(9.34) 



(see IB] for a detailed derivation). It is interesting to observe that this procedure works for 
all the Cartan classes. 



In the non-standard symmetry classes of the Hamiltonian, the density of states is singular 
at the band center (for chiral Hamiltonians) or Fermi energy (for BdG Hamiltonians) . This 
point corresponds to extra symmetries of the Hamiltonian, an issue that was discussed 
in the context of chiral ensembles. The precise form of this singularity depends on the 
symmetry class and, for chiral Hamiltonians, on the parity of A^ (the number of channels). 
Apparently, there is a connection between the anomalous behavior of the density of states 
and the divergence of the localization length (signaling criticality) at the singular point 

m 
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10 Beyond symmetric spaces 



In the previous section we have discussed some of the possible apphcations of the theory 
of symmetric spaces in random matrix models. In this last section we discuss three issues 
which in one way or another go beyond the theory of symmetric spaces developed up to 
now. We will see that the power of the group theoretical methods developed in the first part 
of this review allows us to obtain interesting results also in some cases in which symmetric 
spaces appear not to be useful. At the same time the topics to be discussed represent new 
open directions of research which we hope will lead to future interesting results. 

In subsection 110.11 we will discuss the non-Cartan parametrization of symmetric spaces. 
We will see how it is possible to map exotic random matrix ensembles to non-standard 
(in a sense that will be clear) symmetric spaces. In subsection 110.21 we will discuss how a 
wide set of non-isotropic solutions to the DMPK equation can be constructed by simply 
resorting to the exact integrability of the associated Calogero-Sutherland models. This 
new tool might be useful in overcoming the quasi-lD constraint of the DMPK description 
of quantum wires. Finally, in subsection 110. 31 we discuss the triplicity of (in a certain sense) 
the most general potential in a Calogero-Sutherland model, the Weierstrass potential. We 
will see how this potential in various limits reproduces the components of an arbitrary 
triplet of symmetric spaces. 



10.1 Non— Cartan parametrization of symmetric spaces and S— 
matrix ensembles 

As we have seen in paragraphs 18.3.71 and 18.4. H apparently only the S-matrix ensembles 
labelled by /3 = 2 (these were introduced in references [7H1 [7^1) can be associated to a 
symmetric space. Let us briefly recall the problem. The S-matrix ensembles correspond to 
root lattices of -BC„ type with the following root multiplicities: rrio = (3, = (3 — 2, nii = 
1, where (3 takes values in the set {1,2,4}. Only the (3 = 2 case corresponds to a set of 
root multiplicities associated to a symmetric space. In this case the restricted root lattice 
degenerates into C^, since the multiplicity of short roots = 0. It is instructive to re- 
obtain this symmetric space description from the integration manifold of the random matrix 
ensemble. In the unitary case the scattering matrix is parametrized by (cf. eq. ()8.14|1 in 
paragraph I8.3.4|l 




(10.1) 
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where r, t, r', t' are N x N matrices. One is interested in the eigenvalues {Tj} of the product 
since these determine a variety of transport properties. This imphes that there is a 
hidden symmetry: 



S^S{-1) (10.2) 

with f3, f 4 G f/ (N). Indeed, under the transformation p().2j) . t tv^ and the product tt^ is 
invariant for unitary v^. This additional symmetry defines the space in which the matrices 
S live. It is not SU{2N) but the coset SU{2N)/S(U{N) x f/(A^)), which is exactly the 
symmetric space described by the Cjv root lattice with multiplicities = 2 and mi = 1 
mentioned above. 

We can write the coordinates {Tj} explicitly by using the fact that any 2N x 2N unitary 
matrix can be decomposed as [THIIZSI (cf. eq. ()8.15p ) 



S 




-Vi^l — T 

V2\pT V2VT 



;io.3) 



where f 1, V2, fs, V4 G U{N) (in the presence of time-reversal symmetry fi, fs and V2, are 
related to each other by transposition) and r is a x diagonal matrix which collects 
the coordinates {Tj}. As was obvious to begin with, they are not the radial coordinates 
corresponding to a Cartan subalgebra in the symmetric space of the scattering matrix (cf. 
paragraph I8.3.4|) . 

The same reasoning can be followed in the 13 = 1 and 4 cases. In these cases the cosets 
which one obtains by imposing the gauge symmetry discussed above turn out not to be 
symmetric spaces (paragraph 18.4. l() . 

We can look at these random matrix ensembles as plain circular ensembles in which, for 
well-defined physical reasons, we have chosen a set of parameters (the eigenvalues of the 
matrix tt'^) different from the eigenvalues of the matrix S (i.e., the radial coordinates of 
the symmetric spaces associated to the Dyson circular ensembles). 

Since the standard radial coordinates are related to the Cartan generators, we can consider 
this non-standard choice of "radial" coordinates as a non-Cartan parametrization of the 
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symmetric space of the scattering matrix. As we will see in a moment, the important 
point is that the construction of this non-Cartan parametrization follows the same lines as 
the standard parametrization (identification of an involutive automorphism, separation of 
the generators in even and odd, and so on). As a consequence of this, several properties of 
the standard Cartan parametrization are conserved, and allow us to treat these ensembles 
in a way which is essentially the same as for those which correspond to standard Cartan 
parametrization of symmetric spaces (see the results of [75117^] ). 

Below we will see in detail how this non-Cartan parametrization is constructed. To help 
the reader we will first briefly review the general method (a detailed discussion can be 
found in Appendix A of |H2I) and then look at an example. As we will see, the procedure 
has many similarities to the procedure we used in constructing the restricted root lattice 
in paragraph 15. 2[ and like in subsection 15.31 we will use two successive involutions. 



10.1.1 Non-Cartan parametrization of SU{N)/SO{N) 



Let us begin with a compact^^ symmetric space G/K and denote the corresponding algebra 
subspace G/K. Suppose a is the involutive automorphism that splits the algebra G into 
K © P. Let's now operate on the subspaces K and P = G/K with a second involution r 
exactly like in paragraph 15.31 Then the subspaces K and P are split into even and odd 
parts that were named Ki, K2, Pi, P2 in subsection 15.31 (In reference P2] these are 
the subspaces /Cg, )Co, 'Pe, 'Po-) Let A denote a maximal abelian subalgebra contained in 
P2. That means its elements are odd under the transformations a and r: aAa = —A, 
tAt = —A for A E A (i.e. the elements of A anticommute with the involutions). 

Let M be the subspace of elements of Ki that commute with all the A & A. The subspace 
of G that is invariant under r is Ki © Pi = Gi. The central ingredient in obtaining the 
non-Cartan parametrization is the bijective mapping (p of the manifold Gi/M x A'^ into 
the symmetric space G/K such that (p : {gM, a) —>■ gaK, where elements in the coset space 
are denoted gK (cf. paragraph 12. 3|) . The non-Cartan "radial coordinates" are encoded in 
the matrix a belonging to a connected open subset A'^ oi A = e^. In the tangent space 
this corresponds to a mapping between algebra subspaces [H2] 



Ta : Gi/M X A ^ G/K 

Ta{Z,H) = H + aZa-^\p (10.4) 
^°We thank Martin Zirnbauer for suggesting this possibihty to us. 

^^The fact that we choose a compact space for definiteness is not important and we can equally well 
choose a non-compact one. 
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where the notation |p means the restriction to the subspace P. From the hnear mapping 
Ta we will obtain the Jacobian as a determinant of the matrix expressing the differential of 
(j) with respect to the new basis. Keeping in mind the commutation relations for symmetric 
subalgebras (eq. ()3.2|) of subsection 13. ip and the fact that A anticommutes with the two 
involutions a and r, we easily see that 



[A, Pi] C K2, [A, Ki/M] C P2/A, [A, K2] C Pi, [A, P2/A] C Ki/M (10.5) 

This means that the mapping ad(lna) from the algebra G to itself maps the four eigensub- 
spaces of cxr into each other as follows 



ad(lna) 


Ki- 


^P2 


ad(lna) 


P2- 


-Ki 


ad(lna) 


K2- 


^Pi 


ad(lna) 


Pi- 


^K2 



(10.6) 



Denoting the adjoint action aZa^^ in fll0.4|l with Ad{a)Z and using Ad(a) = exp ad(lna) = 
coshad(lna) + sinhad(lna), we see from equation ()10.5|) that if Z = X + Y is the decom- 
position of Z G Gi/M into parts belonging to the subspaces Pi and Ki/M respectively, 

TaiZ,H) = + cosh ad(lna)X - sinhad(lna)F (10.7) 

(this follows from equations p0.6|) because coshx is an even and sinhx an odd function of 
X, and keeping in mind that in Ta we take the projection on P). 

The Jacobian corresponding to the change of coordinates is the determinant [H21 

JNc{ci) = det (coshad(lna)lpi^P^) det (^sinhad(lna)|Ki/M^P2/A) (10.8) 

and it is obtained as the product of the eigenvalues. The eigenvalues of the automorphism 
ad(lna) are nothing but the restricted roots with respect to the abelian algebra A. We 
therefore obtain the general formula 
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J^^(a) = Yl sinh'""a(lna) [] cosh""^ (3{\na) (10.9) 

where the subscript NC stands for "non-Cartan" and the index + reminds us of the fact 
that the sum is over the positive roots only. Q;(lna) is nothing else than the projection q"' 
introduced in equation ()6.27|) in subsection 16.21 The positive roots have been divided into 
two subsets Rf and R2 in an obvious notation. In case G/K is a compact space, the roots 
a and (3 in ()1().9|1 are purely imaginary. If we then set a = ia' and (3 = i(3' we obtain 

jW(a) = [] sin'^-'a'(lna) J] cos'^^'/3'(lna) (10.10) 

whereas for real a, P (non-compact space) we get hyperbolic functions in eq. ()10.9|) for 
J^(a) (cf. the similar situation in eq. ()6.43p ). For comparison, recall that the Jaco- 
bian corresponding to the standard Cartan parametrization of the G/K space is given by 
equation ()6.30p . 

j(+)(a)= W sin'""a(lna) (10.11) 

where the algebra corresponding to A in this case is a maximal abelian subgroup of the 
whole subspace P. 

Example: Let us now take as an example the (3 = 1 S-matrix ensemble, whose parametriza- 
tion in subsection 18.4.11 falls outside of the Cartan classification of symmetric spaces. We 
will identify it as a non-Cartan parametrization of a standard symmetric space. 

The starting point in this case is the circular orthogonal ensemble of Dyson, i.e. G/K = 
SU{N)/SO{N). Let us assume for completeness that we have a different number of left 
and right scattering channels p = Nl and q = iV^j. In this case the involution r is 
given by Ip^g = lNL,Nf. = diag(lAr^, — Ia^^) (cf. equation ()4.15|) ) with fixed point set Gi = 
SU {Nl) X SU (Nfi) and the rank of a maximal abelian subgroup A is r = mm{NL, Nn). 

We can choose A to be generated by the matrices 
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where Eij is the matrix having a 1 in the ith row and jth cokimn, and zeros elsewhere. 
After finding the corresponding roots hke in the example in subsection 15.21 the radial 
coordinates a' ■ x (leaving out the i) and the corresponding root multiplicities turn out to 
be 



Rf: {Xk±xi) (m = l) {k > I) 

Xk {m = \NL-NR\) (A; = l,...,r) 
: (xfc ± Xi) (m =1) (A; > /) 
2xk (m = 1) (A; = l,...,r) 
Xk {m=\NL-NR\) (A; = l,...,r) 

(note that the total root multiplicities mo = 2, mi = 1, = 2{Nl — Nr) are exactly the 
ones for a BC]^ type restricted root lattice; cf. Table H}. These are the radial coordinates 
in the Jacobian for the non-Cartan parametrization, eq. p0.10|) 



jj^(a) = sin(xi— Xj)sin(xj+Xj) sin^'xfc J]^ cos(xi— Xm)cos(xi+Xm,) cos(2x„) ]^ cos'^Xg 

i>i k l>m n q 

(10.12) 

where z/ = \Nl — Nr\. By making now the variable substitution = sin^2a;fc in this 
Jacobian, we obtain the radial measure in the form 

J^^k^)da = n |T. - T,| n |T,r--^«l-^)/2^T,, (10.13) 

i<j k 

which agrees with equation ()8.69|) for the special case A''^ = Nr and (3=1. As anticipated 
in paragraph l8.4.11 we conclude that equations (|8.69I8.70|) represent the Jacobian for a non- 
Cartan parametrization (albeit expressed in different variables) of a standard symmetric 
space, which for /3 = 1 is the space SU{N) / SO{N). 



10.2 Clustered solutions of the DMPK equation 

As we mentioned previously, there are two major drawbacks in the DMPK approach to 
quantum wires. The first is that the DMPK description only holds in the quasi-one 
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dimensional limit. The second is that the solution discussed in paragraphs 19.4. ll and 19.4.'^ 
does not allow studying the intermediate cross-over region between the metallic and the 
insulating regimes, where no simplifying approximation is allowed. This is true even in the 
simplest case, for /5 = 2. 

In the last few years various generalizations of the DMPK equation have been suggested 
[HT] to avoid the quasi-one dimensional limit. However, as we discussed at the beginning 
of subsection 19.11 in all these generalized equations most of the attractive properties of 
the DMPK equations are lost, mainly due to the fact that the description in terms of 
symmetric spaces is no longer valid. Only a few pieces of information on the expected joint 
probability density of the transmission eigenvalues can be obtained. 

Recently a different strategy has been proposed in [HH" where the DMPK equation is kept 
unchanged, but one looks for a set of special solutions (with non trivial initial condi- 
tions) which break the isotropy ansatz. To this end one uses the exact integrability of 
the Calogero-Sutherland models: recall that the DMPK equation can be mapped into the 
evolution operator of a suitably chosen Calogero-Sutherland model (for a review see P). 
This is a highly non-trivial property which is more general than the underlying symmetric- 
space structure; in fact it holds also for generic integer values of the root multiplicities jH] . 
It is possible to show that as a consequence of their exact integrability, in these models 
- besides the well known symmetric solution - a wide class of non-trivial (but exact) so- 
lutions exists, in which the particles are grouped into clusters. Once the mapping to the 
DMPK equation is performed, the clusters of particles become clusters of eigenvalues. The 
exact integrability of the Calogero-Sutherland model ensures that this asymmetric distri- 
bution of eigenvalues survives in the asymptotic limit, and the remarkable properties of the 
underlying symmetric space allow us to explicitly write down such asymptotic expansions. 

Let us discuss these solutions in more detail. We assume the cluster to be composed of the 
first N' < N eigenvalues. This means: 



|x,-x,|<oo, z,j = l,...,N' [i^]). (10.14) 

|x, ^ oo, i = l,...,N- J = N' + 1,...,N (t^j) 

In the symmetric space framework we can identify the cluster by selecting a subsystem 
of roots associated to the space. Let 11 be the system of simple roots associated to the 
symmetric space X, and 11' a subsystem of simple roots which satisfies the inequality 

W = {ae n| lim < oo} 
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(10.15) 



where = (x, a). At this point there are two possibihties. Since 11 is a Cat type lattice, 11' 
can be either of type Cn (in this case it must also contain the long root, and the ordinary 
roots must be chosen so as to preserve the Z2 symmetry of the lattice), or it can be of 
type An- In both cases one can construct, from the ordinary roots of 11', the differences 
a;° = — Xj+i which correspond to the nearest neighbor distances between eigenvalues. It 
follows from the definition of 11' that these distances must remain finite in the asymptotic 
limit, so that 11' defines a cluster if it is connected, or a set of clusters otherwise. If the 
cluster is of type A^, there is no other constraint and the cluster can in principle flow to an 
infinite distance from the origin (while the eigenvalues inside the cluster are kept at a finite 
distance from each other). On the contrary, if the cluster is of type Cn, the eigenvalues 
are bounded by the lattice structure of 11' itself and consequently remain within a finite 
distance from the origin. In the following we denote the radial coordinates outside the 
cluster by x and the ones inside the cluster by x'. 

The asymptotic expansion of the zonal spherical functions in the presence of such a cluster 
was obtained a few years ago by Olshanetsky in [HH]. It turns out to be a rather natural 
generalization of the Harish-Chandra asymptotic expansion, eq. 1)9.311) : 

E c,i^k)e'^'''''^^irkyix') (10.16) 

r£W/W' 

where W which appears in the coset W/W is the Weyl group associated to the cluster, 
and \l/fc(a;) = C,{x)^k{x) (cf. equation ()9.3ip ). The function Cz{k) is defined by 

c.(A;)= n Ca(A;) (10.17) 

a<^R+/R'+ 

where R'~^ is the set of positive roots associated to the cluster. In pO.Kij) (rfc)' denotes the 
projection of the vector rk on the sublattice 11' and rk is its complement. The function 
C,{x) is given by eq. ()7.9|) : 

^(^) = II I sinh^Xj — sinh^Xjl^ H | sinh2xj|2 (10.18) 

i<j i 

In spite of its apparent simplicity, the expression ()10.16|) is highly non-trivial. Notice 
for instance that the symmetrization with respect to the Weyl coset W/W acts not only 
on the part containing the coordinates x but also on the momenta of the zonal spherical 
function describing the cluster coordinates ^(^rkyi^')- This means that the particles inside 



143 



the cluster do not move independently in a section of the whole space but they "feel" the 
presence of the other particles and are subject to the symmetry group of the remaining 
space. We refer the reader to the recent paper [88 for some explicit examples of this type 
of solutions and for indications of how they could be used to address the two problems 
mentioned at the beginning of this subsection. 

We only mention here that the clustered solutions described in this section may find another 
natural application if one tries to model systems in which the number of open channels is 
reduced by the structure of the wire itself (cf. the wide-narrow-wide geometry of [SO!)- In 
this case, one could consider a configuration formed by a Ctv type cluster (bounded to the 
origin) made of N' eigenvalues and let the remaining N — N' eigenvalues flow to infinity. 
One can choose A^' and s (the dimensionless length of the wire) so as to keep the cluster 
in the metallic regime, while the other eigenvalues are in the insulating regime and do not 
contribute to the wire conductance. 



10.3 Triplicity of the Weierstrass potential 

The Calogero-Sutherland potentials of type IV in eq. ()7.2j) . 



where V{^) is the Weierstrass P-function, is the most general type of potential. The 
Weierstrass P-function is defined as 

where the prime on the sum means we are summing over all pairs (m, n) G except 
{m,n) = (0,0). The P-function is doubly periodic with periods 2uji, 2uj2, which can be 
seen by rearranging the sum. The Weierstrass "P-function is an elliptic function. An elliptic 
function is defined to be a single-valued, doubly-periodic analytic function, whose only 
singularities in the finite part of the complex plane are poles. We will show that as both, 
or just one, of the two periods goes to infinity, we recover the potentials of type I, II and 
III, respectively. This fact was mentioned, but not proved, in P], and a proof for the III 
type potential was briefly indicated in Appendix A of ^T]. Set 
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V{z) 



oo 

E 

m=— oo 



4sin2 (^±f^ 



10.211 



where z is a. complex variable. It is immediately obvious that V{z) is doubly periodic 
with periods 2tt and ij3, has a double pole in each period-parallelogram ai z = 2'Kk + i(3m 
{k, m G Z) and is analytic elsewhere. The function V{z) — \jz^ is analytic in z = and in 
a neighborhood around this point. Expanding it around z = we find 



z=0 



1 °° 1 

T2"S2sinh2 (^) 



(10.22) 



Let ti^fcm = 27rfc + z/3m denote the position of the double poles of V{z). Since V{z^ is 
periodic, analytic except in the poles, and has double poles in w^^rn it must be of the form 



C-2 



+ 



C-1 



Z - Wkm) 



^ Cn[Z - Wkm) 
n=l 



(10.23) 



From the condition that f{z) = V{z) — Ijz^ is analytic and equal to Z^^o (0)-^"/''^' 
and near ^ = we then get c_2 = 1, c_i = by matching terms. 

We now see that the two elliptic functions V{z^ and V{z\'k jT) have the same pe- 
riods, poles, and principal parts at each pole. They then differ by a constant k (see 
P^ . paragraph 13.11), so that c„ = for n > in (jl0.23p . So we have proved that 
V{z) -l/z^ + k = V{z; 7r,if3/2) - 1/z^ where k is defined by eq. (|1().22|1 . since both sides 
have to have the same value at 2; = 0: 















;v(.) - ^; 





{z - 2TTk - i(3my {2TTk + i/3my 



;i0.24) 



or in other words 



Viz; 71,1(3/2)= Yl 



m=—oo 



4sin^ 



/ z+i/3m \ 



12 + ^ 



m=l 



Isinh^ (^) 



(10.25) 
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It is now a simple matter to show that in the hmit /5 — > oo, 



4sm^ f Ij 

that is, apart from a constant ^, ^(-z; 7T,i(3/2) in the hmit j3 ^ oo becomes a potential of 
type III in eq. (TT^ . 

In a completely analogous way one shows that the hyperbolic potential 



V{z)= V -7 ^ (10.27) 

™£^oo4sinh2(^) ^ ' 



is related to the Weierstrass 'P-function by 



OQ -I 1 CX) 1 

P(z;a/2,i7r) = V —, ^ + V — ^ (10.28) 

' ^ 4sinh2(^) 12 i'i2sinh2(f^) ^ ^ 



m=— oo 



As a — > CX), the potential 'V[z) approaches 



^(^) ^ . . (10-29) 
4smh^ (^f j 



so that ^(z; al2,m) becomes a potential of type II in eq. ()7.2|1 . 
Also, if both periods go to infinity, we obtain a potential of type I: 



P(z;a/2,i/?/2) ^ ^ as a.p^oo (10.30) 



^2 



The important consequence of this analysis is that we can see the triplicity of the sym- 
metric spaces as a limiting procedure on the general Weierstrass potential. This gives us 
a framework for interpolating between spaces of different curvature (at fixed root lattice), 
and may have relevant applications in the context of random matrices. Indeed, in the 
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last few years a substantial amount of non-trivial mathematical results has been accu- 
mulated concerning these V{z) type models [HSl- The hope is that some of them will be 
useful for constructing non-trivial generalizations of the known random matrix theories, 
while preserving some of the attractive properties due to the underlying symmetric space 
structure. 



11 Summary and conclusion 



In this review we have discussed the usefulness of viewing random matrix ensembles as 
symmetric spaces. Random matrix theory, that has evolved into an important branch of 
mathematical physics, is used in the description of physical systems with chaotic behav- 
ior, disorder, or a large number of degrees of freedom. The versatility of random matrix 
theories allows for a parameter-free description of an assortment of systems, ranging in 
size from nuclei to mesoscopic conductors. The unifying feature of these systems is chaotic 
behavior effectively resulting in randomness, and the available information is the universal 
statistical behavior of spectra. On the other hand, symmetric spaces are well-understood 
mathematical objects that can be represented as coset spaces G/K of a Lie group G with 
respect to a symmetric subgroup K (or as Lie algebra subspaces). After a general intro- 
duction to Lie algebras and root spaces, we have seen how to construct symmetric spaces 
from semisimple Lie algebras using involutions, and how to identify them with the integra- 
tion manifolds used in random matrix theory. In the process we gave concrete examples of 
all the mathematical concepts that were introduced. We discussed coordinate systems on 
symmetric spaces and identified the spherical radial coordinates with the physical degrees 
of freedom in most matrix models (an exception is the transfer matrix ensembles where 
the physically interesting degrees of freedom, in the Landauer theory, are the transmission 
eigenvalues, which are not the same as the eigenvalues of the transfer matrix). We ex- 
plained how the Dyson and boundary indices of the random matrix ensemble are related 
to the multiplicities of the restricted roots associated to this symmetric space. 

We have seen that the metric on a Lie algebra, defined in terms of the adjoint representa- 
tion, leads to the concept of curvature tensor on the symmetric space, and we have shown 
that the symmetric spaces classfied by Cartan appear in triplets of positive, zero, and neg- 
ative curvature corresponding to a given (restricted) root lattice (and therefore to a given 
set of multiplicities of long, ordinary and short roots). Further we have discussed Casimir 
operators in a general context and their representations in local coordinates (the Laplace 
operators) and given a general formula for obtaining the radial part of the Laplace-Beltrami 
operator in terms of the Jacobian resulting from a transformation from random matrix to 
eigenvalue space. We devoted a subsection and the appendix to their eigenfunctions, the 
zonal spherical functions, that were later used in some of the applications. 
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In section [S] we gave a general introduction to some commonly used random matrix the- 
ories, explaining how universality of the correlation functions leads to general universal 
predictions concerning the statistics of the eigenvalues. The spectral properties that can 
be described in this way are exactly those that are independent of the detailed dynam- 
ics of the system. The only input is the global symmetries. The fact that many of the 
physical systems described by random matrices fall into the universality classes of Cartan's 
classification is a consequence of this fact. The eigenvalue density is determined by the 
matrix potential and by the Jacobian obtained in diagonalizing the random matrices on 
the symmetric manifold. The Jacobian is the origin of spectral correlations in the matrix 
models discussed here. It is completely determined by the restricted root system of the 
underlying symmetric space and by its curvature. 

In this review we also discussed the mapping of Calogero-Sutherland models onto (re- 
stricted) root systems of Lie algebras and symmetric spaces. The mapping is based on 
the fact that the Hamiltonian of these models for certain values of the coupling con- 
stants (determined by the root multiplicities) can be exactly transformed into the radial 
Laplace-Beltrami operator. This mapping allows to obtain several exact results for the 
zonal spherical functions of the corresponding symmmetric spaces (see the discussion in 
the appendix). 

The Dorokhov-Mello-Pereyra-Kumar equation is the differential equation determining the 
joint probability distribution of the transfer matrix as a function of an external parameter. 
The operator appearing in this equation is mapped (by a transformation involving the 
Jacobian) onto the radial part of the Laplace-Beltrami operator representing free diffusion 
on the symmetric space underlying the random matrix ensemble. This can be used to 
solve the Fokker-Planck equation exactly for Dyson index (3 = 2, and to obtain approx- 
imate solutions for the probability distribution and the conductance in other cases using 
asymptotic expansions of the zonal spherical functions. As further examples of applications 
we reviewed how the connection from random matrices to symmetric spaces can explain 
a discrepancy between random matrix theory and the non-linear sigma model regarding 
the magnetoconductance, and further how the orthogonal polynomials and the symmetries 
of the matrix models can be traced directly to the root lattice. In the last section we 
discussed some applications going beyond the Cartan classification, and we showed that 
the Weierstrass potential of Calogero-Sutherland models in its various limits reflects the 
possible signs of the curvature of the symmetric space. 

An interesting project that might be worth exploring in the same spirit as above concerns 
non-hermitean matrix models, i.e. random matrices having complex eigenvalues. Recently 
there has been a boost of activity in this field. Since non-hermitean random matrices are 
useful in a number of contexts, but not much is known about them (see however [94 , 9511^ ). 
a thorough study of them would be extremely important. A classification according to their 
symmetries was attempted in ^H], where 43 spaces were enumerated. However, a study 



148 



of the resulting spaces in the spirit of this work has not yet been performed. Some non- 
hermitean ensembles have only recently been applied to physical systems. 



Exploring what is known about the corresponding manifolds, and how this knowledge could 
be applied in the above mentioned contexts, would appear to be a worthwhile effort. We are 
not aware of the extent to which the properties of such manifolds are known in mathematics, 
and how useful this knowledge might be in the problems involving non-hermitean random 
matrix theories. Let us mention a few physical problems in which non-hermitean random 
matrices are present: 

a) Non-hermitean matrices are important in schematic random matrix models of the QCD 
vacuum at non-zero temperature and/or large baryon density (i.e. finite chemical 
potential). The spontaneous breaking of chiral symmetry is one of the most important 
dynamical properties of QCD, as it shapes the hadronic spectrum. The study of the 
chirally symmetric phase of QCD is one of the primary objectives of heavy ion colliders. 
One of the main problems with including a chemical potential in the QCD action is that the 
fermion determinant becomes complex, which makes lattice simulations extremely difficult. 
By using schematic random matrix models, one can study qualitative features of QCD at 
finite density and/or temperature. 

b) In classical statistical mechanics, non-equilibrium processes can be studied as the time- 
evolution of non-hermitean Hamiltonians. As an example we mention non-hermitean spin 
chains related to a Kardar-Parisi-Zhang equation. Such an equation can describe for 
example interface growth, problems in fiuid dynamics, a driven lattice gas, or directed 
polymers or quantum particles in a random environment j98||. 




c) Non-hermitean effective Hamiltonians appear in the S-matrix description of open sys- 
tems connected to reservoirs. The S-matrix characterizes scattering in an open chaotic 
system like for example a ballistic microstructure pierced by a magnetic fiux j99 j . 

d) Neural networks are described by continuous local variables related through a non- 
linear gain function to a local field that could represent the membrane potential of a 
nerve cell. The dynamics of the network is described by a large number of coupled first 
order differential equations featuring (random) non-hermitean matrices. These contain 
the parameters coupling the output of the j'th neuron to the input of the i'th neuron. 
In the context of neurobiology, the study of chaos in neural networks is relevant to the 
understanding of several features of neural assemblies |10Uj . 

e) Anderson localization in a conductor denotes the phenomenon when the wave function 
of an electron becomes localized and exponentially decaying. As a result the conductor 
becomes an insulator. In conventional conductors localization is known to happen in di- 
mensions d > 2 [d = 2 is the critical dimension). New interesting phenomena may occur in 
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systems with non-hermitean quantum mechanics. In |lUlj . the authors considered parti- 
cles described by a random Schrodinger equation with an imaginary vector potential. The 
model was motivated by flux-line pinning in superconductors. They found that delocaliza- 
tion transitions arise in both one and two dimensions. Such models were discussed further 

in [ma rina.. io4j . 

Even if the suggested direction of research concerning non-hermitean matrix models turns 
out not to be feasible, there are likely to be many more applications for hermitean matrix 
models that can be pursued within the framework of symmetric spaces. 
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A Appendix: Zonal spherical functions 



A remarkable feature of the theory of zonal spherical functions (ZSF) discussed in subsec- 
tion IHISl is that in some special cases exact formulae exist for these functions. This happens 
in particular for all the symmetric spaces with multiplicity of ordinary roots rrio = 2 in 
the classification of Table ^ or Table El The ultimate reason for this becomes clear by 
inspection of the mapping discussed in section [3 The radial part of the Laplace-Beltrami 
operator can be transformed into a suitable Hamiltonian of n interacting particles which 
decouple exactly for rrio = 2. Thus we may expect to be able to write the ZSF's in these 
particular suitable determinant of single-particle eigenfunctions. 

In the simplest case in which the multiplicity of both the short and long roots is zero one 
can write 



(/.,oc[e(gr^Edetse^(^'''') (A.l) 

sew 

where W denotes the Weyl group (which in this case is simply the permutation group, 
while det s is simply the sign of the permutation s), X labels the ZSF and is related to its 
eigenvalue (see eq. ()6.45p ). ^(g) is given by eq. ()7.9|) and (A, q) denotes the vector product. 
In the above expression we recognize, as anticipated, the determinant det e'^'''^ written as 
an alternating sum over the permutation group W = 5'„. Notice, as a side remark, that if 
the curvature of the symmetric space is positive, then A is "quantized" and must belong 
to the dual weight lattice. 

This expression fixes the g-dependence of the ZSF. However, it does not completely fix the 
A-dependence, which is hidden in the proportionality constant. Fixing this constant turns 
out to be a highly non-trivial task. In the case of the An type spaces it can be obtained by 
direct integration of eq. (j6.44j) jH] . The final result (in which we have introduced explicitly 
the parameter a as in paragraph 16.31 so as to describe the three possible curvatures in one 
single equation) is 



ni<fc Uj<k smh[a{qj - g^jj 

The extension to other root lattices is slightly more involved. Due to the presence of the 
long and/or short roots or to the additional reflection symmetry of the lattice, the corre- 
sponding "single particle" Schrodinger equations are more complicated and the solutions 
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are no longer simple plane waves. The problem was solved in 1958 by Berezin and Karpele- 
vick [HS]. The solution they obtained is valid for a generic BCn lattice and, remarkably 
enough, holds for any value of and m/, not only for those related to the symmetric 
spaces (while obviously nio = 2 is mandatory). 



In this case there is no compact way to deal with the three possible curvatures in a single 
formula, so let us look at the three cases separately: 

Positive curvature spaces 



det 



pa./3('cos2(7' 



sin^ ql) 



(A.3) 



where the -P^'^ are Jacobi polynomials with a = {rris + rrii + 1) [3 = {mi — l)/2 and 
7)1 = (2Xi — Tils — 2mi)/4: (recall that in this case A is quantized, thus rji turns out to be an 
integer number). The normalization constant C does not contain any further dependence 
on A and q. Its value is C = 2'^("-i) U?=i ^! (^. + O''"'- 

Negative curvature spaces 



det 



F{ai,bi,c; -sinh^g^; 



Ilj<ki^j 



Ai) ni<fc(sinh g2 -sinh ql) 



(A.4) 



where this time F(ai, 6^, c; — sinh^ g^) denotes the hypergeometric function and the three 
parameters are : ai = (ms + 2mi + 2zA;)/4, bi = (ms + 2m; — 22A;)/4 and c = (ms + m; + l)/2. 
The normalization constant is the same as in the positive curvature case. This is the formula 
that we used in paragraph 19.4. II to solve the DMPK equation in the {3 = 2 case. 

Zero curvature spaces 

In this case the expression is slightly simpler: 



Co 



det [{Xiqj) ^J-y{Xiqj)] 



(A.5) 



where 7 = (m^ + m/ — l)/2 and 
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Co 



(_l)"(n-l)/2 2n(n+m-3/2) + 1/2) ... T{m + 71-1/2) 

l!2!3!...(n-l)! 



(A.6) 



with rh = {nis + m;)/2. 



A.l The Itzykson—Zuber—Harish— Chandra integral 



A remarkable feature of the zero curvature case discussed above is that not only the zonal 
spherical functions themselves are simpler, but also the integral representation from which 
they are obtained drastically simplifies. In fact, in this case (recall that rrio = 2) the 
symmetric spaces coincide with the classical Lie algebras and eq. ()6.44|) can be written as 
a simple integral over the group manifold. For instance in the unitary case, i.e. for systems 
of type An, the integral representation becomes 



0f (g)=|Dt/e^'^(^^«^') (A.7) 

where U G SU{n), A and Q are diagonal matrices: A = diag(Ai, A„) and Q = 
diag(gi, 

On the other hand, if we inspect the a —>■ limit of eq. ()A.2|) we find: 



(0) . . ^ l!2!...(n-l)!Esgi^det.e-(-^-'^) 



In this last expression we recognize the det e'^** written as an alternating sum over the 
symmetric group W = S'„, while the two products in the denominator can be interpreted 
as Vandermonde determinants 



A(A) = n(A^ - A.) (A.9) 

i<j 

Moreover, we easily see that eq. ()A.7|) holds unchanged if A and Q are two generic her- 
mitean matrices with eigenvalues {Aj} and {qi}, respectively. At this point, equating these 
two expressions for 0^°-' we find the remarkable relation: 
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which is the well-known Itzykson-Zuber-Harish-Chandra (IZHC in the following) integral. 
This result was originally obtained by Harish-Chandra in |lU5j . It was later rediscovered 
in the context of random matrix models by Itzykson and Zuber |l(J6j and fully exploited 
by Mehta in |107j . It plays a major role in several physical applications of random ma- 
trix ensembles, ranging from the exact solution of two-dimensional QCD to the study of 
parametric correlations in random matrix theories. 

It is not difficult, using equation ()A.5|) . to generalize eq. ()A.10|) to other Lie groups, i.e. 
to other root lattices. This general result was already present in the original paper by 
Harish-Chandra |lU5j . A recent review with detailed expression for SO{n) (both for even 
and odd n) and Sp{n) can be found for instance in appendix A of |lU8j . 



A. 2 The Duistermaat— Heckman theorem 

The remarkable elegance and simplicity of the IZHC integral suggests that there should be 
some deep mathematical principle underlying this result. Indeed, in the last few years it 
has been realized that the IZHC integral is a particular example of a wide class of integrals 
which may be solved exactly with the saddle point method provided one sums over all the 
critical points (and not only over the maxima). In the IZHC case this means that one has 
to sum over all the elements of the Weyl group (i.e over all the permutations in the unitary 
case) and not only over the ones with positive signature (as one would do in a standard 
saddle point approximation). 

The rationale behind this remarkable result is the Duistermaat-Heckman theorem |109j . 
It states that the saddle point is exact (provided one sums over all the critical points) if 
the integration is performed over an orbit with a symplectic structure. As a matter of 
fact, it was later recognized that this is only an instance of a wider class of results which 
are known as localization theorems. Their common feature is that they can be used to 
reduce integrals over suitably chosen manifolds to sums over sets of critical points. A 
more detailed dicussion of this beautiful branch of modern mathematics goes beyond the 
scopes of the present review. To the interested reader we suggest reference |110j , where a 
comprehensive and readable review on the localization formulae and some their physical 
applications can be found. 
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